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SUMMARY 

A mathematical model and associated computer program (flSbas) to simulate a twin-tailed high 
performance fighter airplane (McDonnell Douglas F/A-18) are described. The simulation is 
intended to siroptwt advanced control law research fw application to high peiformance aircraft. 

The modeled flight envelope is extensive, allowing investigations in the high-angle-of-attack 
portion of the flight re^me. Included in die simulation mathematical model are the nonlinear, six 
depee-of-freedom rigid-body equations, an engine model, and simple sensor and actuator models. 
A simplified form of the F/A-18 digital control laws (version 8.3.3) is implemented, including only 
the Auto FUq) Up (AFU) flight mode used for air combat maneuvering. Aerodynamic forces and 
moments arc calculated from a wind-tunnel-derived database using table look-ups with linear 
interpolation. The aerodynamic database has an angle-of-attack range of - 10® to +90° and a sideslip 
range of -20° to +20°. Elastic deformation effects are incorporated in a quasi-static-elastic manner, 
with no dynamic simulation of the elastic degrees of fieedom. 
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1. INTRODUCTION 


The mathematical model and associated computer program (fISbas) to simulate a twin-tailed high 
perfonnance fighter airplane (McDonnell Douglas F/A-18) arc described. The simulation is 
intended to simport advanced control law research for triplication to high perfonnance aircraft 
The modeled rhght envelope is extensive, allowing investigations in the high-angle-of-attack 
portion of the flight repme. Included in the simulation mathematical model are the nonlinear, six 
de^ce-of-fneedom rigid-body equations, an engine model, and simple sensor and actuator models. 
A simplified form of the F/A-18 digital control laws (version 8.3.3) is implemented, including only 
the Auto Flap Up (AFU) flight mode used for air combat maneuvering. Aerodynamic forces and 
moments are calculated from a wind-tunnel-<terived database using taWe look-ups with linear 
interpolation. The aerodynamic database has an angle-of-attack range of -10® to -^90° and a sideslip 
range of -20® to +20®. Elastic deformation effects are incorporated in a quasi-static-elastic manner, 
with no dynamic simulation of the elastic degrees of freedom. In the engine model, the throttle- 
command^ steady-state thrust level and the dynamic response characteristics of the engine are 
based on airflow rate as determined from a table look-up. Afterburner dynamics arc switched in at 
a threshold based on the engine airflow and commanded thrust. The control surfaces are driven by 
simple transfer functions with constant rate and position limiting. 

1. 1 Comments on ACSL 

The simulation program written in the Advanced Continuous Simulation Language (ACSL). The 
ACSL simulation system, described in ACSL (1987), consists of a special purpose high level 
language, a uanslator, and various libraries to satisfy the commands available in the language. 
ACSL functions according to specifications established under the auspices of Simulation Councils, 
Inc. in 1967 [Format Subcommittee (1967)]. In a typical application of the ACSL system, a 
simulation model is ^t constructed using the special ACSL language. A u^slator is invoked that 
converts a model written in ACSL to Fortran 77 code. Fortran subroutines can be called from the 
smulation model with certain restrictions. The Fortran code is compiled and linked with the ACSL 
libraries. The resulting executable program will prompt for user input and produce plots and 
printed output. 

1.2 Simulation Development History 

In 1986, a team was established at the NASA Langley Research Center to address some of the 
conuiol-related problems of designing highly integrated aircraft. This team was composed of 
members of the Aircraft Guidance and Control Branch and the Aeroservoelasticity Branch at LaRC 
and was dubbed the Functional Integration Technology (FIT) team. The focus of the FIT team's 
effwt was highly maneuverable aircraft with minimal structural weight. It was anticipated that 
such aircr^t would require active strucUiral control. The original focus vehicle was the X-29 
research airplane. About 6 months after the FIT team was formed, the focus vehicle was changed 
to the F/A-18. Focusing on an F/A-18 allowed the FIT team to leverage off NASA's High Angle- 
of-Attack Research Vehicle (HARV) program. Because of the HARV program, NASA would 
own a preiHtxluction F/A-18 aircraft that would be used in an extensive flight test program. The 
X-29, on the other hand, was a joint program over which NASA had less control. 

One of the first tasks of the FIT team was to develop a simulation of an F/A-18 that incorporated; 

(1) nonlinear, six de^e-of-freedom, rigid body dynamics, (2) linear elastic modes, and (3) 
unsteady aerodynamics. At that time the team had access to a batch simulation computer program 
(nutp 8) and documents written by the McDonnell Aircraft Company of the McDonnell Douglas 
Corporation (MDC) for the F/A-18. The mdcflS simulation was a complete, nonlinear, six 
degree-of-fiteedom batch simulation that ran on a Digital Equipment Corporation VAX computer. 
The mdcflS simulation included the 8.3.3 flight control system, hinge-moment calculations, and a 
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full aerodynamic database. A considerable portion of the mdcflS simulation was devoted to failure 
modes of the propulsion and actuation systems. The FIT team also had access to a real-time 
simulation computer program of an F/A-18 that executed on a CDC Cyber 175 and was written by 
Unisys personnel for use with the Differential Maneuvering Simulator (DMS) at LaRC. The 
LaRC-DMS real-time simulation {dmsflS) was derived from the MDC simulation and documents 
provided by MDC. The FIT team chose to focus on dynamics issues rather t|iM failure modes. 
'There was also a retjuirement to interface closely with the LaRC real-time facilities and to begin 
setting up procedures whereby the transfer of specifications and code between researchers and the 
real-time simulation software development staff could be streamlined. The original FIT F/A-18 
simulation [Arbuckle, et al. (1987) and Buttrill, et al. (1987)1 was a batch simulation written in 
AC!SL that included 6 nonlinear degrees-of-freedom and 20 elastic degrees-of-freedom, but no 
control system. The rigid-body aerodynamic model and the engine model of the original FIT F/A- 
18 simulation were largely derived from the LaRC DMS real-time simulation {dmsfIS). 

Another research focus at LaRC is agility at high angle-of-attack and low-to-mo^rate dynamic 
pressures. To support researchers working in this area, the fully elastic simulation developed by 
the FIT team was modified. The elastic dynamics were removed. The code that simulated the 

8.3.3 flight control software in the LaRC-DMS real-time simulation was obtained and interfaced 
with the ACSL batch simulation. This report describes the resulting simulation code (fISbas). 

The initial simulation code (flHhas) developed to support the Aircraft Guidance and Control 
Branch 's high-angle-of-attack research was frozen so that this document could be written. The 
flSbas simulation models a preliminary-design version of the thrust-vectoring system proposed for 
the HARV program. The preliminary-design thrust vectoring system employed two paddles per 
engine nozzle, whereas the actual system subsequently installed on the HARV demonstrator has 
three paddles per nozzle. The thrust vectoring capability of \hcfl8bas simulation will only be 
described in a limited fashion as required to explain parameter lists. The simulation model without 
thrust-vectoring can be recovered by setting the logical, LTHVEC, to false. 

One outgrowth of the flSbas simulation was the flSharv simulation. The flSharv simulation 
represents the Aircraft Guidance and Control Branch's best current model of the HARV 
demonstrator. The mathematical model for flSharv simulation is documented in a variety of 
internal memoranda and will be more formally described in a future report. 

1.3 Code Development Philosophy 

ACSL allows one to call Fortran subroutines from within ACSL source code. Since ACSL source 
code gets translated into Fortran before compilation, this is not surprising. If one chose, the 
calculation of all the state derivatives could be done in a Fortran subroutine. The derivatives could 
be returned as outputs from the subroutine and integrated by ACSL. ACSL could plot and record 
variables in the parameter list. The equations of motion would then be a black box. Linear 
analysis on the ACSL-integrated states could be performed. Alternatively, if (me had an existing 
Fortran simulation, the entire Fortran simulation could be made into a subroutine and treated as, 
from the perspective of ACSL, a discrete process that returned outputs at current time, t, based on 
inputs at a previous time, t-h. The states of the Fortran subroutine would be hidden from ACSL 
and none of the linear analysis capability of ACSL could be used. 

The goals adopted by the authors for building the simulation were: (1) all integrations arising from 
continuous processes would result in ACSL-defined states, and (2) the equations of motion, to the 
extent possible, would be coded in the ACSL language. Goal (2) was compromised slighdy in that 
the calculation of aerodynamic forces and moments and engine steady state performan(» are done 
in Fortran subroutines. The authors had available software that greatly facilitated creating code to 
implement large table look-up type databases characteristic of aerodynamic force and moment 
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build-ups and the ACSL c^abilities in this area are limited. Since not everyone has ACSL, coding 
the aer^ynamic and engine static models in Fortran served to make major portions of the 
simulation model more accessible to the research community. The continuous dynamics of the 
simulation include actuator and sensor dynamics, en^e dynamics, and the rigid-body equations 
of motion. Riysically discrete dyn^ics, such as a digital flight control system, are not coded in 
ACSL. These dynamics are cod^ in Fortran subroutines and accessed by the simulation through 
subroutine calls which occur at discrete time increments using the DISCRETE block feature of 
ACSL. DISCRETE blocks are explained in the next section (Section 2). 
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2. SIMULATION STRUCTURE 

The general structure of an ACSL simulation is described in section 3 of the ACSL reference 
manual [ACSL (1987)]. Section 2 of this report summarizes key elements of section 3 of the 
ACSL manual and will cast the specific structure of this simulation in the context of the discussion 
and figures of the ACSL reference manual. 

ACSL imposes some limitations and programming rigor on its users. ACSL works best when the 
state derivatives arc explicitly defined. If absolutely required, the IMPLICIT operator provi^ by 
ACSL can be invoked and the resulting simulation will attempt to solve numerically the implicit 
formulation, but this would result in unacceptable computing performance in all but the simplest 
models. Implicit loops arc best solved by inserting fast filters or solving them algebraically. 

ACSL creates a large common block (ZZCOM) that contains all the simulation variables. Included 
in ZZCOM arc three arrays: the state array, the state derivative array, and the state initial condition 
array. A state is the output of an integration. A state variable in the state array cannot also be in the 
state derivative array. 

2.1 Overview 

The ACSL code structure shown in figure 2.1 is somewhat generic yet also contains the salient 
features of the FI 8 simulation. The DYNAMIC...END block of figure 2.1 results in the main 
simulation loop. The main simulation loop is indicated by the lower shaded box in the overview 
flow diagram of figure 2.2. In the main loop, if the termination conation, T > TSTP, is not met, 
the program will integrate forward in time over one communication interval, CINT {ACSL 
0987)]. The communication interval determines how often the code in the dynamic block that is 
not part of cither the derivative block or the discrete blocks gets executed. The default value for 
CINT in tYitflSbas simulation is 0.025 (1/40) seconds. 

The maximum integration step for the DERIVATIVE block is defined by MXSTP in \hcfl8bas 
simulation and has a default value of 0.00625 (1/160) seconds. The integration algorithm is 
selected by setting lALG. A variety of integration schemes are offered (lALG = 1 through 7) and 
arc described in ACSL (1987). The default choice lorflSbas simulation is IALG=4, a Runge- 
Kutta second-order method. The Runge-Kutta second-order method was recommended by the 
author of ACSL for most mechanical systems and has been found to work satisfactorily. No 
systematic investigation of integration methods has been made for this simulation. Using the 
Runge-Kutta second-order method, two derivative evaluations arc made to advance time over one 
integration step, MXSTP. 

The DISCRETE blocks arc executed according to their INTERVAL or SCHEDULE statements. 
DISCRETE block namel in figure 2.1 is an example of an interval block that gets executed every 
Tname seconds. DISCRETE block name2 in figure 2.1 is executed whenever the logical condition 
defined in its SCHEDULE command becomes true. This type of block might be used to reverse 
the velocity of a ball bouncing on a table. In ^tflSbas simulation, the digital control law 
DISCRETE block, FCS, is an INTERVAL block and is called every TFCS seconds. The default 
value for TFCS is 0.025 seconds. The choice of a value for TFCS is discussed in section 9 of this 
repOTt. When different blocks are scheduled for execution at the same time on the event list, such 
as when a DISCRETE block interval is an even multiple of MXSTP or some other DISCRETE 
block interval, the wder of execution is determined by order of appearance in the m^l definition 
code It is not required that the DISCRETE interval period be a multiple of the maximum 
integration size for the DERIVATIVE block. An INTERVAL DISCROT block sets up a 
time barrier on an event list. ACSL will shorten the last integration step to insure all simulation 
variables arc synchronized at the correct interrupt time. 
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PRCX}RAM 

INITIAL 

Includes siaiemcnts to be executed at lime zero after the nintime command "START" has 
been entered. Values in the the state IC array have not yet been transfened to the state 
anay. 

END 

DYNAMIC 

DERIVATIVE 

Includes the statements that dcHnc the continuous system dynamic model. The 
elements of the state derivative array, the variables to be integrated, arc denned here. 
A stale is the output of an integrator. 

SCHEDULE na/nc2 .XP.wir2 
END 

DISCRETE namel 

INTERVAL Tname = number 

An INTERVAL scheduled DISCRETE includes .statements that arc executed every 
Tname seconds. 

END 

DISCRETE namel 

Includes statements to be executed when the scheduling conditions calculated in the 
DERIVATIVE bkx;k (in this case) for DISCRETE /ia/ne2 become .true., .i.e. when 
variable var2 cros.scs zero in the positive direction. 

END 

The remainder of the DYNAMIC block includes statements that define auxiliary variables 
that are required for output only. Thc.se .statements arc executed every CINT seconds. 
CINT is the default name for the communication interval. 

TERMT( T.GE.TSTP) 

END 

TERMINAL 

Includes statements to be executed when the input logical variable to the MACRO 
TERMT becomes .TRUE. 

END 

END 


Figure 2. 1 . Outline of Typical ACSL Program 

The paths in figure 2.2 associated with the variables LTR, ITFLG, LWRTR and LDBTR arise 
from code installed to trim the simulation. This path does not arise from any of the typical ACSL 
constructs found in figure 2. 1 . The trim logic is discussed in more detail in section 2.2 and figure 
2.4 of this report 
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1TFLG = 1 


Call subroulinc LOAD. LOAD reads in trim 
input fik, ”F18TR”, and calls subroutines that 
initialize aero and engine data arrays stored in 
common blocks. 


(1) Set acuialor stale ICs to conliol surface 
position commands determined by the trim search. 

(2) Set engine state IC's to commands from trim search. 
Insures that engines are in steady state at l?=0. 

(3) Compute steady-stale sensor oulpuU. 


lNmAL...END 

codeblock 


rTFLG=2 


Call F18FCS to execute 
control law subroutine 
and gencralc actuator 
commands from pilot 
inputs and sensor data. 


Initial i/c integration routine by transferring the 
sUle IDC's to the states and calling the 
derivative evaluation code. 


Execute code in 
DYNAM1C...END block 
not contained in 
DERIVATIVE or 
DISCRETE blocks. 


Call TRIM subroutine to generate 
next set of trial trim ICs. Return 
ITFLG to indicate: 

(2) trim attempt ongoing 
(>2) trim ended. 


Execute various output options 
according to values of ITFLG, 
LVDEF, LWRTR and LDBTR 


Execute 

DERIVATIVE ...END 
code block. 


Integrate over the 
communication 
interval, Cl NT. Call 
DERIVATIVE and 
DISCRETE code 
blocks as required. 


Write variables on 
the OUTPUT and 
PREPAR lists to 
units DIS and PRN. 


Execute 

TERMINAL-END 

block 


Write variables on the 
OUTPUT and PREPAR lists 
to units DIS and PRN. 


Figure 2.2. Overview of Simulation Flow 
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Rrom execuiivc , 
ACSL> START < 



Call subroutine I,OAD. IX)AO 
reads in trim input file, ”F18TR", 
and calls subroutines that initialize 
aero and engine data arrays stored 
in cemmon blocks. 



(1) Set actuator state ICs to control 
surface position commands 
dcicrmined by the trim search. 

(2) Set engine state ICs to commands 
from trim search. Insures 
cftgiftes in steady stale M U=0. 

(3) Ccmpule steady-stale sensor 
outputs. 




Read in file ’’F18TR" to obtain 
state values for previously 
calculMed trim solutioiv 





Read in aero data and basic 
engine data arrays. Call 
TYIIiNG to initialize thrust 
vectoring engine daU arrays. 



Call CNSTRN to calculate ICs for {u,v,w) based 
on ( a, ^ V, p, n/) . Kinematics for the 
following 3 trim cases CrCASn=1.2,3) arc 
calculated; 

(1) wings level, 1 g flight 

(2) coordinated turn, p= MIJDGTR (bank angle) 

(3) pull-up, nz=GLOAD (load factor) 


Call FI Hies to execute control 
law subroutine and generate 
actuator commands from pilot 
itipuu and senior data. 


Initiali/e tiH^aUan routine by 
transferring the state IDCs to 
the states and calling the 
derivative evaluation code. 


Execute 

DERIVATIVE...END 
code block. 


Call ITUM lubroutirK to 
generate next set of trial 
trim ICs. Return FFPLO to 
indicate; 

(2) Trim attempt ongoing 
(>2) Trim ended. 


lixecule various output options 
according to values of I'lTl^G, 
LVDIT, LWRTR and LDBTR 


Execute ccxic in 
DYNAMJC...END block not 
contairKd in Df^tlVA'riVlt 
or DISCRimi block.. 


I.TR.and. 

(rm.G^2) 


T^TSTP 


Integrate over the 
communication 
interval, CINT. Call 
DliRIVA'nVIJ»id 
DISCRirrEcodc 
bkKks as required. 


Write variables on the 

OU'mn'md 

PKEPAR lists to units 
DIS and PRN. 


Evaluate Discrete 
Code Blocks 




DISCRETE FCS 


DISCRETE FDT 


DISCRIHR ABR 


Engine ( 
AR I ( r 


DISCRETE ABI. 


Execute 

13*RMINAE...liND 

block 


Write OUTPUT «id 
PREPAR list variables 
louniu DIS and PRN. 


Return to ACSL 

executive: 

ACSL> 



7 












In figure 2.3, two sections of the simulation flow shown in figure 2.2 are expanded. The actions 
pioduced by calling subroutine LOAD in the IN1TIAL...END code block are detailed and the four 
discrete blocks in the simulation are shown. The discrete blocks that model the flight control 
system (PCS) and provide a mechanism for making comparisons with flight data (FDT) are 
invoked by INTERVAL commands. The discrete blocks associated with the engine afterburner 
models (ABR & ABL) are invoked by SCHEDULE statements. 

2.2 Finding a Trim Solution 

ACSL has a TRIM command that works with purely continuous models. To use the ACSL TRIM 
command, a subset of the simulation states (outputs of integrations) is selected using the command 

ANALYZ CLEAR',’RELEAS'={ state list). 

The ACSL trim function works by adjusting the selected states until their state derivatives are zero. 
A combination a Newton step and a steepest descent step is used (section 5.2.4 of ACSL (1987) 
to adjust the states. Without some modification, this is not what one generally wants in an airplane 
trim problem. For example, a simple longitudinal level-flight trim problem might be: 


Find, 
such that 


|u) = |5e,T, a,0 r 


( 2 . 1 ) 


{y)=[d, V, q. H|’ = {0} (2.2) 

for fixed V and H, where 8 g, T, H, a, 0, V, q are elevator deflection, thrust, altitude, 
angle-of-attack, pitch angle, total airspeed and pitch rate, respectively. 


In this case, elevator deflection and thrust are control inputs, and altitude, angle-of-attack, pitch 
angle, total airspeed and pitch rate are states. Trying to find ot, V, q and H such that {y)={0) for 
fixed, arbitrary, values of 65 , T and 0 is almost always fruitless. To use the ACSL trim capability 
to solve the trim problem described in equations ( 2 . 1 ) and ( 2 . 2 ), one could add the following 
equations to the simulation model: 

5e = ki q (2.3) 

t = k 2 V -t- k 3 H (2.4) 


When attempting to trim, kj, k 2 and k 3 would be set to something nonzero. The states {5g, T, a, 
6 , q) can then be released and the ACSL command "ANALYZ TRIM" invoked. ACSL will then 

drive { fig, T, 6 t, 0, q ) to zero and in the pnKess, by virtue of equations (2.3) and (2.4), assure 

that V, q, H are zero. Equations (2.3) and (2.4) mimic the idea of a trim button or elevator trim- 
wheel. CJeneralizing this procedure to handle additional trim cases, such as trimming in a steady 
level turn or in a symmetric pull-up, could get cumbersome. In addition, it is not clear whether the 
ACSL trim feature could be made to work with a discrete-block control law in the loop. For these 
reasons, the authors chose to implement the alternative trimming strategy shown in figure 2.4. The 
trim strategy employed in flSbas is not without drawbacks, and a strong case can be made for 
retaining the ACSL structure and search algorithms as much as possible. 
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From executive 
ACSL> START 


ITFLG - 1 


Set first pass Rag 
for trim loop. 


Call subroutine LOAD. LOAD 
reads in trim Input file, T18TR", 
and csdts subroutines that 
initialize aero and er^gine data 
arrays stored in common blocks. 


INITIAL..END 
code block 


Trim Logic 


(1) Set actuator state IC's to control 
surface position commands 
determined by the trim search. 

(2) Set engine state ICs to commands 
from trim search. Insures 
engines in steady state at t»0. 

(3) Compute steady-state sensor outputs 



Call F18FCS to execute control 
law subroutine ar>d generate 
actuator commands from pilot 
inputs ar>d sensor data. 


Trim inner loop that allows - 
implicit loops k> settle out. 
Eg; y-f(xy) 

Implicit loops are aeated 
when LFCS is engaged. 


Call TRIM subroutine to 
ger>erate next set of trial C*s, 
Return ITFLG to indicate; 

(2) trim attempt ongoirtg 

(3) trim successful 

(4) trim failed 


ITFLG-4 


Initialize Integration routine by transferring 
the state IC's to the states and calling the 
derivative evaluation code. 


ITR< ITRMX 


Call DEBUG. 
Sal TSTP « 0. 


LWRTR 


f Write 
dictionary 
listing 


Write new trim 
solution to F16TI 
Set LWRTR falsi 



Execute code in 
DYNAMIC... END block not 
contained in DERIVATIVE or 
DISCRETE blocks. 


LTRand. ^ 
(ITFLG < 2) 


T>TSTP^ 


Execute 

TERMINAL...END 

block 


Write variables on the 
OUTPUT and PREPAR 
lists to units DIS artd PRN. 


Integrate over the 
communication interval, 
CINT. Call DERIVATIVE 
and DISCRETE code 
blocks as required. 


Write variables on the 
OUTPUT and PREPAR 
lists to units DIS and PRN. 


Return to ACSL 

executive: 

ACSL> 
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Details of the logic used to trim the simulation are shown in figure 2.4. Subroutine TRIM defines 
an array called XTRIM and an array called YTRIM. The first NXTR elements XTRIM are 
varied until the the first NXTR elements of YTRIM are within some tolerance of zero. The 
allowable elements of XTRIM and YTOIM are shown in tables 2.1 and 2.2, respectively. 


T^bIe2T^^clectabl^nemtersoOC^W^^n^ 


1 

MACHTR 

Mach number, n.d. 

2 

BETTR 

sideslip angle, p, radians 

3 

ALFTR 

angle of attack, a, radians 

4 

PIC 

body axis roll rate, rad/sec 

5 

QIC 

body axis pitch rate, rad/sec 

6 

RIC 

body axis yaw rate, rad/sec 

7 

THETR 

Euler pitch angle, radians 

8 

PHITR 

Euler roll angle, radians 

9 

PSITR 

Euler heading angle, radians 

10 

GAMTR 

longitudinal flight path angle, y, radians 

11 

DTVL 

left engine thrust vectoring nozzle deflection, degrees 

12 

DTVR 

right engine thrust vectoring nozzle deflection, degrees 

13 

PCATR 

lateral stick deflection (pilot commanded aileron), inches 

14 

PCSTR 

longitudinal stick deflection (pilot commanded stabilator), 
inrhp^ 

15 

PCRTR 

rudder pedal force (pilot commanded rudder), lbs 

16 

DPSYTR 

throttle position, symmetric, [31,127] 

17 

DPASTR 

throttle position, antisymmetric, [96,-96] 

18 

DSSYTR 

stabilator position, symmetric, degrees 

19 

DSASTR 

stabilator position, antisymmetric, degrees 

20 

DASYTR 

aileron position, symmetric, degrees 

21 

DAASTR 

aileron position, antisymmetric, degrees 

22 

DRSYTR 

rudder position, symmetric (both t.e. in), degrees 

23 

DRASTR 

rudder position, antisymmetric (both t.e. left), degrees 

24 

DNSYTR 

leading-edge flap position, symmetric, degrees 

25 

DNASTR 

leading-edge flap position, antisymmetric, degrees 

26 

DFSYTR 

trailing-edge flap position, symmetric, degrees 

27 

DFASTR 

trailing-edge flap position, antisymmetric, cfegrees 

28 

CSB 

commanded speed brake, degrees 

29 

MUDGTR 

bank angle, p, degrees 

30 


open for future expansion 
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Table 2.2. Selectable members of YTRIM arra 



UD 

d/dt x-t)ody component of inertial velocity, g’s 

2 

VD 

d/dt y-body component of inertial velocity, g’s 

3 

WD 

d/dt z-body component of inenial velocity, g's 

4 

PD 

d/dt body-axis roll rate, rad/sec^ 

5 

QD 

d/dt body-axis pitch rate, rad/sec^ 

6 

RD 

d/dt body-axis yaw rate, rad/sec^ 

7 

GAMZR 

= GAM - GAMTR = y - ypR, radians 

8 

PHIZR 

= PHI - PHITR = y-yj-R, radians 

9 

THE 

Euler pitch angle, radians 

10 

LAMBDA 

lateral flight path angle, X, radians 
aerodynamic ay, g's 

11 

FVRA(2) 

12 

CLERR 

= Q.RFSF - CLTR = actual - target lift coefficient, n.d. 

13 

14,15 

CMRFSF 

steady flow pitch coefficient, n.d. 
open for future expansion 


of XTRIM and YTRIM ^ elected by setting the first NXTR elements of the IXSEL 
and lYSEL arrays. A typical 4x4 longitudinal open-loop trim case could be set up as indicated in 

figure 2.5. In this example, Mach and altitude are selected a priori and the variables {a, 0, 
power setting) are varied until {ii, w, q, y-YrR ) are driven to zero. 


^sy’ 


A(L^iL> NXTR = 4 

ACSL> HIC=10000.,MACHTR=.5, GAMTR=0.0,LFCS=.F 
ACSL> IXSEL=3,7,16,18 
ACSL> IYSEL=1,3,5,7 
ACSL> START 


Figure 2.5. Trim example 1 

Trim ex^ple 2 in figure 2.6 is similar except that a trim with the flight conttx)! system engaged is 
sought. The pilot longitudinal stick command is substituted for stabilator position in the list of 
variables to be varied. 


aO^L> NXTR = 4 

ACSL> MACHTR=.5, GAMTR=().0,LFCS=.T 
ACSL> IXSEL=3,7,16,14 
ACSL> IYSEL=1,3,5,7 
ACSL> START 


Figure 2.6. Trim example 2 

Achieving trim requires that all state derivatives are zero, not just the states associated with aircraft 
motion. The states associated with actuator, sensor, and control law dynamics must be trimmed. 
For the continuous actuator and sensor states, one approach would be to include the.se additional 
s^tes explicitly in the nonlinear trim problem. This would serve to unnecessarily increase the 
dimension of the trim search. In addition, the trim architecture would become dependent on 
actuator and sensor transfer functions. For the transfer-function-based dynamics typically 
asociated with actuator and sensors, the calculation of the required state initial conditions to achieve 
trim for a given input is simple. Therefore, when in the trim loop, the dynamics of the actuator. 
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sensor, and control system are trimmed locally, not as part of the nonlinear search. Integrators are 
turned off ai^ steady state values for all filters are calculated and passed through. Actuator and 
sensor state initial conditions are set so that the stale derivatives are zero for the given inputs. 
However, local trimming of actuator and sensor dynamics creates an implicit loop in the global trim 
problem when the flight control system is engaged (LFCS set true - see section 9.2). For example, 
suppose trim pitch stick position (PCSTR) is one of the admissible XTRIM variables. A change in 
PCSTR produces a change in the control surface commands, which instantaneously produces 
(since Ae actuator dynamics are suppressed) a new contnil deflection. The change in control 
deflection produces a change in the vertical accelerometer output, n^^. Since njj^ is a feedback, a 

change in n^^ produces an additional change in the control surfaces, which produces new, slightly 
different n^^. The trim search assumes the following functional relationship: 

{YTRIM} =/{ XTRIM) 

The effect of the implicit loop is to make { YTRIM ) dependent not only on | XTRIM ) , but on the 
past values of (YTRIM). The degradation the functional relationship will impede or prevent 
convergence. 

The somewhat ad hoc solution adopted by the authors was to implement an inner trim loop to 
converge the implicit loop (see figure 2.4). The inner loop increments a counter, ITR. The 
program loops through the control and flight dynamics code with constant trim driver and initial 
conditions values until ITR = ITRMX. The default value for ITRMX is 6. 

Subroutine TRIM is called after the inner trim loop completes ITRMX passes. The trim routine 
perturbs the control commands and flight states specified for XTRIM and checks for convergence 
of the YTRIM values. When convergence is achieved, ITFLG is set to 3. If convergence fails, 
either because the maximum outer iteration count (MXITR) is exceeded or the sensitivity matrix 
becomes singular, ITFLG is set to 4. ITFLG will be set back to 1 at the first pass of each run (as 
shown in figure 2.4), and trimming will occur each time the program is started at the ACSL 
prompt, unless the variable LTR is set false. If LTR is false, the program will pass once through 
the initial section and go directly into the main program loop, bypassing the trim loop. 

Trim convergence is determined in the ACTRML subroutine, which is called Irom the TRIM 
subroutine. The variable TRTOL (default = 0.00005) sets the tolerance criterion for convergence. 



3. 


MASS AND GEOMETRY DATA 

3.1 General Arrangement 

The general arrangement of a single place F/A-18 is depicted in figure 3.1. This report and the 
MDC documents listed in the References locate points on the airframe using fuselage station (FS), 
buttock line (BL), and water line (WL) coordinates. The (FS,BL,WL) coordinates arc defined in 
table 3.1. At the plane of symmetry, the BL coordinate is zero. 

Table 3.1. (FS,BL,WL) co()rdinates^ffigurejJ___ 

FS = positive aft, inches 

BL = positive out the right wing, inches 

The F/A-18 is a single place, mid-wing, high performance twin engine fighter capable of carrying a 
variety of armament loadings. For the F/A- 1 8, the primary mission is that of "Fighter Escort". 

The Fighter Escort mission loading is two wingtip mounted Si^winder missiles, two fuselage 
mounted Sparrow missiles, and the M61 A1 gun. Each engine is a General Electric, F404-GE-400 
rated at 16,100 lbs of uninstalled, static, sea level thrust. The total internal fuel capacity is 
10,860 lbs while 6,732 lbs of fuel can be carried in external tanks. The design load factors are 
+7.5/-3.0 g’s. The design sinking speed at touchdown is 25.1 ft/sec. 

A moderately swept wing and a highly swept leading-edge extension (LEX) form the prima^ 
lifting surface. Leading-edge and trailing-edge flaps on the wing arc used for carrier operations, 
subsonic maneuvering, and cruise. The flap deflections arc scheduled with M^h number and 
angle of attack to provide improved drag and stability characteristics. All moving horizontal tail 
surfaces are mounted behind and below the wing to provide positive stability at high angles of 
attack. Twin vertical fins arc both canted out and toed out (figure 3.1). 

Pitch control is provided by the horizontal tails. Lateral control is accomplished with a blend of 
ailerons, differential flaps, and differential horizontal tail. Directional control is provided by a 
rudder on each vertical fin. A speedbrake is also available for deceleration control in the 
maneuvering configuration. (MDC A7247, p2-l( 

3.2 Sign Conventions 

While the (FS,BL,WL) coordinate system is u.sed in various MDC reports and the MDC simulation 
to locate positions on the aircraft, the equations of motion make use of a right-hand body-fixed 
frame aligned with the (FS,BL,WL) and typically called the body frame. The body frame origin is 
at the center-of-gravity (c.g.) of the aircraft and is oriented so that Xg is positive forward, Y 3 is 
positive out the right wing, and Z 3 is positive down. For example, a sensor package that is located 
at (FS(j, BL<j, WL 0 ) would have the following coordinates in the body frame. 




- (FSo - FScg) ’ 


> 

BL0 - BLcg *' 



- (WL<y - WLcg) 


where (FScg, BLcg, WLcg) are the coordinates of the c.g. in the (FS, BL, WL) axes system. 
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The body frame, along with the stability frame, control surface sign conventions, and coefficient 
sign conventions, is depicted in figure 3.2. The stability frame is denoted by the subscript "S" in 
3.2. A rotation of a of the stability frame about the Ys axis brings the stability and body frames 
into alignment. The roll, pitch and yaw aerodynamic moment coefficients, Q, and Cp, are 
defined about body frame axes, X 3 , Yg, and Zg, respectively. The lift, drag, and side force 
coefficients, C^, Cj) and Cy. are defined as positive along the -Zj axis, the -Xj axis and the +Ys 
axis, respectively. Note that Yj is aligned with Yg. 



The control surface sign conventions, depicted in figure 3.2, apply to each member of a control 
surface pair and are described in table 3.2. 


^bleT2^^^m^^Su^ceSign^bnve^M« 


5n 

leading-edge flap 

positive leading edge down 

5f 

trailing-edge flap 

positive trailing edge down 

5a 

aileron 

positive trailing edge down 

8 h 

horizontal tail 

positive trailing edge down 

6r 

rudder 

positive trailing edge left 
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FS 687.54 
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WL 100.00 
Thrust Centerline 


This page left intentioiiuily blank. 



3.3 Dimensional Data 

The airplane component dimensional data are given in table 3.2 and pages 2.6-2.S of MDC report 
A7247, vol I. 


Table 3.3 




Dimensional Data 


Total^Splan^^^ 

Net Wetted Area (minus engine nozzles) 
Overall length 
Oerall Height 

Wing 


Area, Sr^f 
A rea (Actual) 

Wetted Area (Including launchers and 
aileron actuator fairings) 

Span, bREF 
Aspect ratio, AR 
Chords 

Root Chord (Theoretical) 

Root Chord (Exposed (g) 54.0) 

Tip Chord 


2028 ft2 
56.0 ft 
15.3 ft 

400 ft2 
401 .7 
562 ft2 

37.42 ft 

3.5 

15.86 ft 
13.19 ft 
5.52 ft 


MAC, Crff 

Leading Edge (L.E.) Sweep 

c/4 Sweep 
Taper Ratio 
Dihedral 
Geometric Twist 

Incidence 
Airfoil Section 
Airfoil Thickness at 
XW 56.876 
XW 145.390 
Tip Chord 
Leading-Edge Flaps 
Type 
Area 
Span 

Deflection (+ leading edge down) 


11.52 ft 
26.7° 

20 ° 

0.35 

-3° 

0.0° at Xw 145.39 
4.0° L.E. down at tip 
0 . 0 ° 

Modified NACA 65 A with sharp L.E. 

5.0 % chord 

3.5 % chord 

3.5 % chord 

Plain 
48.4 ft2 total 
13.8 ft per side 
[0.0°, 34°], maneuvering 
[12°, 34°],Takeoff and Landing 
t 3° differential 


Trailing-Edge Flaps 
Type 


Area 

Span 

Wing Span Range Covered (percent) 


Chord 

Deflections (+ trailing edge down) 


Single Slotted w/Drooped Shroud 
61.9 total 
8.72 ft per side 
19-68.6 % b/2 
30 %c 

[-(-17°, +30°l Takeoff 
[-(-17°, ->-45°], Landing 


PRECEDING PAGE BLANK NOT FILMED 
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Tabic 3.3 Airplane D imensional Data (continued) 

Shroud 


Area 
Chord 
Hinge Line 
Deflection (+ down) 

Aileions 

Type 

Area 

Span 

Ptrcent Wing Span 
Deflections (+ trailing edge down) 

Shroud 

Area 
Chord 
Hinge Line 
Deflation (+ down) 
Leading Edge Extension (LEX) 

Planform Area (Slots Closed) 
Wetted Area 
Leading Edge Sweep 
Incidence (forward of inlet) 

LEX Intersection with Wing 

Horizontal Tails 

Exposed Area 
Wetted Area 
Aspect Ratio 
Taper Ratio 
L.E. Sweep 
c/4 Sweep 
Dihedral 
Span, b^Ep 
Total Span 
Chords 

Root Chord (Exposed) 

Tip Chord 

MAC, Ci^-p 

Tail length (c/4 to Cht/4) 

Airfoil section 
Airfoil Thickness at 

R(X)t 

Tip 

Hinge Line 


13.9 total 
7%c 

68.5 % c 
[10*,+30’] 

Single Slotted w/ Drooped Shroud 

24.4 total 
5.68 ft per side 
68.8 -100% b/2 
I -25’, +45’), Takeoff and Landing 
f-25’, +25*], Maneuvering 

5.6 total 
7%c 
68.5 % c 
MO’,+30’] 

56.0 

210.0 ft2 
43° 

6 ’ 

F.S. 403.9 
B.L. 54.0 
W.L. 108.2 

88.1 ft2 

176.0 ft2 
2.4 
.46 

47.2’ 
42.8’ 
- 2 ’ 
14.67 ft 

21.6 ft 

8.23 ft 
3.79 ft 

6.28 ft 

16.81 ft 

Mtxlified NACA 65A with Sharp L.E. 
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6 % Cj.pi’ 
2 % C|.|Y 

13 % Cfpp 



Tlwizontall Tails (continued) 

Deflections (+ trailing edge down) 

f-24", +8“], Symmetric 
1-24', +10.5*1, Max 


Vertical Tails 

Area 52.0 ft^ each 

Wetted Area 1 04.0 ft2 each 

Aspect Ratio 1-2 

Taper Ratio 0 -'^ 

L.E. Sweep 41.3° 

c/4 Sweep 35.0* 

Cant (Tip out) 20* out 

Toe Out 1 * 

Panel Span 7-92 ft 

Chords 

Root Chord (Expo.sed) 9.42 ft 

Tip Chord 3.75 ft 

MAC, CvT 6.99 ft 


Tail length (c/4 to cvx/4) 
Airfoil section 
Airfoil Thickness at 


10.18 ft 

Mtxlified NACA 65A with Sharp L.E. 


Root 


5 % Cyx 


Tip 3 % cvT 

Rudders 

Area 111 . ft^ each 

Span 5.21 ft 

Percent Chord 20 % c 

Deflection - 30° 

Fuselage 

Length 53 ft 

Maximum Width 7.6 ft 

Wetted Area (includes nacelles, canopy, and nozzles) 890 ft2 

Sf>eedbrake 

Planform Area 1 3.9 ft^ 

Span 2.5 ft 

Chord 5.57 ft 

Maximum Deflection (relative to ML) ^ 



3.4 Aerodynamic Reference Dimensions 

The aerodynamic reference dimensions are given in table 3.4 below and in MDC report A7247, 
page 3-2. 


Table 3.4. Aerodynamic Reference Dimensions 

Wing Area, Sref 

Wing MAC, CREF (L.E.MAC @ FS 423.99) 

Wing Span, 1 >ref 

HcMizontal Tail Moment Arm, l^j (c/4 to Chj/4) 

Horizontal Tail Volume Coefficient, 

Vertical Tail Moment Arm, lyr (c/4 to cvt(4) 


138.275 in = 11.523 ft 
37.42 ft 

201.76 in 

.321 
122.1 in 


Vertical Tail Volume C!oefricient, Vyj 


.0707 


3.5 Weight, Center of Gravity, and Inertia Data 

The weight, center-of-gravity location, and mass distribution data for the default simulation 
configuration and for a fighter escort configuration are shown in Table 3.5. The fighter escort 
description, along with several others, is found in figure 2.3-1 of MDCT report A7247. The default 
configuration represents an early estimate of an FI 8 HARV configuration after the addition of 
thrust vectoring hardware. Of all the configurations listed in figure 2.3-1 of MDC report A7247, 
the data fw the fighter escort configuration with 60 percent internal fuel most closely matches the 
default values. 

Table 3.5. Summary of Simulation Weight, Center-of -Gravity, and Inertia 

""^"'''""'''’^^™"^^^^*c!g^ocanons***^'^*''Moment^n3ftoducts^nnerta^^^ 
Weight FS WL Ixx lyy ^zz ^xz 

(lbs) (in) (slug - ft^) 


Default 

~imr 

455.0 

102.8 

23OO0 

151293 


-2971 

Fighter Escort 
(60% internal 
fuel) 

31665 

457.3 

101.6 

22337 

120293 

138945 

-2430 


The Ixz product of inertia in table 3.5 is defined as the integral over the volume of the aircraft of the 

product (xz) (body frame ctxirdi nates) weighted with the mass density, a, of the aircraft. The 
negative of this quantity is what appears as the (1,3) and (3,1) elements of inertia matrix. Thus, 

Ixz =nJoxzdV (3.2) 

The product Ixz typically positive for aircraft with vertical fins but is negative in this case. 
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4. EQUATIONS OF MOTION 

The derivation of the equations governing the motion of a rigid airplane over a flat, nonrotating 
Earth has been performed many times in the literature and is not likely to be improved upon here. 
Readers interested in the derivation of said equations are referred to Etkin (1972) and Roskam 
(1979). In the flSbas simulation, mass and inertia properties are constants and do not vary with 
time If desired, the mass and inertias could be made functions of fuel load according to data 
presented in MDC-A7247 and MDC-A8575. Gravity is assumed to be aligned with the +Z axis of 
a local (Earth) right-handed reference frame fixed at sea level. Gravity is assumed constant over 
the airplane volume, which results in gravity acting as a force at the center of mass and producing 
no moments. In fact, gravity is constant with altitude with the default value being the sea level 
value. A reference frame fixed to the body of the airplane (b^y frame) is the most convenient for 
expressing the applied aerodynamic and thrust loads. As indicated in figure 3.2, the body frame 
origin is at the center of gravity and the frame is oriented .so that +Xb is out the nose, +Yb is out 
the right wing, and +Zb is out the belly. The airplane is assumed symmetric about the plane 
formed by the Xb and Zb axes, so that the I^y and the ly^ products of inertia are zero and are not 

included in the simulation equations. 

4.1 Translational Equations 

Let the translational inertial velocity of the airplane expressed in body frame components be given 
by [u,v,w]’ and the rotational inertial velocity of the airplane in body frame components be given 
by [p,q,r]’. The equations governing translational motion can then be written as. 


r • 'i 


m 

0 

0 " 


• 

u 


" 0 

-r 

q " 


m 

0 

0 ■ 

[“1 

0 

m 

0 

< 

• 

V 

^ + 

r 

0 

-p 


0 

m 

0 

V- 

_ 0 

0 

m _ 




- -q 

P 

0 _ 


_ 0 

0 

m _ 

lw> 



fFxl 


rF 1 


rF 

^x 

= 

Fy 

► + < 

Fy 

► + < 

Fy. 



G 

[FzJ 

A 



where the subscripts "G","A", and "E" denote gravitational, aerodynamic and thrust-induced 
forces. The components of the gravitational force are given by. 


fFxl 


*xz 

r - sin(O) 

fy 

^ = -mg^ 

^yz 

> = m g cos(0)sin(<|)) 

i^zj 

G 

.^zz. 

Lcos(0)cos(<I))j 


where II 1„, UJ' is the third column of the direction co.sine matrix describing a transformation 

^ XZ XtZi 

from local frame to btxly frame components and 0 and <I» are the Euler pitch and roll angles, 
respectively. The aerodynamic force components are given by. 
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(4.3) 


fPxl 


fCx] 

fy 

► = qS ^ 

Cy. 

[FzJ 

A 

Cz. 


The thrust-induced forces are defined in equations (6.8) of section 6. 


4.2 Rotational Equations 

The body frame rotational equations can be written in a form identical to equation (4.1), with the 
understanding that moments due to gravity are zero. 


d 

dt 


Iv 0 -1 


*x 
0 I 


-I 


xz 


y 

0 1 


xz 

0 


Z “■ 



0 -r q ' 

r 

[P 1 



br 

r 0 -p 



. -q p 0 . 



r Iv 0 -I 


0 I 


-I 


xz 


y 

0 


xz 

0 


L 



fMxl 


fMxl 


r^xi 

My 

► + < 

M 

► + ^ 

My 

My 

G 

[MzJ 

A 

M 


and where the aerodynamic moments are given by. 


fMx' 

pc, ' 

.M 

^ = «ls]cC^ ^ 


m 

A U C J 


(4.5) 


The thrust-induced moments are defined in equations (6. 10) of section 6. 

4.3 Combined Accelerations 


If the body frame components of translational and rotational velocity are combined into a 6x1 
generalized velocity vector. 





q 


^ T j 

then equations (4.1) and (4.2) can be combined and written as. 


(4.6) 


|M] fV) + |f2](M| {V| =|F) +{F) +{F) 

* 1 O A E 

where the inertial mass is given by. 


(4.7) 
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[Mlj = 


m 0 0 0 0 0 

0 m 0 0 0 0 

0 0 m 0 0 0 

0 0 0 L 0 

X Xi 

0 0 0 0 L, 0 


0 0 0 0 U 


and the matrix representing body frame rotation is. 


" 0 -r q 0 0 0 “ 

r 0 -p 0 0 0 

-q p 0 0 0 0 

|Q1 = 

0 0 0 0 -r q 

0 0 0 r 0 -p 

- 0 0 0 -q p 0 - 

In order to achieve an explicit definition of the airplane accelerations, equation (4.7) must be 
modified. Observe that the time rate-of-change of angle of attack in still air, given by. 


- 2 . 2 

U'^ + 

is clearly a function of body frame accelerations. In the F/A-18 aerodynamic build-up, as is the 
case for almost all such aei^ynamic models, the coefficients of lift and pitch contain contributions 
by a nondimensionalized form of the angle-ofiattack rate-of-change 


Cx = ■■ - C,. 


where the symbol ’X’ in equation (4.1 1) can be replaced by 'L' (lift) or 'm' (pitch). In the F/A-18 

simulation, Cv. is a function of Mach, altitude, and angle-of-attack. The implication of equations 
a 

(4.10) and (4. 1 1) is that the generalized aerodynamic force column vector appearing in equation 
(4.7) should be written as. 


(F)^ = {F} I . -|MI (V) 
{V)=0 


where. 
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(4.13) 


[MIa 


a{v) 


(Ml represents an increment to the inertial mass matrix of equation (4.7). If we further define a 
A 

generalized inertial force as. 


{F}j = -|Q)|M|j(V) (4.14) 

and combine equations (4.7), (4.12) and (4.14), we get the equations of motion as implemented in 
the simulation. 


where the total mass matrix is: 


|M1(V) 


IF) 

* 'l+G+A+E 


|M1 = |M|| + (M|a 


(4.15) 


(4.16) 


Equation (4.15) must be solved each simulation time step for the unknown accelerations in ( V) . 
Recall { V) = [u,v,w,p,q,rl' (equation 4.6). If the mass matrix, [M], were constant, then the 

inverse, [M]’*, could be formed at time zero and a simple left multiplication of the right-hand-side 

forces by [M]*^ would produce the accelerations. The aerodynamic part of the mass matrix, [M]y^, 
is a function of angle of attack, Mach number, altitude, and dynamic pressure and is therefore time- 
varying. Forming the inverse of [M] each time step would be prohibitive. T^e approach used in 
the simulation documented in this report is to solve (4.15) using Gaussian elimination. In a typical 
application of Gaussian elimination, row and column switching (pivoting) is performed to insure 
that the largest possible elements arc used as divisors. Both partial and full pivoting s^tegies 
exist. In a general application, if the diagonal elements arc used as divisors without pivoting, the 
divisors may become zero. For most airplanes, and for any configuration of the F/A-18, the 
inertial mass matrix, [M)j, is sufficiently diagonally dominant and [Ml^^ is sufficiently small 
throughout the flight envelope, that the diagonal elements of [M] can be used as divisors with 
impunity. A custom 6x6 equation solver that employs no pivoting strategy is used in the 
simulation to solve (4.15). The custom solver is implemented in subroutine XRBVDC. In order to 
improve efficiency, the use of DO-loops in XRBVDC was minimized. The simulation was tested 
using both the custom solver and a general purpose subroutine that employed a full pivoting 
strategy (the most robust form) with double precision arithmetic. The custom solver gave identical 
answers to 5-6 decimals and achieved a substantial reduction in total simulation run time. 

An alternate approach to the solution of (4.15) is to use a series approximation to the inverse of 
[M]. This approach is valid if the time-varying part of the mass matrix is small. Thus if, 

[Ml = (MjllI-t-AM) (4.17) 

then, 

[Ml * = (I + AMf’lMj!*' » lI-AMl|Mjl * (4.18) 
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if [AM] is small. The inversion of [M^] could be done once in the initialization portion of the 
simulation. 

4.4 Body Orientation Parameters 

Three equivalent sets of parameters to describe the orientation of the body reference frame with 
respect to the the Earth frame are maintained in the simulation. The three sets are the Euler angles 

(\|f, 0 , 4 )» the quaternions (eQ, e^ e 2 , 63 ), and the elements of the direction cosine matrix l^y, 
..., 1^). These parameters are cdculat^ by two separate methods according to the logical switch, 
LEULER. If LEULER is selected, then the three sets of orientation parameters are determined by 

integrating the time derivatives of the Euler angles, \|/, 0, 0. If LEULER is not selected, then all the 
orientation parameters are determined from integrating the time derivatives of the quaternions 
(Euler parameters). The quaternion method was implemented first and avoids the problem of a 

singularity that occurs in the Euler angle derivatives when the pitch angle, 0, is 90°. The 
singularity is associated with 0 because of the order of application of the Euler angle rotations in 
transforming from the Earth frame to the body frame, that order being heading (y), pitch ( 6 ), and 

roll (4). If the 0 = 90° region is avoided, the two methtxls should produce the same answer to 3-4 
decimals. The Euler angle integrations were added later so the Euler angles would occur as states 
in the simulation and therefore linear analysis using the Euler angles as state variables could be 
easily perframed. Time derivatives are calculated and integrated for both the Euler angles and the 
quaternions at all times in the simulation. The effect of LEULER is to select one set of integrated 
outputs and ignore the other. Equations 4. 19 through 4.3 1 below can be found among or derived 
from the equations in Appendix E of Gainer and Hoffman (1972). 

4.4.1 Euler Angles 

The time derivatives of the Euler angles are given by; 


<J) = p + (cos(0)} ' (sin((t>)sin(0)q + cos(<|))sin(0)r) 
0 = cos(4>)q - sin(<|))r 
\|/=|cos(0)} (sin((]))q + cos(<l))r) 




(4.19) 


In the simulation implementation of (4.19), the following expression is used for cos(0) where 
division by cos(0) is called for. 




101 


if lcos(0)l < e = 1 0 ' ® 


(4.20) 


cos(0) otherwise 


so that division by zero when the pitch angle is at 90° is avoided. States created by integrating the 
derivatives arc subscripted with "x" to denote the fact that they may be ignored later depending on 
the choice of LEULER. 
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(4.21) 


^ 

y X = J Vdt + V 
0 




1C 


0x = J Odt + 0j, 


0 = J ^dt + <|> 


1C 


4.4.2 Quaternions 

The derivatives of the quaternions are given by: 


r . 


®0 “Cl “C2 -63 

Co 



• 


®i ®o -63 C2 

J Cl 



• 

C2 


C2 C3 Cq -e| 

• 

V e3 J 


_®3 -C2 Cj eg ^ 


0 

P 

q 

r J 


Define 


e»eo + lci+je2 + ke3 = 



^ ^ and w = ip + jq+ kr = 


0 1 
P 

q 

r J 


where die following rules of multiplication govern the i, j, k basis vectors. 

Ij = -jl = k, jk = -kj = i, ki = -ik = j, and ii = jj = kk = -1. 

Note that the amiplex conjugate of e, the quaternion vector, is given by e = Cq - iej 
Observing that the product ew is given by, 

-pei-qe2-rc3 
+ i[peo - qc3 + rej] 

+ j[pe3 + qeo-reil 
+ k[-pc2 + qej + rcol 

Equation (4.22) can be rewritten as, 

e = ew or e = we where e = 


fcol 


fcol 

Jc, 

\ e = ^ 

■®1 !► 

®2 


-C2 

IcsJ 




(4.22) 


-je2-kc3. 


(4.23) 


(4.24) 
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where left multiplication by e is equivalent to left multiplication by the 4x4 matrix of equation 
4.22. The components of the quaternion rate, e, are then integrated to produce quaternion states. 


JOv 


= J Cj dt + ej.^ (j = 0,3) 


(4.25) 


The quaternions are used as the generating set from which the other two equivalent set of 
(orientation parameters are deriv^. The quaternions in their unnormalized form are calculated 
according to the value of LEULER. The additional zero subscript on the quaternions serves to 
indicate their unnoimalized status. 


Cn_ = 




% 


eo„ = 


' 2 o 


— 


'^0 


\g„ 0„ d)„ 0„ <))„ 

) cos(-^)cos(-y)cos(-^) + sin( -^)sin(-^)sin(-^) 

LEULER=.T. 

(4.26.a) 


LEULER=.F. 

r ¥« o„ Vx Ox 

1 cos(-^)cos(-^)sin(^) + sin(^)sin(-^)cos(^) 

LEULER=.T. 

(4.26.b) 

riO 

LEULER=.F. 

^ cos(^)sin(-^)cos(^) + sin('^)cos(-^)sin(^) 

LEULER=.T. 

(4.26.C) 


LEULER=.F. 


Vx Ox <1>X Vx Ox . <|)x 

1 sin(-^)cos(-^)cos(^) + cos(-^)sin(-^)sin(^) 

LEULER=.T. 

(4.26.d) 


LEULER=.F. 



The quaternions arc then normalized. Note that if \|/, 0, <[) are small, then to first order, \j/ = 63 , 0 = 
Cj. 0 = Cj, and Cq = 1 . The norm of the quaternion vector is formed. 


and then used to scale the quaternion components: 

ej = ejy/lel 0=0.3) 


1/2 


(4.27) 


(4.28) 
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While not required for the following development, it is interesting to observe that: 


riel^l 
0 
0 
0 


> = 


= ee 


4.4.3 Direction Cosine Matrix 

The direction cosine matrix elements and the final Euler angle can then be derived from the 
normali29Bd quaternions: 



where in expanded form, the first equation of (4.30) is given by. 


Cq Cj 62 63 

-Cj Cq 63 -C2 

"^3 

— ^2 1 ^0 


eo -63 -e2 <3 
®3 

63 -62 Cj Cq 



0^ 

1 

< > 

0 

()> 


Finally, the Euler angles can be calculated directly from the elements of the direction cosine 


V = 



0 - sin { *xz } 

f = lan'M 

I zzj 


> 


(4.29) 


(4.30) 


(4.31) 


matrix. 


(4.32) 
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5. AERODYNAMIC MODEL 
5.1 Aerodynamic Data Sources 

The aeiodynamic characteristics of the F/A-18 arc discussed in detail in MDC A7247, vol. 1.2. 
The low a report, (MDC A7247, vol. 1), describes the force and moment build-up for thrce basic 

configurations (maneuvering, carapult a^^^pptoach) in the a^ ■ 

The maneuvering configuration is described in section 3 of MD ^ t _ a7?47 

sneed/hieh lift configurations (approach and catapult) arc desenb^ m secttwi 4 of MDC A7247, 
5 SlMDC A7247fwl. 2 describes the force and moment build-up for the maneuvering 

configuration in the angle-of-attack range of +40° to +90°. In the MDC bateh 
fdisc^^ in section l72 of this report) four separate force and moment build-ups, 
aDDioach” "uike-off ', and "high" are implemented. The first three configurauons of the MDC 



7, vol. z. me mancuYcniig wiuiguiauv/-. 

* confiEuration cciuations of the MDC batch i 

Sightly SfcrcnL The MDC A7247, vol. 1 "maneuvering" build-up h^ 
while the MDC batch simulation code "clean" build-up has incremente due to 8®^ P 5^ . 

The simulation aerodynamic model described in this report is a ^bination of the maneuvenng 
S^fiSSrSE MDC A7247. vol. 1.2 and in the MDC 

few small extensions extracted from the high-lift configurations. In addition, MDC report A85 /b 
COTtLns background material, including a history of configi^tion (tevelopment, a history of the 
substantiating wind-tunnel testing, and propulsion system characteristics. 

As reouired by the simulation equations of motion discussed in section 4 of this rcpoit, 
sien^namic force and moment coefficients arc calculated in two parts. The portion of the total 
coeffiaent that is an explicit function of airplane velocities or control surfece positions, the steady 
fl^wSs cicuSSi^^^^^^^^ SFAEiJrf. The part of the total coefficient th^ is propomonal 

to aimlatie acceleration, the unsteady part, is calculated in subroutine USAERRF,_pie steady flow 
resuk is the sum of the static components plus the effect of b^y angul^ rates. The unsteady part 
includes the effect of the time rate of change of angle of attack or sideslip. 

The moment coefficients and input velocities arc defined about the aerodynamic reference center, 
(a.r.c.), which is defined in Table 5.1. 


Fuselage Station (FS), in. 

458.56 

(+aft) 

Buttock Line (BL), in. 

0.00 

(+ right) 

Water Line (WL), in. 

100.00 

“P> 


The (FS/BUWL) coordinates define a body-fixed axes systein. TJe 

eouations of motion is centered at the center of gravity, aligned with the 

Krsign convention of +x forward, + y out the right wing, and +z down (figure 3-1). 

There are two logical variables that affect the coefficient calculations in the simulation ^fcn^d i^^^ 

If the quasi-static-elastic logical. (LQSE). is true, then all the flex/ngid 
SSIS Sd that were in the original model remain in effect. If LQSE is 
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false, then all the flex/rigid coefficients are set to the values they approach as air density, and 
therefore deformation, goes to zero. Setting LQSE false causes the aerodynamic forces and 
nioments appropriate to the wind-tunnel-shape, which is not necessarily the jig-shape, of the 
airplane to be cdculated. The LQSE option allows the use of the same code in a simulation where 
the flexible modes are actively modeled. The second logical is the real-time-equivalent flag (lrtE). 
If LRIE is true, then various differences that exist between the aerodynamic model used in the 
LaRC real-time simulation (dmsflS) and the ACSL simulation described in this report are 
removed. These differences are discussed in section 5.5.2 "(Leading-Edge Flap Effects (Y,l,n)"), 
section 5.5.5 ("Rudder Effects"), and section 5.5.9 ("Effect of Airplane Asymmetries (Y,l,n)") of 
this report. 

'rae difference between the two aerodynamic models arises from changes made to iheflSbas 
simulation math model that were not duplicated in the LaRC real-time simulation. The first of these 

changes was the inclusion of a rolling moment increment, that is a function only of angle 

of attack, ot, and nonzero for ot > 40°. The effect is to intrcxluce a rolling disturbance as the 

simulated airplane moves through the 40° to 60° angle-of-attack region. The term arises from nose 
asymmetries and is discussed on pages 4-20,21 of MDC A7247, vol. 2 and section 5.5.9 ("Effect 

of Airplane Asymmetries (Y,l,n)") of this report. The dmsflH simulation has the term ACi 

‘asm 

zeroed out. The second change is the inclusion of a nonzero estimate for (lift coefficient per 

"n 

unit leading-edge flap) for the case of 40° < a < 45° and leading-edge flap position different from 
the scheduled value. This change is discussed in section 5.3.4 ("Leading-Edge Flap Effects") of 
this report. The third difference is the inclusion of rudder toe-in effects on drag, pitch, and lift in 
the ACSL simulation. Rudder toe-in can induce significant pitch moment and is an appropriate 
effect to include in a simulation suppcMiing agility research. Symmetric rudder effects on pitch are 
discussed in section 5.3.9 ("Rudder Effects") of this report. 

In the reconstruction equations, the total aerodynamic force and moment coefficients are built up 
fewn a l»ge collection of intermediate coefficients. An attempt has been made to adopt symbolic 
conventions that indicate the source and purpose of the intermediate coefficients. Typically, the 
total force ot moment coefficient consists of a "basic quantity" plus actional coefficients. The 
"basic contribution" is denoted by a zero subscript and is a function of Mach and either angle of 
attack or sideslip. For the "basic contribution", the controls are at a nominal position (typically 
zero) and the airplane is undeformed. Additional coefficients that represent increments to another 

coefficient (basic or otherwise) will begin with 'A'. Increments due to 'x' are the result of 
evaluating a tabular function where 'x' was an input. Coefficients that resemble classic stability 
derivatives and are to be multiplied by the appropriate control or nondimensional quantity will have 

no 'A' symbol. 

5.2 Control Surface Sign Conventions 

The control surface sign convention used in the aerodynamic model follows the convention 
established in MDC A7247, vol. 1 . The terms "stabilator" and "horizontal tail" are used 
interchangeably in various MEXT documents. MDC stability and control reports that describe the 
aerodynamic model (MDC A7247, vol. 1,2), tend to use the term "horizontal tail". MDC reports 
describing the flight control laws (MDC A4107, vol. 1 and MDC A7813, vol. 1,2,3), use 

"stabilator". "Stabilator" will be used in this document. In addition, the symbols 5j^ and 
used interchangeably in MDC A7247, vol. 1 to refer to "horizontal tail" deflection. 
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®Hr.l 

^R.L 

^R,L 

Sf 

Sn 

A5h 

A5f 

A5„ 

5rt 

«R 

^f/A 


^^^^_Jabl^5^2^ConOTolJDeflection^Si^ 

= right or left stabilator (horizontal tail), + trailing edge down 
= right or left aileron, + trailing edge down 
= right OT left rudder, + trailing edge left 
= right or left leading-edge flaps, + leading edge down 
= right or left trailing-edge flaps, + trailing edge down 

= K*Hr*8hl) , symmetric (average) stabilator 
= ^ symmetric trailing-edge flap 

= ^ ^^or symmetric leading-edge flap 

= 5hl‘V’ differential stabilators (horizontal tails) 

= 8f - 8f , differential trailing-edge flaps 

JL^ lx 

= differential leading-edge flaps 

~ ^ { ^Rr 1 ’ ™dder (both t.e.in) 

^’^^symmetric rudder (both te.left) 

= \ {^fR-"^AR + \+SAL} , combined t.e. control 


5.3 Lift and Pitch (L,m) 

The reconstruction equations that are used to compute the total lift and pitch coefficients are almost 
identical and are combined in a common form in equation (S.l). To the basic lift and pitch 
reconstruction equations (MDC A7247, vol. 1, section 3.1.1) the following have been added; (1) 
landing ^ar effects from the MDC batch simulation code (lines 805-906 of subroutine 
AR018CN), and (2), rudder toe-in effects extracted from section 4.1.10 of MDC A7247, vol. 1. 
Replacing "X" in the equation below with "L" or "m" as appropriate will give the lift and pitch 
equations, respectively. The total lift/pitch coefficient is given by. 



where the steady flow lift/pitch coefficient is given by. 


(5.1) 
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^Xsf " 2 


+ R. 


FLT 


'L5, 


ACxj^^ + ACxg^^ 


+ Cvc 5rn + ACv^ 8n +i ^^XSjrj^f 


'X6n"ni 


X5n-n 


+ ACvo + ACyc S| rj 

^5 s 3 ^Olg ^Or ^Orj 


ACy +(Cy +ACyl^ 

^stores I ■^q Xqj 2V 


(5.2) 


while the unsteady component of lift/pitch, i.e., proportional to nondimensional rate of change of 
angle of attack is given by, 



+ ACy. 


(5.3) 


The individual components of the lift and pitch reconstruction equations, (5.2,5.3), are defined in 
sections (5.3.1) through (5.3.1 1). These components are for the most part calculated by 
performing linear interpolations according to tabulated functions. The variables used as input to 
die interpolation functions are normally limited to prevent extrapolation. Thus once the actual state 
or control variable reaches the limit of the data, no further change occurs in the coefficient. The 
input variables and the limits (if applied) are shown for each function look-up. 

5.3.1 Basic Airframe Coefficient (L,m) 

The basic airframe lift/pitch coefficient for the rigid airplane with control surfaces undeflected is 
given by. 


C 


[ Fcxo<«f^?2 ’ ^^^2^ a < 40°; 5n = 0° 
|fhcxo(« 140-M|'2') a>40°;5„ = 34° 


(5.4) 


where the notation is defined as follows, 

( the function F^xo defined by data tables 
is evaluated at an angle of attack and Mach number 
that are limited according to the brackets. 

At low angle of attack (a < 40'), Cx^ refers to a zero stabilator, retracted flaps configuration. At 
high angle of attack (a > 40*), basic lift/pitch, Cx^, refers to a configuration with zero stabilator 
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and a leading-edge flap deflection of 34°. The elastic increment to the basic airframe lift/pitch 
coeffrcient is: 




^DCXOF^^l^’ ^1.6’ ^ 


60000 


0.0 


LQSE=.T. 


LQSE=.F. 


(5.5) 


In figures 3.1. l-I through 3.1.1-22 of MDC A7247, vol. 1, the coefficients Cx^ and are 
defined and combined for presentation into one coefficient. 

S.3.2 Stabilator Effects (L,m) 

Increments due to left and right stabilator (horizontal tail) deflections are. 


AC 



Jf^DCXDH -24 *^1.2^ a <.40 

I^Phdcxdh (“I 40 ’ rI .2 ^ a > 40 


(5.6) 


The elastic/rigid stabilator effectiveness factor is. 



f^FRCXDH ( *^i^2 ’ 0 ^ LQSE-.T. 

1.0 LQSE=.F. 


(5.7) 


5.3.3 Aileron Effects (L,m) 


The lift/pitch increments due to left and right aileron deflections are. 


AC 



pDCXDA(“f^ ’ ^Al/r^ -25 ’ .2 ^ 

PhE)CXDa(®[40 ’ ^Al,r^-25 ’ .2 ^ 


a <40° 
a > 40° 


(5.8) 


The elasticMgid aileron effectiveness factor has been replaced by a flight-test derived quantity, a 
"flight/rigid" factor. This factor is larger than unity in the "rigid" limit because it reflects the 
configuration change that occurred between 1979 and 1984 in which the aileron span was 
increased outboard to the wingtip, (MDC A8575, page 1-4). The data in Fjx^xda f^HDCXDA 


are 


FLT 

for the small aileron (MDC A7247, vol. 1 , page 3-357). Note that R is used in both the 


lift (Cj^j) and pitch equations. 


FLT 



.>FFCLDA(Ml'a.Hp|“““®) LQSE=.T. 

.1.25 LQSE=.F. 


(5.9) 
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5.3.4 Leading-Edge Flap Effects (L,m) 

The assumption in MDC A7247, vol. 1 ,2 and the MDC batch simulation code is that the leading- 
edge flaps are commanded according to the maneuvering flap schedule and will remain fixed at 34° 
if angle of attack is above 40°. If angle of attack is above 40°, Cx^, defined in equation (5.4), 
includes the effect of the nonzero leading-edge flaps. Accordingly, , the lift/pitch per degree 

of leading-edge flap, is not defined in MDC A7247, vol. 2 and is set to zero in the MDC batch 
simulation code if a > 40°. As a result, when the simulated airplane transitions from the low 
angle-of-attack database (a < 40°) to the high angle-of-attack database (a > 40°), a potential 
discontinuity can arise. If the leading-edge flaps are at 34° when the simulation goes through 
tt=40°, no discontinuity in Cj^ occurs. However, if the leading-edge flaps are not on schedule, 
then a discontinuity in C^^ exists. To remove this potential discontinuity in the flSbas simulation, 

is not set to zero abruptly at (X=40°, but rather is ramped to zero as angle of attack goes from 

40° to 45°. Furthermore, when angle-of-attack is above 40°, C^j. is multiplied by the difference 

between leading-edge flaps and 34°, the position at which the high angle-of-attack Cv da»a was 
taken, instead of the total 6^. Thus, 



j f^CXDN ^ 

I (^)li Pcxdn<“=40”, Ml ) 

lo.o 


a <_ 40 

a > 40; LRTE=.F. 
a > 40; LRTE=.T. 


(5.10) 


where 


8 „, = 


Sp a <_ 4 0 
5„-34 a > 40 


The elastic increments to the leading-edge flap sensitivity for lift and pitch are. 


and 




M - .6\i 


.2 

0.0 




60000 


) LQSE=.T. 
LQSE = .F. 


AC^j. = ^ 

mSn 


^(^);,FDCMDNF(«Pi.M(y,Hp|« 
0.0 


60000 


) LQSE=.T. 
LQSE=.F. 


(5.11) 


(5.12) 


(5.13) 
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The factor ^ serves to ramp in the full value of as Mach goes from 0.6 to 0.8.. 

Similarly. ramps in the full value of AC^g by going from zero to one as Mach goes 

from 0.2 to 0.6. 

5.3.5 Trailing-Edge Flap Effects (L,m) 

The lifr^tch per degree deflection of symmetric trailing-edge flaps is. 


'X8f 


r40 

CXDF ^1-4 

0.0 


(al7. Ml * 2 ’ ) a <40° 
a > 40° 


(5.14) 


The elastic increment to the trailing-edge flap sensitivity is, 

^DCXDFf(°^I 4^ ’^f .6 a<40; LQSE=.T. 


^^X5f “ ' 


(5.15) 


0.0 


otherwise 


5.3.6 Speed Brake Effects (Lytn) 


The increments due to speed brake deflection are slighdy different for lift and pitch. The lift 

g 

increment is to be multiplied by normalized speed brake deflection, For the low angle-of- 
attack case, the pitch increment is itself a function of speed brake and no multiplier is needed. 


AC 


lsb 


(lf> 


DCLSB ( ® ^ 


(5.16) 


HDCLSB ^ ^ ^ 


AC 


"™SB 



DCMSB ( ’ ^SbIq" ) a ^ 40* 


1.6 


r60 


(a|2® ,MfV ) a > 40* 


1.2 


HDCMSB ^ ^*140 


(5.17) 


.2 
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5.3.7 Stores (L,m) 

The increment due to stores is discussed in MDC A7247, vol. 1, section 3.1.8, for the cases of 
Escort, Escort Overload, and Interdiction mission loadings. The quantity is currently set to zero in 
the fISbas simulation described in this repon. 


ACy = 0.0 
■^stores 


(5.18) 


5.3.8 iMnding Gear Effects (L,m) 


Increments to lift, pitch, and drag due to landing gear deployment are discussed on pages 4-75 and 
4-76 of reference MDC A7247, vol. 1 . In MDC A7247, vol. 1, landing gear increments are 
included only for the high-lift (take-off and approach) configurations. In the MDC batch 
simulation, landing gear increments are includi^ in the lift, pitch, and drag build-up for the 
maneuvering ("clean") configuration. The increments used in the MDC batch simulation for the 
maneuvering configuration are different than those shown on 4-76 of MDC A7247, vol. 1 . The 
increments used in ihefISbas simulation described in this report are those of MDC batch 
simulation. These increments are to be multiplied by normalized landing gear deployment. While 
studies to date have had no requirement for the aerodynamic effects of the landing gear being 
included, there is virtually no extra computational cost as the code calculating the coefficients is 
bypassed if the gear is retracted. 



FcxlG<“I-2) “<40° 

0.0 a > 40° 


(5.19) 


5.3.9 Rudder Effects (L,m) 

The increment to lift due to antisymmetric rudder is included in the high lift configuradons 
described in section 4. 1 . 10 of MDC A7247, vol. 1 , but not in the maneuvering configuration. The 
lift increment is set to zero in the fISbas simulation described in this report. A nonzero pitch 
increment is included in the flHhas simulation, 


AC, - = 0.0 

l8r 


(5.20) 


AC 


mSR 




f^DCMDR ( • I^rII 0^’ ^ 

f^HDCMDR ( ® I 40 ’ I^rI ^ 0^’ ^^^.2 ^ 


a < 40° 

a > 40° 


(5.21) 


Symmetric rudder is used in approach and landing, but is not part of the nominal maneuvering 
configuration model. While the pitch channel of the flight control system makes no use of 
symmetric rudder in the auto-flap-up (AFU), maneuvering mode, future agility-oriented research 
control laws may employ such a command. Therefore, increments in lift/pitch due to symmetric 
rudder (toe-in) have been added to the maneuvering configuration by using the increment 
developed for die high-lift configurations and discu.ssed in section 4. 1.1 Oof MDC A7247, vol. 1. 
The increments due to symmetric rudder for the flHhas simulation are. 
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ACxjj — ^ 

Srt 


*^DCXDRt( ’ ^RtI -15’ ^ ^ 40 

0.0 a > 40° 


r45 


(5.22) 


Since the rudder toe-in coefficient was developed for the low speed/high lift configurations, there 
is no flex^gid increment or effectiveness factor (page 4-5 of MDC A7247, vol. 1). 

5.3.10 Pitch Rate (L,m) 

The pitch rate derivative for lift/pitch is, 

fFcXQ<“l-4"- ) a £.40° 


Cx — 1 

Aq 


Phcxq ( ®I40 ’ ^ 


1.1 


40‘ 


ACx = 

Xq 


The elastic increment to the pitch derivative for lift/pitch is, 

[pDCXQR(Ml^^HplT^) L0SE = .T. 

lo.o LQSE-.F. 

5.3. II Angle-of’Attack Rate (L,m) 

The angle-of-attack rate derivative for lift/pitch for the rigid airplane is, 

^eXAD ( 4 ’ .2 ^ a ^ 40 

f^HCXAD ( ®l40 ’ .2 ^ a > 40 


Cv • — 

^«o 


= < 


while the elastic increment to the angle-of-attack rate derivative for lift/pitch is. 


AC 


X& 


ni-a\i 
V 3 /o 


DCXADF'^*'''1.2 ’"P‘0 


0.0 


LQSE=.T. 
LQSE=.F. 


(5.23) 


(5.24) 


(5.24) 


(5.25) 
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5.4 Drag 


Drag is not discussed in MDC KV2A1 , vol. 1 for the maneuvering configuration. The total drae 
coefficient equation is taken from the MDC batch simulation and is repeated below, 

+ + +Co 5__8„, + Coj^5f 




'RT 


(5.26) 


5.4.1 Basic Drag Coefficient 
The basic drag coefficient is given by. 


^CDO 


( a(^j^2 ’ ) a<40' 


, 2 . 


Cd = 


'HCDO 


12 

1 90 


‘ .2 


l.l 

.2 


( «l40 . Ml V ) a > 40° 


(5.27) 


When in the low-angle-of-attack regime, an additional drag increment due to lift is calculated. 
However, does not represent zero-lift drag. Cq^ has the characteristic parabolic shape of the 

standard drag pol^and already includes the basic effects of lift on drag. An additional drae 
increment due to lift is defined in the MDC batch simulation as ^ 


AC 




^DCDCL ( ClL. 5 . M('2 ) a <. 40° 
0.0 a > 40° 


(5.28) 


small corrections to Cp^ and can be positive or negative. In the MDC batch 

simulation comments, is described as a correction "from tunnel drag to full-scale" and is 

4 ? explicit formulation of the equations of motion 

(section 4.3 of this document). Cl is not readily available in the flHhas simulation. Therefore, in 

the //^ftos simulation implementation of (5.28), AC^^ is a function of Cl instead of Cl- The 

difference in aircraft response induced by using Cl^^ instead of Cl in equation (5.28) is negligible. 

5.4.2 Stabilator Drag Increment 

The drag increments due to left and right stabilator (horizontal tail) deflection are given by: 


AC 




- < 


»H 


L.R 


f^DCDDH ( “1^12 ’ ^Hl - 2 / ’ ^1^2 ^ tt ^ 40° 


[*^HDCDDH ( ®1?() ' g| , Mf' 2 * ) a > 40° 


(5.29) 
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S.4.3 Leading-Edge Flap Drag Derivative 
The drag per degree of symmetric leading-edge flap is: 



fPcDDN(“l4 • ) 



^CXDN , M[ 


. 2 . 

.2 ■ 


a <. 40° 
a > 40° 


5.4.4 Trailing-Edge Flap Drag Derivative 
The drag per degree of symmetric trailing-edge flap is. 



/cDDF<«l^i. Ml‘2' ) 

0.0 


a ^ 40° 
a > 40° 


5.4.5 Speed Brake Drag Derivative 

The drag increment due to full speed brake deflection is. 






»SB 

1^60 


60 


DCDSB 


(all?, M[V ) 


1.2 


.2 


HDCDSB 


(a(^?,M[*,") 


1.2 


40 


.2 


a < 40° 
a > 40° 


5.4.6 Landing Gear Drag Increment 

The drag increment due to full landing gear deployment i.s. 


AC 


DlG 


f^DCDLG ^ ® ^ a ^40° 

0.0 a > 40° 


(5.30) 


(5.31) 


(5.32) 


(5.33) 


5.4.7 Rudder Drag Increment 

As discussed in the lift/pitch term build-up, a drag increment due to symmetric rudder (toe-in) has 
been added to the maneuvering configuration by using the increment developed for the high-lift 
configurations and discussed in section 4.1.10 of MDC A7247, vol. 1. Recall that the low 
speed/high lift configurations are considered low-angle-of-attack configurations. As discussed in 
section 4. 1.2.9, no elastic/rigid adjustments are defined for the rudder toe-in coefficients in MDC 
A7247, vol. 1. 
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0.0 


FcDDRT(“l"?0'«RTl’!’5.8f/Al“) “ 


a > 40* 


5.5 Side-Force, Roll, Yaw (Y,l,n) 

As in the case of lifl/pitch, the s^e reconsmiction equatiw 

yaw. The appropriate equation is recovered by substituting Y (side force), (roll) , o 
(yaw) for "X" in equation (5.34) below. 

Cxrf = Cx„ + AC^«P ] +ACxjj[^ ] +Cx^A8„ 


* CxA8f*«f ^ ''^X8 ^ 


+ '^*^X8x‘*Cx5^'*Cx8r ■" '*Cxa8H*^XA8H**” * 


+ ACx 


% *^CXp f Xr + ^^lr)[2W 


5. 5. 1 Basic Airframe Coefficient (Y,l,n) 

The basic airframe coefficient for the rigid airplane with controls at zero is given by. 


k Cx^ + (l-k)Cx^^ 


a <40° 


where. 


^cxoi ^ ® ^-10 ’ ^^.2 ^ 


^CX02 ( P ’ ® l-io ’ ^h.2 ^ 


M-l.l ,1 

and k = — j — Ip 

5.5.2 Leading-Edge Flap Effects (Y,l,n) 

If sideslin is nonzero, then symmetric leading-edge flap deployment pr^uces a change in^e 
la^dSSTco^fncienK. The incement .0 full 25 de^ deploymeu. of symmetnc 
leading-edge flaps in the presence of nonzero sideslip is given by. 


40 




^ "j^pDCXON ( “U ’ *^^.2 ^ 

. 0.0 


a <40° 
a > 40° 


(5.37) 


where 

reconstruction 
is given by: 


is to be multiplied by a normalized leading-edge flap position in the appropriate 
equation. The change in total coefficient per degree of differential leading-edge flap 



= 0.0 


(5.38.a) 



f^ClFDDN ( ^l’ ^^2 ^ a <_ 40 
0.0 a > 40° 


(5.38.b) 



0.0 


(5.38.C) 


where: 



LQSE=.T. 

LQSE=.F. 


(5.39) 


5.5.3 Trailing-Edge Flap Effects (Y,t,n) 

Increments due to full 20 degree deployment of symmetric trailing-edge flaps in the presence of 
nonzero sideslip are, 



^ ^DCXDF ( IPII() ’ 

.0.0 


a < 40° 
a > 40° 


(5.40) 


where as in the case of the leading-edge flaps, ACxg^ is multiplied by normalized trailing-edge flap 
position in the appropriate reconstruction equation, and the factor kg^ is defined as. 


, 1.1-M 


li 


(5.41) 


The change in total coefficient per degree of differential trailing-edge flap is given by, 

Cv^Sf = 0.0 (5.42.a) 
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f 



FciFDDF(«f>SL(r”-. 

P'cifddf(«14'Hsl=6o«)o,mi22) 


lo.o 


a < 40°; LQSE=.T. 

a < 40°; LQSE=.F. 
a > 40° 


(5.42.b) 



0.0 


5.5.4 Aileron Effects (Y,l,n) 

The increments due to left and right aileron deflection arc given by. 


AC 


X8 


Al.R 


^ *^DCXDAL ( *^Al,r1^25 ’ *^^^2 ^ ^ 

^ ^HDCXDA ( ®l40 ^1^25 ’ ^^.2 ^ ® ^ 


(5.42.C) 


(5.43) 


where k = 1 fw the left aileron and k = -1 for the right aileron. 
The elasticMgid ratio for aileron effectiveness is. 


rE = 

f^FRCYDA 


LQSE = .T. 

CV8A 

.1.63 


LQSE=.F. 

_ J 

f^FRCLDA 


LQSE=.T. 

V ■ 1 

.1.25 

LQSE = .F. 

- i 

[f^frcnda 


LQSE = .T. 

'^"Sa I 

..25 

LQSE = .F. 


(5.44.a) 


(5.44.b) 


(5.44.C) 


The elastic/rigid ratios assume values other than one in the rigid limit for the same reasons as 
discussed in section 5.3.3 of this report. The tHbular look-up wind-tunnel data was taken before a 
configuration change that increased aileron area. 


5.5.5 Rudder Effects (Y,l,n) 


Increments due to antisymmetric (both t.e.left) rudder deflection arc given by: 



I ~ Fdcxdr ( , |5r|[ 5,« , ) a < 40° 

f^HDCXDR ( “I40 ’ ^r) “ ^ “^0° 


(5.45) 
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The flexAigid lado far the effect of antisymmetric rudder is: 



, PpRCXDR ( ^1^2 ’ ) LQSE=.T. 

.1.0 LQSE=.F. 


(5.46) 


An additional effectiveness factor that incorporates the effect of symmetric stabilatm* on the 
increment due to antisymmetric rudder at high angle of attack is given by. 



*^HRXDRH^ “ ^40 • ^ ® ^ "X" = "n 

1.0 otherwise 


(5.47) 


5.5.6 Stahilator Effects (Y,l,n) 

The derivative with respect to a differential stabilator, which is itself a function of angle of attack, 
symmetric stabilator, and Mach, is given by. 



^CXDDH ( ’ ^1^2 ) ® ^ '^O® 

f^HCXDDH ( ®l40 » ^f .2 ^ ® ^ 

k. 


(5.48) 


The elasdcAigid ratio on the differential tail derivative is defined as. 


c u 1 60000^ 

*^FRXDDH^‘'^1.2 ’ ”pl 0 ^ 


LQSE=.T.; "X"=’’Y"or"r 


“cxaSh ' 


^FRXDDh(^[^ 2 ’ o”) EQSE=.T.; "X"=’’n” 

1.0 LQSE=.F. 


r 60000 _r40x 


5.5.7 Dynamic Derivatives (Y,l,n) 

The derivatives due to nondimensionalized roll rate are defined as. 



"Fcxp («[4 , M[^2 ’ Hpl 0 ®^^’®) a <40°; "X’’="r; LQSE=.T. 
Fcxp(a (4 , Mf ^2 > Hp=60000) a <40°; "X"=’T';LQSE=.F. 
Fcxp(a |4 , M[^ 2 ) a <40°; "X"="Y"or”n" 

/hcXp(«I 42’M[‘2^) a > 40° 


(5.49) 


(5.50) 


while the effect of nondimensional yaw rate is given by. 


43 


(5.51) 


C 


Xr 


^CXR ( ® [ -4 ’ .2 ^ ® ^ 

f^HCXR ( ®l40 ’ ® ^ 


The increment to yaw damping due to elastic deformation is given by, 

[f^FDCXRF ( *^1^2 ’ ) LQSE=.T. 


ACx = 

Xr 


(5.52) 


lO.O 


LQSE=.F. 


A separate flex/rigid ratio on roll damping is calculated for side-force and yawing moment. The 
flexible effect (Mt roll moment due to roll rate is included in the function used to calculate Cj^. 


c ( x/ir2 H 1 60000 

*"FRCYP ^ ‘'^1.2 ’ "pl 0 


'cy. 


lO.O 


) LQSE = .T. 
LQSE=.F. 


R" = ^ 

Cnp 


FfRCNp( Ml^2*“lo‘) LQSE=.T. 


(M=.2,a|p^) LQSE = .F. 


FRCNP 


(5.53) 


(5.54) 


Rp = 1.0 all conditions 
*P 


5.5.# Speed Brake Effects (Y,l,n) 


(5.55) 


The speed brake increments shown below are due to full, normalized, speed brake deployment and 
arc multiplied by sideslip angle in the reconstruction equations. 


60 J 


*^DCXSB ’iPlto^ .M[ .2) 

^HDCXSB^lPlIo^ ’^1.2 ^ ® ^ 


r2.v 


a < 40' 




(45-a\ 

1 

fSsa) 

[ 5 / 

0 

[60 j 


(5.56) 


V. 


S.S.9 Effect of Airplane Asymmetries (Y,l,n) 

A rolling moment is induced at zero sideslip and symmetric controls for angle of attack between 
45° and 60°. This rolling moment is attributed (section 4.1.2 of fMDC A7247, vol. 2]) to 
asymmetries in the nose. The asymmetric increment is given by. 
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AC, 


^asm 


^ 0.0 

^ PhCIASM ( ® I 40 ^ 

I 0.0 

lo.o 


a <40° 

a > 40°; "X" = ’T’;LRTE=.F. 

a > 40°; "X" = ’T';LRTE=.T. 
a > 40°; "X" = "Y"or”n" 


(5.57) 


5.6 Force and Moment Calculations 

In this section wc describe how the total aerodynamic forces and moments at the c.g. are 

The aerodynamic model has been worked up about an aerodynamic reference center. 
While it is common practice to transform the outputs of the aerodynamic model (force and moment 
coefficients) from the aerodynamic reference center to the c.g., it is less common to tr^sform the 
inputs to the aerodynamic model. Both transformations are performed in the /78has simulation in 
the interest of completeness. The practice of transforming the inputs has been ^n by the authors 
in a simulation model for the X-29 and in section II of Gainer and Hoffman (1972). 

First the velocity of the air is subtracted from the inertial translational velocity of the airplane to get 
relative velocity (with respect to the air mass) at the c.g. The body frame components of velocity 
with respect to the air mass are given by, 



(5.58) 


where [Uwg v^g w^g]' is the velocity of the air due to winds, turbulence, and gusts. 
The velocity components relative to the air mass at the a.r.c. are given by. 


%1 

fUal 

■ 0 -r q‘ 


''arf 

^ ^ Va U 

r 0 -p 

]yarf^ 

WarfJ 

[wj 

. -q p OJ 

IzarfJ 


(5.59) 


where [xarf Yarf ^arf]' locates the a.r.c. in the c.g.-centered body frame. Rotational velocities at 
the c.g. can induce small changes in translational velocities at the a.r.c. Incidence angles ( 0 ^ and 
pjf), airspeed, and Mach number, are then computed from relative velocities at the a.r.c. 

ttrf = arctan( Waj.f / Uaj.j. ) ^ 


Prf = arcsin( ) 


Vrf = '^ H 


2 2 2 

arf ''an- 


M,r = V,f/a 




(5.60) 
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Dynamic pressure (q) is left referenced to the c.g. The dynamic pressure in lbs/ft2 is given by, 

q = (5.61) 


where p is the density in slugs/ft3. A density altitude, Hp. is calculated from the air density, p. 

Hp is the altitude in feet that is appropriate for p in the 1962 Standard Atmosphere tables. Air 
density can be varied as an independent variable if required (e.g., hot day conditions). 

5. 6. 1 Steady Flow Outputs 

The quantities defined in equations (5.60) and (5.61) are then combined with the body frame 
rotational rates (p,q,r), pressure altitude, the control surface deflections ({5}), and the reference 

dimensions (c,b), to serve as inputs to the aerodynamic model described in sections 5.2 through 
5.5. Outputs of the aerodynamic model include the force and moment coefficients for steady flow 
conditions. The process is characterized by equation 5.62. 


^IlsF Aero( «rf . Prf .^rf , Mrf, q, Hp, (5}, p, q, r, c ,b) (5.62) 


wlrere "n" = "L”, "D", "Y", "1", "m", or "n". 
axis components via, 


The lift and drag coefficients are converted to body 


Cxjjf = - cos(arf)CD^j. + sin(arf)CL^^ 
Czsf = - s*n(arf)Co^j. - cos(arf)CLj.|: ^ 
The steady flow aerodynamic forces in the body frame are given by. 



|V)=o 


rcxsf 

l^ZsfJ , l^Zsf, 


The steady flow aerodynamic moments in the body frame are given by, 


(5.64) 


(M) I . 

|V) = 0 


M 

= ■1 M 

M 


Xsf 

Yj.f> qS^ 
Zsf 


b 


® "^aif Varf 



c ^mjjf 

b ^nsf> 

> + 

Zaif 0 -Xarf 

- -Yarf Xarf 0- 


Fy f 

1-3 


(5.65) 


As a simple sanity check, observe that a positive implies the a.r.c. is ahead of the c.g., and a 
positive force applied at the a.r.c. induces a negative pitch moment. My. 
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5.6.1 Unsteady Flow Outputs 


Wc seek now to define the 6x6 apparent mass matrix defined in equation (4.1 3) (repeated below) 
that arises from the Cl^ and stability derivatives. 

a {F) . 

[Ml. (4-13) 

a [V) 


First, in a process similar to equation (5.63), transform the a stability derivatives from stability to 
body frame components at the a.r.c.. 


= * cos(arf)Ci)^ sin(arf)CL^ 
Cz^ = - sin(arf)CD^ - cos(arf)CL^ 


(5.66) 


In iheflSbas simulation, as is typically the case, is zero. Define some intermediate terms that 
have units of mass and mass-length. 


m, = q S 


f - \ 
c 


l2VrfJ 

2 2 
“arf ^ 

J 


.slugs 


m 

u 

= 

slugs 

m 

w 


slugs ^ 



slug-ft 


(5.67) 


(5.68) 


For notational convenience, let V and (o denote the 3x1 column vectors associated with 
translational and rotational velocity, respectively, 



The apparent mass matrix is given by, 


IMIa 


a I FI A 

a IV) 


[Myvl [Mvo)l 
I M(fl V 1 f _ 


(5.69) 


(5.70) 


where. 
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A 

Wrf = 


(5.71) 


[Mvvl = 


m (vrf 0 -m 
u u 


0 0 .slugs 


m Wrf 0 -m 


[Mvcol = 


-myyarf mu(^arf+<Vri2aif) my^riyarf 


-m^yarf 


m^Wrfyarf 


-m^yarf 


,slug-ft 


fMtovl = (^^uZarr^-m^Xar^^-m Jwrf 0 ,slug-ft (5.73) 


-niu^rfyarf 


muyarf 


(5.71) 


(5.72) 




Afy^ 


(myZ.^rm^x^m^)(xarffWrfZarf) (-m^Zarf+m^Xarrmq)^rfyarf 
m^y^ -mu(xarr+Wri^„f)yarf 




,slug-ft7 


(5.74) 
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5.7 Implementation of the aerodynamic model 

As described in section S.l, the aerodynamic fence and moment coefficients are calculated in 
subroutines SFAERRF and USAER^. In sections 5.7.1 and 5.7.2 the input/output variables of 
the subroutines SFAERRF and USAERF, respectively, are described. 

5.7.1 SFAERRF - Input/Output List 

The subroutine SFAERRF produces the total force and moment coefficients about the aerodynamic 
reference center for steady flow conditions (constant angle of attack and sideslip). The names of 
each input/ou^ut parameter as they appe^ in both the subroutine (SFAERRF) and in the fISbas.csl 
code which SFAERRF is called, are listed in tables 5.3 and 5.4. The density altitude is used 
as an input to the flexAigid multipliers and increments. Ground clearance was used for ground 
effect c^culations. The calculations associated with ground effect have been removed, but the 
parameter list was left unchanged in the interest of compatibility. 




Symbol 

Definition 

SFAERRF 

.csl 

a 

angle of attack, degrees 

ALPDEG 

ALDGRF 

p 

sideslip angle, degrees 

BETDEG 

BEDGRF 

M 

Mach number 

MACH 

MACHRF 

P 

Body frame roll rate, rad/sec 

P 

P 

q 

Body frame pitch rate, rad/sec 

Q 

Q 

r 

Body frame yaw rate, rad/sec 

R 

R 

Hp 

Altitude above sea-level in a standard atmosphere 
equivalent to a given density, feet 

ALT 

HRF 

Hgci 

(^und clearance, feet 

HGCL 

HGC 

q 

Dynamic pressure, Ibs/ft^ 

QBAR 

QBAR 

C 

Reference aerodynamic chord, feet 

CREF 

CWRF 

b 

Reference aerodynamic span, feet 

BREF 

BWRF 

V 

Total airspeed at the aero reference center, ft/sec 

VTRF 

VTRF 

Shr 

Right stabilator (horizontal tail), positive trailing edge 
down (+ t.e.d), degrees 

DHTR 

DSR 


Left stabilator (horizontal tail), + t.e.d, degrees 

DHTL 

DSL 

5ar 

Right aileron, + t.e.d, degrees 

DAR 

DAR 

5al 

Left aileron, + t.e.d, degrees 

DAL 

DAL 

^Rr 

Right rudder, + t.e. left, degrees 

DRR 

DRR 

Srl 

Left rudder, + t.e. left, degrees 

DRL 

-DRL 

^^R 

Right trailing-edge flap, + t.e.d, degrees 

DTFR 

DFR 

K 

Left trailing-edge flap, + t.e.d, degrees 

DTFL 

DFL 

SnR 

Right leading-edge flap, + l.e.d, degrees 

DLFR 

DNR 

Sol 

Left leading-edge flap, + l.e.d, degrees 

DLFL 

DNL 
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Table 5.3. Input parameters for subroutine SFAERRF (continued) 


Symbol 

l^«nltlon' 

SFAERRF 

.csl 

^SB 

Speed brakes, 0-60, degrees 

DSBK 

DSB 

^LG 

Normalized landing gear, 0*1, full down 

DLG 

DLG 

logical 

(.T.) 

E)etailed print of intermediate variables on 
unit=20. 

LX>EBUG 

LMBUG 


(.F.) 

No action 



logical 

(.T.) 

RexAigid multipliers and increments calculated 

LQSE 

LQSE 

logical 

(.F.) 

Aero for rigid aiiplane, flex/rigid -> 1 . 

(.T.) 

Real time equivalent aero 

LRTE 

LRTE 


(.F.) 

Enhanced aerodynamic model 




Table 5.4. Output parameters Tot subroutine SFAERRF 


Symbol 

Definition 

SFAERRF 

.csl 


Total drag coefficient at a.r.c. evaluated at ( 6t= 0) 

CDREFO 

CDRFSF 

^Y.sf 

Total side-force coefficient at a.r.c. evaluated at ( a= 0) 

CYREFO 

CYRFSF 


Total lift coefficient at a.r.c. evaluated at ( o?= 0) 

CLREFO 

CLRFSF 

^Isf 

Total roll coefficient at a.r.c. evaluated at ( o= 0) 

CIREFO 

CIRFSF 


Total pitch coefficient at a.r.c. evaluated at ( a= 0) 

CMREFO 

CMRFSF 


Total yaw coefficient at a.r.c. evaluated at ( a= 0) 

CNREFO 

CNRFSF 
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S.7.2 USAERRF - InputlOutput List 


The subroutine USAERRF returns the lift and pitch derivatives with respect to the nondimensional 
rate-of-change of angle of attack. The names of each input/output parameter as they appear in both 
the subroutine (USAERRF) and in the flSbas.csl code from which USAERRF is called, are listed in 
tables S.S and 5.6. 


'^bl^^SjIngutjgarameto^^o^ubr^^ 



Definition 

SFAERRF 

.csl 

a 

angle of attack, degrees 

ALPDEG 

ALDGRF 

M 

Hp 

Mach number 

Altitude above sea-level in a standard atmosphere 

MACH 

MACHRF 

logical 

equivalent to a given density, feet 

(.T.) Detailed print of intermediate variables on 

ALT 

HRF 

unit=20. 
(.F.) No action 

LDEBUG 

LCCBUG 

logical 

(.T. ) Flex/rigid multipliers and increments calculated 
(.F.) Aero for rigid airplane, flex/rigid -> 1 . 

LQSE 

LQSE 


Tabl^^6^0utpu^»aramete^^o^ub^^ 


Symbol 

Definition 

SFAERRF 

.CSl 



CLADRF 

CLADRF 

Cma 


CMADRF 

CMADRF 




6. ENGINE THRUST MODEL 


The engine model takes inputs from the pilot throttle and an autothrottle command, uses current air 
data (altitude, Mach number, and dynamic pressure), and computes the engine forces acting on the 
airframe for the equations of motion described in section 4. This engine model is referred to as the 
"(3T Half’ model in MDC documentation and correspondence (MDC internal memorandum 338- 
5250, dated 5-Dec-80, from R.O. Michael and B.R. Williams), and apparently is a refined model 
(circa 1980) based on installed engine performance flight tests. It should be noted that the authors 
focused on implemendng this model as provided, in order to facilitate comparison with an existing 
re^-time simulation model. The authors are not able to explain some aspects of the architecture of 
this model -- these will be noted as appropriate in the text below. 

The model from pilot/autopilot inputs to thrust detennination for a given engine is represented 
graphically in figure 6.0. Since the airplane has two engines, this model is implemented in code by 
parallel sets of ACSL equations and external subroutine calls for the left and right engines. 
Although the autopilot input is provided for future use, in the current simulation this input, CAT, is 
set to zero. Therefore, PLA, the total throttle command, is equal to CP, the pilot throttle 
command. 



Figure 6.0. Engine Model Block Diagram 


6.1 Throttle Dynamics 

PLA is input to the Throttle Dynamics, a straightforward .servo model with throttle position, DPL, 
as its output. This servo model includes a limited integrator, which functions as a normal 
integrator as long as its output remains between the indicated limits {31, 130}. If the integrator 
output equals or exceeds these limits, the state derivative DPLD is set to zero within the limited 
integrator macro and the integrator output is appropriately limited. The limited integrator block 
shown in figure 6.1 is the implementation provided by ACSL’s built-in limited integrator 
command. This structure is shown for graphical convenience; in the ACSL simulation model, a 
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custom macro using IF condition tests is actually used. Either implementation should be 
equivalent; the custom macro was developed by the authors to retain control over the variable 
names of states and state derivatives [ACSL (1987)]. 



Figure 6.1. Throttle Dynamics 


6.2 Static Performance 

As shown in figure 6.2, throttle position (DPL), along with altitude (H), Mach number (MACH), 
and an internal engine dynamics state (VAEO) are input to the static peiformance module (Fortran 
subroutine ENGl). VAEO is an artificial state created to break an algebraic lo(^ fw the static 
performance motile. In the original specification, VAE is an input into the static performance 
module, and VAE is dependent on output THRC. VAEO is created by passing VAE through a fast 
first order lag and is discussed further in section 6.5, Core Dynanucs. 

The outputs of the static performance module are engine thrust command (THRC), fuel flow 
(V^. idle thrust (TIDLE), maximum military power thrust (TM^), minimum afterburner thrust 
(TAUGMN), maximum afterburner thrust (T AUG MX), and various gains and parameters for the 
engine dynamics (TESTV, GKl, PI, P2, TAU). THRC and WF are each computed as a table- 
lookup function, of H, MACH, and DPL. There are two table-lookup functions for THRC; one 
f(M* military power throttle settings and one for afterburner throttle settings. 



■■■■■■■ 





Static Performance 





Core Thnist Conunand 


npi 

' W. 


A 


LfrLf 





— ^ THRC 

MACH 


Table 




H 


Look-Ups 

Augmented DPL >106.5 





Thrust 

GKl. TAU, 




Command 

^ P1.P2TMIL 




(includes afterburner) 

TESTV 




TAUGMN 


■■■■■■■■ 



Figure 6.2. Static Engine Performance 


The authors cannot fully explain the internal architecture of the static performance module — the 
switch between Ae Core and the Augmented (afterburner) Thrust Command lookup tables occurs 
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at a pPL (PLA after dynamics) of 106.5*, rather than the expected military/afierbumer throttle 
position boundary of 88-9 1 * (PLA). ^ 

npLE, TMIL, TAUGMN, and TAUGMX are computed for appropriate fixed throttle positions 
using the same table-lookup function as for THRC. TESTV, GKl, PI, P2, TAU are computed as 

table-lookup functions of ^^ . 

THRCis used by two sets of engine dynamics running in parallel: Afterburner Dynamics and 
Core L^namics. This will be described further below. GKl. TAU, TESTV. and TMIL are used 

TAUGMN is used in the Discrete Event Logic (see below). 
Wh, I IDLE, and TAUGMX are not currcndy used elsewhere in the ACSL simulation model. 

6.3 Discrete Event Logic 

As dcsOTbed alwve, the engine thrust command (THRC) is provided to two sets of engine 
dynamics running in parallel: Afterburner Dynamics and Core Dynamics. These dynamics 
provide a single thrust output that is used as the total net thrust for a given engine. The Discrete 
Event Logic (figure 6.3) is used to determine when to switch from one set of dynamics to the 
other, and implements code: 

1) to select the appropriate thrust output from eithtr Afterburner Dynamics or Core 
Dynamics; 

2) to enable a smooth transition between the two dynamic systems by updating 
appropriate state variables - this is accomplished in ACSL using DISCRETE 
blocks (one per engine) called by the ACSL SCHEDULE operator [ACSL (1987)]. 

The authws n^ that the Discrete Event Logic uses a thrust comparison test that is inconsistent 
with the ^itch test used in the static performance module to select the output of the thrust table 
lookup. The authors cannot explain why this is so; nor can we explain why this parallel dynamics 
^hitecture wm cht^n. It would seem more physically intuitive to implement an Afterbi^ner 
Dyn^cs module that outputs an incremental thrust to be added to the Core Dynamics thrust to 
yield total net thrust. 
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6.4 Afterburner Dynamics 

As shown in figure 6.4, THRC is one of three inputs to the Afterburner Dynamics. The physics of 
the MDC model are such that when the thrust command (THRC) exceeds the afterburner thrust 
(TAB), TAB integrates up at the rate ABUL. Conversely, when TAB is greater than THRC, TAB 
integrates down at the rate ABLL. 

ABUL and ABLL are computed based on current values of altitude (H) and dynamic pressure 
(QBAR). The lower rate limit, ABLL, is computed as a function of H, 

ABLL = .0275 H - 5500 (6.1) 

The uppo* rate limit, ABUL, is computed in two steps. The first step is to compute an intermediate 
value as a function of QBAR, 

ABUL* = 1.7QBAR2 - 12.8 QBAR + 3000 
ABUL is then computed by limiting ABUL* as, 

r 900 ABUL* < 900 

ABUL = 5 ABUL* 900 < ABUL* < 4000 
I 4000 4000 < ABUL* 

In the ACSL implementation of this system, the error signal THRC minus TAB is multiplied by a 
large gain (200) to yield ZZAB. ZZAB is pas.sed through a rate limiter and is appropriately limited 
to ABUL or ABLL to yield RTAB. RTAB is then integrated to yield TAB. Synchronization of the 
Afterburner Dynamics with the Core Dynanucs is achieved by updating TAB with the last value of 
core engine thnist, TT, when transition from Core Dynamics to Afterburner Dynamics occurs. 

This has the effect of "restarting" the afterburner dynamics from the thrust value at transition. 


(6.2) 


(6.3) 



Figure 6.4. Afterburner Dynamics 
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6.5 Core Dynamics 

As shown in Hgure 6.5, THRC is one of two primaiy inputs to the Core Dynamics, which 
implements a second order dynamics approximation to a more complicated, higher order dynamic 
model of the turbomachinery used for batch analysis of the engine by its manufacturer. The other 
five inputs are variables that have been described previously: TESTV, GKl, PI, P2, TAU. 

The error signal VAC is computed by subtracting VAF from THRC. VAC is multiplied by the 
gain GKl and is passed through the rate limiter bounded by limits RTEUL and RTELL to yield 
VAE. The rate limits RTEUL and RTELL are the engine spool-up/-down rate limits and are 
defined as CONSTANT variables in the ACSL simulation model. Therefore, although they have a 
"default" value, these limits can be altered at runtime (ACSL (1987)] by the user. The default 
values arc ±5500. 


Core Dynamics 


THRC— ( 


RTEUL 


RTELL 


.02 S + 1 


VAF 1 



Reset VAF=TAB 
if AB->Core 

Reset 

Vm^VAH 

if 

AB->Core 


AB->Core 


VH1 1 
S 


0.0 Reset VIIDI^Oij 
AB->Core AB->Core 

VI 


I ^ 2°^ Order 
VHID Dynamics 


1 

t 

i 


VAEO 


THRR>0, 


< TESTV 


PI.P2 


THRRO 


VAII remains constant until updated by TMIL 
VAF when conditions test true. 



Figure 6.5. Core Dynamics 


VAE is input to a low-pass filter to generate the artindal state VAEO as described earlier to break 
the algebraic loop for die static performance module. The time constant for this filter is sufficiently 
small (50 rad/sec) to insure that all engine dynamics are unaffected by the filter. VAE is input to an 
integrator, yielding VAF. Synchronization of this integrator with the Afterburner Dynarrucs is 
achieved by updating VAF with the last value of afterburner engine thrust, TAB, when transition 
occurs from Afterburner Dynamics to Core Dynamics. This has the effect of "restarting" this 
integrator with the thrust value at transition. 
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VAF is the majOT component of the core engine thnist, represented by a rate-limited first order lag. 
VAF is used to form the error signal VAC as previously described. VAF is tested against 
TAUGMN in Discrete Event Logic to determine whether to initiate a transition between Core 
Dynamics and Afterburner Dynamics. VAF is supplied to the "Nonlinear Switch" [MDC internal 
memo 338-5250]. The "Nonlinear Switch" will te described below. VAF is summed with VAK 
to determine core engine thrust, TT. 

The "Nonlinear Switch" can only be active when the core dynamics are being used. Inputs to this 
switch are VAF, VAC, TMIL and TESTV. At each simulation time step for which the condition 


abs(^)<testv 


(6.4) 


is true, the "Nonlinear Switch" assigns VAH to the current value of VAF. Since the "Nonlinear 
Switch" is the only mechanism for defining VAH, it remains constant when the above condition is 
not true. TMIL is a function of Mach number and altitude; TESTV is a positive-valued function of 
VA^TMIL, varying proportionally with the magnitude of VAE. Therefore VAH only changes 
when VAC is large relative to VAE (big change in throttle command). 

Subroutine ENG2, coded in FORTRAN, is only called when the core dynamics are being used. 

VAF 

The input is an internal engine dynamics parameter The output is a gain for the second 

order dynamics block (GK3). 


VHl is subtracted from VAH and the result is divided by TAU, yielding VHID. VHID is 
integrated to give VHl. Synchronization of this integrator with the Afterburner Dynamics is 
achieved by setting VHl equal to VAH when transition occurs from Afterburner Dynamics to Core 
Dynamics. By "zeroing" the input signal to this integrator, this part of the dynamics is not active at 
transition. 


VHID is input to the 2^ Order Dynamics block, a transfer function with the following Laplace 
expression. 


VAK ^ (P2)2(GK3) 

VHID s2 + 2(P1)(P2)S + (P2)2(GK3) 

The VAK term provides the small "ballooning" effect for the engine described in MDC internal 
memo 338-5250 (1980). VAK is summed with VAF to give the total thrust for the Core 
Dynamics, TT. 

6.6 Total Engine Thrust 

The final total thrust for each engine, TH, is determined by the position of the AB/Core switch. If 
THRF ^ 0, then TH is equal to the output of the Core Dynamics, TT; otherwise, TH is equal to the 
output of the Afterburner Dynamics, TAB. 

6.7 Calculations for Equations of Motion 

The equations of motion for the simulation, described in section 4, require the 3-axis body frame 
components of force and moment due to the thrust from the engines. Therefore, operations are 
performed on the final total thrust for each engine, TH, and then added together to yield these 
required force/moment components. Note that engine gyroscopic moment effects are not modeled. 
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In keeping with overall assumptions that an x-z plane of symmetry exists, the engines arc assumed 
to be located symmetrically with respect to the body frame x-z plane. However, the thrust axis of 
each engine is not parallel to the x-axis. Each engine is canted out slightly (1.98*) so that the exit 
nozzles arc closer to the airplane centerline than the front of the engines. The pitch or elevation 

angle of each engine is zero. Engine cant is measured by an azimuth angle ('I^eng) represents 
rotation about the body frame z-axis. Given this information, the body-axis components of thrust 
for each engine arc computed as. 


TH,;L 

— THl 

'FHyjL - THl sin(Tgng^) 



TH,;L 

= 0 

TH^.r = THr COs('FengR) 

► 

(6.6) 

THy;R 

= TH|^ 

TH,.r = 0 




where the subscripts "L" and "R" denote left and right engines, and where. 




and 'FengR= 1-98* 


(6.7) 


The force quantities for the equations of motion can now be computed. 


Fxe = + TH,;R 

^yE ^ + '^^y,R 

Fze = + TH,.r 


(6.8) 


The engine thrust is assumed to be applied where the thrust axis intersects the "engine/airftame 
interface”, defined in the engineering drawing, figure 3.1. As can be seen from figure 3. 1 , the 
location of the "engine/airframe interface" for the right engine is defined in the FS/BITWL 
coordinate frame as (687.5, 18.9, 1(X).0). The interface point for the left engine is (687.5, -18.9, 
100.0). This information is specified in the variables FSENG, BLENGR, BLENGL, and WLENG in 
CXINSTANT statements in the fIShas simulation. The vector that locates the engine/airframe 
interface point in the c.g.-centered body frame is computed in units of feet for use in calculating the 
moments about the c.g.. 


xerf 

ycrfR 

yerft 

Zerf 


(FSeng ~ FScg) / 12 
(BLengR ~ BL^g) / 12 

(BLengE “ BLcg) / 12 


/-\x/i - wi 


(6.9) 


The calculations of (6.9) arc performed once in the initialization phase of the simulation. Given 
this information, the moment quantities for the equations of motion can now be computed. 


Mxe = - ZerfFy^ + ycrfRTlly;R + ycrfLTHz;L ^ 

Myg = Zerf — Xgrf Fjje r (6.10) 

Mze = -(ycrfR THx;R + ycrfL TH^-l) + Xerf Fy^J 


58 



7. SENSORS 


The sensor code was taken from the LaRC real-time simulation {dmsflS). These models can also 
be found in subroutine FCSENS of the MDC simulation of the F/A-18. A table of sensor 
characteristics can be found in figure 7-7 of (MDC A7813, vol I, rev A, page 7-16] and selected 
parts arc repeated below as Table 7.1. In the dmsfIS simulation, the senstn* calculations and 
dynamics arc performed in the same subroutines that model the control laws. In keeping with a 
policy of modeling all continuous dynamics in the derivative portion of the ACSL simulation, 
sensor MAQiO's were written to implement the sen.sor models, and the sensor code was stripped 
out of the cmtrol law subroutines that are called from the FCS discrete block. 


Table 7. l.a . Rate Gyro and Accelerometer Parameter Bounds 


Parameter 

PitchA'aw 
Rate Gyro 

Roll Rate Gyro 

Normal 

Accelerometer 

Lateral 

Accelerometer 

Range 

t 60 7sec 

t 3(X) Vsec 

t lOg 

t2g 

Threshold 

.005 Vsec 

.01 Vsec 

.002 g 

.002 g 

Linearity 

1 % full scale 

1 % full scale 

1 % full scale 

1 % full scale 

Null 

.5 Vsec 

1 .5 Vsec 

.05 g 

.01 g 

Hysteresis 

.05 Vsec 

.30 Vsec 

.005 g 

.001 g 

Gradient Tolerance 

t 5 % 

t 5% 

+ 5% 

t5% 


Table 7. 1 .b. Pedal Force and Pressure Measurement Parameter Bounds 


Parameter 

Pedal Tran.sducer 

Static Pressure 

Dynamic pressure 

Range 

1 134 lbs 

10,38] in HG 

f0,55] in HG 

Threshold 

.5 lbs 

.0005 in HG 

.001 in HG 

Linearity 

2 % full scale 

inci in gradient 
tolerance 

incl in gradient 
tolerance 

Null 

0.5 % full scale 

0.09 in HG 

0.105 in HG 

Hysteresis 

1.0 lb (0 to 60 lbs) 

4.0 lb (60 to 134 lb.s) 

incl in gradient 
tolerance 

incl in gradient 
tolerance 

Gradient Tolerance 

t5% 

1 0.3 % 

+ 0.3 % 
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The sensor outputs required by the subroutines that model the flight control laws are given in Table 
7.2. The pedal force is an autonomous input to the simulation. 


Tabl^j2^ensorOmguts 


Symbol 

Variable 

Name 

Definition 

«s 

ALFS 

angle of attack, degrees 

"ys 

AYS 

lateral acceleration, positive right, g’s 

"zs 

AZS 

vertical acceleration, positive up, g's 

Ps 

PS 

body axis roll rate, ’/sec 

qs 

QS 

body axis pitch rate, ’/sec 


RS 

body axis yaw rate, ’/sec 


QCIS 

compressible impact pressure, PSF 


PSTS 

static pressure, PSF 


The sensor values used by the flight control system model can be either the ideal values or values 
computed using various sensor models. The variable ISENS determines which type of sensor data 
is sent to the flight control subroutines. The ISENS variable was made integer type in order to 
accommodate multiple sensor models. The simulation currently calculates two sets of sensor 
outputs. The first set consists of the variables created by adding a "0" to the variable names in 
table 7.2, i.e. (ALFSO, AYSO, ..., PSTSO), and represents ideal sensors. The second set 
consists of the variables created by adding a "1" to the list in Table 7.2, i.e. {ALFSl, AYSl, ..., 
PSTSl ), and represents the sensor models used in the MDC simulation. At one time, a third, 
more ctnnplex, set of sensor models was implemented in the AGCB simulation to produce the 
variable set { ALFS2, AYS2, ..., PSTS2). The 1SENS=2 option was removed to reduce 
complexity and because it was not entirely integrated with the simulation trimming algorithm, but 
could be reinstalled if the need arose. The final output variables in Table 7.2 are determined from 
the sensor model sets according to the value of the integer switch, ISENS, as shown in figure 7.1 . 


IF (ISENS=0) THEN 
ALPS = ALFSO 
! 

PSTS = PSTSO 
ELSE 

ALPS = ALFSl 

# 

PSTS = PSTSl 
END IF 


Figure 7.1. Code block that selects sensor model 
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The sensor equations appear in both the initial section of the program code and in the derivative 
section. The initial section contains full code for the computations. The derivative section makes 
use of macros which appear at the top of the ACSL program as ACCELl, ALGBTl, PITOT, and 
RGYROl. 

Ideal static pressure (PSTATC) is computed from the aircraft altitude (H, Height, in the main loop, 
HIC, Height Initial Condition, in the initial section) through the subroutine call to PSTAT. It 
inteipolaies for the static pressure value from an single independent variable (H) ATMAT62 data 
table. Ideal compressible impact pressure (QC) is computed from the ideal static pressure 
(PSTATC) and the Mach number (MACH, for main loop, MACHTR, MACH TRim, for the initial 
section) through the subroutine CPRESS. It interpolates for the compressible impact pressure 
from a two independent variable (MACH and PSTATC) data table. 

Sensed static pressure (PSTSl, ISENS = 1) is computed from the ideal pressure divided by a 
scaling factor, SPSFXD. SPSFXD is computed as a function of Mach number and sensed angle 
of attack (ALFSl). The QCDA,C?CDB...CiCDH variables which appear in the scaling equation are 
a function of Mach number and are defined as constants in the initial section of the program code. 
The sensed compressible pressure (QCISl, ISENS=1) is also determined as a function of Mach 
and angle of attack, where the static pressure is added to the ideal compressible pressure before the 
computation and then removed afterwards. A final computation is then made based on a limited 
angle of attack. These computations model the Pitot tube and static probe mechanics, where the 
measured value, as determined by pressure differentials, does not exactly reflect the true air 
pressure values. These pressure values are used in the flight control system for air data 
scheduling. The air data scheduled gain values were designed to be scheduled from measured 
pressure values and not the true or ideal pressures. 

Sensed angle of attack is computed from the true angle of attack with a gain and bias based on air 
flow. The bias, (ALFB), is determined by the landing gear position, if the landing gear is retracted 
(the normal position for any runs) ALFB=1.9, extended ALFB=1.2. The gain, GKAOA, (where 
it is divided into the true alpha minus the bias) is a function of the Mach number. The measured 
angle of attack is limited to -14 and 56 degrees. 

7.1 Ideal Static Pressure and Compressible Impact Pressure 

Static pressure, Pj, is calculated as a function of altitude in feet by the subroutine PSTAT. The 
function. 


PS, Ibs/ft^ = Fpsxat(H) 

implemented in subroutine PSTAT is plotted in figure 7.2. 
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Altitude (feet x 1000) 

Figure 7.2. Static Pressure as a Function of Altitude 

Subroutine CPRESS calculates compressible dynamic pressure as a function of Mach and static 
pressure, Pg, as follows: 

Qc = (Qc/Ps)Ps 

where the ratio Qq/Ps » function of Mach number and is given by figure 7.3. The subroutines 
PSTAT and CPRESS are called from the DERIVATIVE block. 



Figure 7.3. Ratio of Compressible to Static Pressure as Function of Mach 
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7.2 Ideal Sensor Set (ISENS=0) 


®s0 degrees 

ALFSO 

"yso "y 

AYSO 

"zso = - "z - 1 . g's 

AZSO 

PsO = (180/tt)p, Vsec 

PSO 

qsO = (180/7c)q, Vsec 

QSO 

’’sO = (180/7t)r, Vsec 

RSO 

Qc ,0 = ‘3c, PSF 

QCISO 

Ps,() = Ps. PSF 

PSTSO 


where, 

Py = {v+ru-pw-glyj,J/g 
Hz = {w + pv-qu-gl^zl/g 


7.3 Accelerometers (ISENS=1) 

The calculations for the lateral and vertical accelerometers are performed in the macro ACCEL 1. 
The lateral accelerometer output, AYSl, is positive for an acceleration to the right and the vertical 
accelerometer output, AZS 1 , is positive for an upward acceleration. Both AYSl and AZS 1 have 
units of g's. Lever ^ effects are calculated and the 1 g bias removed from AZSL Both outputs 
are bounded according to Table 7. 1 .a. Hysteresis, bias, quantization and misalignment are not 
modeled in the MDC sensor set. The MDC simulation, mdcflS, has comments that refer to a 

transfer function of G(s) = (^220)+l approximates the accelerometer dynamics, but was 

deleted from the MDC simulation. The flSbas simulation does not model accelerometer dynamics 
either. The location of the accelerometers in the body frame in feet is given by, 

xa = -(FSa-FScg)/12. 

ya = (BLa-BLcg)/12. 

za =-(WLa-WLcg)/12. 
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where Xa, ya, Za represent distances in feet forward from the c.g., right of the c.g. and below the 
c.g., respectively, and where, 


FScg 

= FSCG, Fuselage Station c.g. 

= 455.0 inches 

Bl^g 

= BLCG, Buttock Line c.g. 

= 0.0 inches 

Wl^g 

= WLCG, Water Line c.g. 

= 102.8 inches 

FSa 

= NYLOC(l) 

= 306.75 inches 

BLa 

= NYLOC(2) 

= 10.5 inches 

WLa 

= NYLOC(3) 

= 91.3 inches 


The lateral and vertical accelerations with lever arm effects included and 1 g bias removed are, 
nys = ny + ((pq + f)xa - (p^+r^a + (qr-p)za)/g 


n^^ = -1 - Hz - ((pr-q)xa + (qr + p)ya - (p2+q2)za)/g 
The ISENS=1 output accelerations are bounded according to table 7- la. 

”ysi ""ys^5z 

nzsi -"zsi-io 


(AYSl) 

(AZSl) 


7.4 Rate Gyros (ISENS=1) 

The outputs of the rate gyros are simply the body frame roll, pitch, and yaw rates in units of 
degrees/sec with limits applied as inthcated in table 7. 1 a. The calculations are performed in macro 
RGYROl as follows. 


, 180 . ,300 „ 
PsI “ ( P) 

(PSD 

n 

, 180 , ,60 
9sl -( 9)1-60 

(QSi) 

n 

f 180 . [60 _ 

hi -( r) 1^) /sec 

7t 

(RSI) 
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7.5 Angle of Attack Measurement (ISENS=1) 

The angle-of-attack probe model is implemented in macro ALFBTl. The "probe" quantity is 
calculated by passing the angle of attack at the c.g. through a filter representing the probe 
dynamics. 


“p “ (0.073 s + r)“ 

The local probe angle of attack is adjusted for the induced upwash due to pitch rate, 

«p, = ttp - ( 180 /tc)*^ 

where Ip = 24 feet. From purely geometric consideration, one would expect Ip to be the distance 
from the center-of-gravity to the angle-of-attack probe where Ip has units of feet and is positive if 
the probe is forward of the c.g. A probe located at a fuselage station (FS) of 1 67 inches is 
consistent with Ip = 288 inches and a nominal c.g. location of FS^g = 455.0 inches, placing the 
probe between the nose (FS=60.50) and the leading edge of the cockpit (FS=188.0) (see figure 
3.1). However, one must wonder why the probe dynamics are not applied to the total local 
measurement that includes pitch rate effects instead of to the angle of attack at the c.g. At this 
point, the only answer is that this is the way it is done in subroutine FCSENS of the MDC 
simulation. 

An angle-of-attack probe bias and gain inverse are calculated. The bias is given by, 

{ 1 .9 degrees if gear up 
1.2 degrees if gear down. 

The gain inverse is a function of limited Mach number, 

GKot = FGKAOA( M l 


where the function Fqk/^q/^ is defined by figure 7.4. 


The final limited sensor output is then. 


«S1 


Pi Pbias 
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GKa 


-14 


(ALFSl) 
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0.68 r 



Figure 7.4. Inverse of Sensed AOA Gain 

7.6 Air Data Parameters (ISENS=1) 

The calculation of indicated static pressure and indicated compressible dynamic pressure requires 
eight parameters, (^Dx, where x={ A,B,...,H } . The QCDx arc 21 element arrays and the current 
value for each QCDx parameter is determined as a tabular function of Mach. XHM 1 is the 
similarly sized independent variable array of Mach numbers with a maximum value of 1.65 and a 
minimum of 0. A search is performed to determine the integer Ixhm ^HMl (Ixhm) ^ 

Mach < XHM1 (Ixhm+ 1)- T'*’® ^XHM'th elements of the QCDx arrays arc used in the following 
calculations. No interpolation is done on the QCDx variables. In the notation below, Qcdx “ 
QCDx(Ixhm)- 7^® QCDx parameters as functions of Mach are plotted in figure 7.5.a-h. 

The calculation of the sensed static pressure for ISENS=1, shown in equations 7.1 

through 7.3. 


SpSFXD ' QcdbM - QcdA ■ {Qcdb^^ + QcDcJ^S] 

PSs;! = Ps / SpSPXD 

PSTSl = Psj,, (7.3) 


The calculation of indicated compressible dynamic pressure is shown in equations 7.4 through 7.8. 

Pt =Qc Ps (7.4) 

Pt 


Pt. = 


1.0 - Qcdf^^ ■ Qcde ■ (Qcdh^ ^ Qcdg)®S] 


(7.5) 
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An uncorrected initial value for indicated compressible dynamic pressure is calculated as follows. 


Qcx = Pti + PSsl 


(7.6) 


The initial value for indicated compressible dynamic pressure is then adjusted for angle-of-attack 
effects per MDC memo 34 1 -624 1 dated "23Feb8 1 


Qcsi = QcxH -(0.04|a(” -301)2} 
QCISI 


(7.7) 

(7.8) 
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8. ACTUATORS 

The simulation described herein includes dynamic models for the following control surfaces, 

(1) stabilators (right and left), 

(2) rudders (right and left), 

(3) ailerons (right and left), 

(4) trailing-edge flaps (right and left), 

(5) leading-edge flaps (right and left), 

(6) speed brake, and 

(7) thrust vectoring vanes (2 per engine, 4 total). 

The stabilators, rudders, ailerons, trailing-edge, and leading-edge flaps are hydraulically powered 
and are referred to in both MDC documents and this report as the prim^ controls. The single 
speed brake is positioned between and forward of the vertical fins and is driven at a constant rate of 
24 degrees per second. The lower and upper limits of the speed brake deflection are 0° and 60°. 

The primary controls and the speedbrake are found on the basic F/A-18 airframe. 

The thrust vectoring vane actuation systems modeled in this simulation reflect a preliminary design 
of the thrust vectoring system installed on the F- 1 8/H ARV and do not agree with the final 
installation. The documentation of the final system will be undertaken in a separate repon. The 
thrust vectoring vanes are each modeled as first order transfer functions with a steady state gain of 
one, a time constant of (1/30) seconds, rate limits of 80 7s and position limits of +/- 30°. 

8.1 Linear Dynamics - Primary Controls 

The actuation systems for the primary controls in their normal mode of operation are discussed in 
detail in pages 6-1 through 6-40 of MDC A7813, vol II. The servo loop block diagrams for the 
primary controls are given in figures 6-2, 6-4, 6-6, 6-8 and 6-10 of MDC A7813, vol II. These 
block diagrams are repeated in a simplified form in this document as figures 8.1 through 8.5. The 
high order models represented in these figures describe main ram response with external load. 

The transfer functions for the five primary actuators (stabilator, rudder, aileron, leading-edge flaps, 
and trailing-edge flaps) are presented in Table 8.1 through 8.5. Each table contains three transfer 
function models of the corresponding actuator type; (a) a high order model; (b) an intennediate 
order model; and (c) a low older model. For each actuator type, transfer function (a) is derived 
from the high order block diagrams in figures 8. 1 through 8.5. Transfer functions (b) are 
discussed in pages 6-18 through 6-24 of MDC A7813, vol II. The (b) transfer functions are 
second order for all except the stabilators, which are fourth order. Transfer functions (c) consists 
of the lowest order models, the first order transfer functions that arc used in both the MDC 
simulation code and the DMS real-time simulation at LaRC. The (a) transfer functions can be 
realized by linearizing the nonlinear block diagrams of figures 8.1 to 8.5. Comparisons of gain 
and phase response of the (a), (b), and (c) transfer functions for the primary controls are shown in 
figures 8.6 through 8.10. The low order models agree pretty well in phase with the intermediate 
and high order models out to about 10 radians/second. 

An initial form of the ACSL simulation included elastic dynamics (see section 1.2), and in that 
simulation the (b) transfer functions were implemented. When the elastic dynamics were removed 
and time histories were compared with results generated by the DMS simulation, significant 
differences in aircraft responses due to pilot pitch and yaw doublets resulted. These differences 
were much reduced when the stabilator and rudder actuator models were changed to the first order 
models. The ACSL simulation now uses the first order actuator models. 
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STABILATOR SERVO LOOP BLOCK DIAGRAM 



COMMAND 



RUDDER SERVO LOOP BLOCK DIAGRAM 

Figure 8.2. Rudder Servoactuator Block Diagram 


COMMAND 



AILERON SERVO LOOP BLOCK DIAGRAM 


Figure 8.3. Aileron Servoactuator Block Diagram 
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ELECTRO- 



TRAILING EDGE FLAP SERVO LOOP BLOCK DIAGRAM 




3 


LEADING EDGE FLAP SERVO LOOP BLOCK DIAGRAM 




Table 8.1 . Linear Stabilator Actuator Models 


S(s) 

6c(s) 




+ 1 ) 


((^)2+2(.359)(3|3)+l) ((]f3)2+2{.637X-^)+l)((35p2+2(.MX^^ +^><1098 ^ 

((^) 2 . 2 (. 068 X^Ki) 


(a) 




Table 8.2. Linear Rudder Actuator Models - 


5(s) 

5c(s) 




(^*1)(^*1)((5|3)^-XWX5J3K1) 

1 

1 




(b) 


(c) 


(a) 


(b) 


(c) 


Table 8.3. Linear Aileron Actuator Models 


5(s) 

Sn(s) 


M03' 


•’”<356 *" 

>1)(4;,+')«75/*2(.789x4).1) 

(a) 

1 

(b) 

(:^)2+2(.59X:^)+1 


1 

(c) 
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Figure 8.6. Gain and phase response of the (a), (b) and (c) transfer functions for the stabilator 






Figure 8.7. Gain and phase response of the (a), (b) and (c) transfer functions for the rudder 
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Figure 8.8. Gain and phase response 


of the (a), (b) and (c) transfer functions for the aileron 
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Figure 8.10. Gain and phase response of the (a), (b) and (c) transfer functions for 

the leading-edge flaps 
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8.2 Nonlinearities > Primary Controls 

The signiflcant nonlinearides in the actuators include a minimum resolution in the digital command 
due to 12 bit digital-to-analog conversion, hydromechanical thresholds, rate limiting, and nonlinear 
gearing. The first three of these types of nonlinearities are given in Table 6.2. 1-1 of MDC A7813, 
vol n, and are repeated in Table 8.6 below. 


Table 8.6. Actuator Nonlinearities 


Actuator 

Minimum Digital 
Cbmmand (degrees) 

Threshold (fraction of 
full) 

No Load Rate Limit 
(degrees/second) 

StabilatOT 

0.0083 

0.10 

40 

Rudder 

0.0146 

0.05 

56 

Aileron 

0.022 

0.05 

100 

Trailing-Edge Flaps 

0.0127 

0.05 

18 

Leading-Edge Flaps 

0.0106 

0.15 

18 


The main ram slewing rate limit is proportional to (Pj - where Ps is the supply pressure 

(3000 PSI) and P] is the pressure on the resisting side of the piston due to external load. The 
maximum slewing rate occurs when P] is zero. When Pj equals Pg, the actuator stalls. 

In the ACSL simulation, of the possible nonlinear effects, only the no load rate limits are modeled. 
In both the MDC simulation and the DMS simulation, the no load rate limit used for the rudder 
actuators is 61 degrees/second instead of 56 degrees/second as shown in Table 8.6. Another set of 
specifications for actuator models can be found on page 12 of MDC report A8449 and is repeated 
below as Table 8.7. The no load rate limit for the rudder actuators is given as 61 degrees/second in 
Table 8.7. The 61 degree/second limit is used in the flHhas simulation. 


Table 8.7. 

MDC F/A-18 Flight Hardware Integration Simulator 

Actuator 

Rate Limit 

Position Limit 

Mechanization 

Stabilator 

+/- 407s 

-247+10.5° 

isl order lag, x = 1/30 s 

Rudder 

4-/-617s 

-3(),+30° 

151 order lag, x = 1/40 s 

Aileron 

+/- 100°/s 

-25°,+45° 

151 order lag, x = 1/48 s 

Trailing-Edge Flaps 

+/- 187s 

-8°,+45° 

15l order lag, x = 1/20 s 

Leading-Edge Flaps 

+/- 187s 

-3°,+34° 

151 order lag, x = 1/20 s 
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For further discussion of gearing nonlinearities and plots of the effect of hinge moment on 
maximum actuator rate, the reader is referred to section 6.2 of MDC report A781 3, vol II. 

8.3 Speedbrake 

The speedbrake system is described in section 8.7 of MDC A78I3, vol I, Rev A. The speedbrake 
actuator is operat^ by a selector valve which ports hydraulic pressure to the main ram for 
extension or retraction. The selector valve is electrically controlled by a three position switch, 

OFF, RETRACT, and EXTEND, on the right hand throttle. The EXTEND position is a 
momentary. The OFF command holds the speedbrake in its current position. When the 
speedbrake is fully retracted, both "OFF" and "RETRACT" switch positions command the selector 
valve to the retract position. 

The speedbrake rate is controlled by a flow regulator and restrictor. "The speedbrake extension 
time is limited to 2.5 seconds or approximately 20 degrees to 30 degrees per second (no load)" 
(MDC A7813, p8-14]. Maximum surface deflection is 60 degrees. In x\\cfI8bas simulation, the 

speedbrake no-load rate limit, SBNLRL, has been set to 60/2.5 = 24 °/s. No provision is made 
for reduced rate due to aerodynamic hinge moment. The maximum driving hinge moment available 
at 60 degrees open with the full 30(X) psi of hydraulic pressure available is 296,000 in-lbs. 

8.4 Thrust Vectoring Vane Actuators 

The four actuators in the preliminary vane actuation system are each modeled with first order lag 
transfer functions with a time constant of 0.05 seconds, a rate limit of 80 Vsec, and position limits 
of +/- 30*. Again, this vane configuration of two vanes per nozzle is a preliminary version. The 
actual system implemented on the HARV FI 8 has three vanes per nozzle. 
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9. FLIGHT CONTROL SYSTEM (FCS) 

T\)Rfl8bas simulation implements the Control Augmentation System (CAS) of the inner loop 
F/A-18 Flight Control System derined by the 8.3.3 production programmable read-only memory 
(PROM) set The the^ of operation and design considerations of the CAS is described in section 
9.1. The implementation issues associated with the fIShas simulation (as opposed to the actual 
aircraft) are discussed in section 9.2. 

9.1 FCS Theory of Operation 

The inner lo(^ flight control law theory of operation for the F/A-18 is described in detail in section 
16 of MDC report A78 1 3, Vol I, rev A and section 1 of MDC report A4107, vol I, rev J. The 
inner loop control laws are divided into three channels; longitudinal, lateral, and directional. A 
preliminary thrust vectoring option (LTHVEC=.T.) exists in the code wherein additional pitch and 
yaw commands are generated by the control laws for use by the mixer/predictor. The 
mixer/predictor resolves the pitch and yaw command into actual vane commands for proposed 
thrust vectoring h^dware. As indicated in section 1 .2 of this report, the thrust vectoring option in 
the fISbas simulation represents a preliminary design version that has been superseded and will not 
be described in this document. 

The inner loop CAS is a multi-rate digital controller. The actual flight computers update at a rate of 
160 Hz. The inner loop calculations are performed at alternate frames of the flight control 
computer, so that the fastest inner loop of any of the controllers is 80 Hz. The details of the 
control law discretization can be found in MDC report A4107, vol I, rev J. All the sampling rates 
are fast enough that the continuous description offered here will suffice for analysis. 

9. 1. 1 Longitudinal Auto Flap Up CAS 

Three modes of leading-edge and trailing-edge flap management are performed in the pitch inner 
loop control laws of the F/A-18; (1) auto flap up (AFU), (2), half flaps, and (3) full flaps. The 
AFU mode is for up and away flight. The half and full flap settings are used in take-off and 
landing, respectively. In the fISbas simulation, only the auto flap up mode is supported. A block 
diagram of the pitch inner-loop control augmentation system in the AFU mode is presented in 
figure 9.1 (a-d). Figure 9.1 is a simplified version of "FIG 1 /OFF VlO.l" on page 1-13 of MDC 
report A4 107, vol I.revJ. Paths in "FIG 1 /OFF VlO.l" having to do with outer loops, the 
direct electric link (DEL) reversion mode, and HALF and FULL flap management modes have 
been deleted. Also omitted from figure 9.1 are blocks defining stick dynamics, trim functions, and 
the g-limiter logic. The definition of the pitch gains (PKxx) in terms of pitch functions (FUNxx) are 
given in table 9.1. Table 9.2 defines the longitudinal filters used in figure 9.1. The pitch functions 
required to define quantities in both tables 9.1 and 9.2 are defined in figure 9.2 (a-o). The 
following quote is taken from pages 1-2 through 1-4 of MDC report A4107, Vol 1, REV. J. Italics 
are comments inserted by the authors. Grammatical hyphens have been added. 

Primary control of the F/A- 1 8 aircraft longitudinal motions is accomplished through the 
CAS. The Direct Electric Link (DEL) and mechanical system provide backup for the pitch 
CAS. Management of the leading and trailing edge flaps is provided through the use of a 
cockpit switch which calls for either AUTO, HALF, or FULL flaps. A block diagram of 
the longitudinal control system is presented in Figure 1 {figure 9.1 of this document). 
Lon^tudinal functions and logic diagrams arc given in Figure 1.1 through 1.44 (figure 9.2 
of this document). Longitudinal digital filter characteristics are summarized in Table 1-1 
(Table 9.2 of this document). 
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Auto Flap Up (CAS On) - The pitch CAS uses a pilot commanded longitudinal stick 
position input as a cwnmand to the CVVS. The forward path gain is air data scheduled 
(Function 32A) to yield a uniform initial pitch acceleration response for sha^ inputs. The 
CAS feedback parameters are a blend of air data scheduled pitch rate (Functions 40 and 
68), naroal acceleration, and angle of attack. Pitch rate and nonnal acceleration feedbacks 
give improved pitch dynamic characteristics and load factor control in the mid to high 
dynamic pressure portion of the flight envelope. Air combat maneuvering characteristics 
and increased stick-force-per-g cues in the low to mid dynamic pressure flight regime are 
provided by the air data scheduled pitch rate feedback. Angle-of-attack feedback provides 
additional increased stick force cues for low speed high angle-of-attack air combat 
maneuvering. Roll rate multiplied by yaw rate is fed to the longitudinal control system 
(Function 107) to reduce the effects of inertial coupling. 

The pilot commanded stick position, pitch rate, normal acceleration, and angle-of-attack 
feedback signals are iterated at the rates shown in Figure 1 (figure 9.1 of this document). 

The pitch CAS forward patch is comprised of an air data scheduled gain (Function 32A), a 
forward loop integrator, and the stabilator servo and power cylinder. The aircraft response 
(feedback) is compared to the maneuver command; the resulting difference drives the 
servo/power cylinder to reduce the maneuver error to zero through the forw^d loop 
integrator. Consequently, the aircraft is automatically kept in steady state trim (hands off) 
since any aircraft normal acceleration different from the commanded normal acceleration is 
reduced to zero by the action of the integrator. 

The notch structural filter located in the feedback paths attenuates any undesirable motions 
which are sensed by the motion sensors due to aeroelastic bending of the aircraft structure. 
The longitudinal forward loop gain (Functions 32A) and pitch rate feedback gains 
(Functions 40 and 68) are scheduled with air data. A supersonic (air data scheduled time 
constant. Function 22) lead-lag filter is employed in the pitch rate feedback path. The 
forward loop integrator gain is air data scheduled (Function 12) and angle-of-attack 
activated (AOA less than 33 degrees). This implementation is necessary to maintain 
adequate loop gain margins, phase margins, and handling qualities throughout the flight 
envelope. Scheduling pitch rate feedback gain with air d^ results in increased feedback 
gain with decreasing dynamic pressure or increasing altitude. This increases the stick force 
cues for mid to low speed flight conditions and provides a pitch attitude type of control law 
for tracking since pitch rate feedback is nulled via the forward loop integrator. 

The stall margin network introduces an aircraft nose-down command proportional to angle 
of attack (AOA) ateve 22 degrees AOA. The effect of the stall margin p^ion of the 
control laws is to increase the stick force cues for high angle-of-attack maneuvering. 

A speedbrake-to-stabilator interconnect (Function 100) is incorporated to alleviate g 
transients for speedbrake extension/retraction. The interconnect command is on/off in 
response to sp^brake solenoid valve command. The magnitude of the input command is 
air data scheduled. The interconnect command is lagged trough a first order low-pass 
filter with the time constant scheduled with air data and as a function of speedbrake 
extension or retraction. The speedbrake extend time constant is less than the speedbrake 
retract time constant. 

Flap Ctontrol - Auto Flap UP - With the flap switch in AUTO, the leading and trailing edge 
maneuvering flaps are scheduled with angle of attack and air data (Functions 24 and 27) to 
optimize performance, to improve the high angle of attack characteristic, and provide load 
alleviation at elevated g's. Rap position limits are computed as a function of air data inputs 
(Functions 25, 28 and 29) to prevent the application of excessive air loads on the flaps. 
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Table9J^Htd^^a^ 


Function 

AFUMode 

Spin Mode 

TK5 

Fli*F32A + F116 

0.0 

PK4 

1.5*F32A 

0.0 

PK15 

F32A 

0.7 

PK16 

3.5*F32A 

0.0 

PK17 

F32A 

0.0 

PK18 

F68*F32A + F40 

0.0 

PK19 

8.5*F32A + F12 

0.0 

PK20 

1.0 

0.0 

PK21 

F68*F32A + F12 

0.0 

PK22 

F32A*F12 

0.0 

PK23 

-22.0 

NA by virtue of PK17 

PK26 

0.0 

0.1 

PK29 

FI 07 

0.0 




Filter 

Transfer Function 

Where Used 

P2 

1 

ii 

Xn = .015(1+FUN22) 
Xd = .015 

Pitch Rate 

P5 

1 

Xd = 0.04 

Normal 

Acceleration 

Tq S + 1 

P8 

(S/C0)2 + 2^n(s/“) + 1 

Cn = 0.03 
Cd = 0.70 

CO = 60.0 (r/s) 

Notch 

(s/ 0))2 + 2Cd(s/(D) + 1 

Pll 

1 

Xd = 0.39 

Leading- 
Edge Flaps 

Xd s + 1 

P12 

1 

Xd = 0.79 

Trailing- 
Edge Raps 

Xd s + 1 

P15 

1 

Xd = 0.60 for extension 

0.49 <Xd <2.0 for 
retraction 

Speed brake 

Xd s + 1 

P16 

1 

Xd = 0.174 

NZAF 

Xd s + 1 
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Gain 



Qc(psi) 


(c) Pitch function 22 - Fader on the supersonic compensation 


Auto Flap Up (a, Rj) 


Scheduled 

position 

(deg) 



a (deg) 

(d) Pitch function 24 - Trailing-edge flap schedule 
Auto Flap Up (Qc) 

Figure 9.2 Continued. 
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(e) Pitch function 25 - Trailing-edge flap schedule 
Auto Flap Up (Qc) 



(0 Pitch function 27 - Leading-edge flap schedule 
Auto Flap Up (a, Rj) 

Figure 9.2 Continued. 


92 









(i) Pitch function 32A - Longitudinal forward loop gain schedule 

Auto Flap Up (Qc) 



(j) Pitch function 40 - Pitch rate feedback gain schedule 
Auto Flap Up (Pjj, Qc, Ri, STORES) 

Figure 9.2 Continued. 





(k) Pitch function 68 - Pitch rate feedback gain schedule 
Auto Flap Up (Qc) 



Qc (PSF) 


(1) Pitch function 107 - Longitudinal inertial gain schedule 
Auto Flap Up (Qq) 

Figure 9.2 Continued. 
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(m) Pitch function 1 15 - longitudinal nonlinear stick gradient 
Auto Flap Up (Q(.) 


P3^973PSF 
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(n) Pitch function 1 16 - Pitch nonlinear integrator gain schedule 
Auto Flap Up (Rj, Pj.,) 

Figure 9.2 Continued. 







9.1.2 Lateral Auto Flap Up CAS 

IlK roll fonctioos required ro *fiire quanntres in v<5"i!reV. L Italics 

1 Qtannii rrmtml i«! accomolished bv a contiol-by-wire system which comm^ds the ailerons, 
differential stabilators, differential trailing-edge flaps, and differential le^ng^ge fl^s. 
The ctxitrols consist of Control Augmentation System (CAS) and Distal and An^og D^t 
Electric Link (DEL) functions. A block diagram of the l^eral ctmtrol system is ^own m 
Hraire 2 (Figure 93 of this report). Lateral functions and logic dia^:am ^ presen^ in 
Figure 2.1 through 2.20 (Figure 9.4 of this report). Lateral digital filter charactensocs are 

summarized in Table 1-3 (Xoble 9.4 of this report). 

The mU CAS aucments the pilot's roll stick inputs and provides roll japing tlrough Ae 
aileron, differential stabilator, differential trmling^ge flaps, and 

flans The center stick position command signal is summed with Ac roll rate feed^k (p) 
to^rovide Ac closed loop control of Ac ailerons, Afferwoal stabilators, Afferenu^ 
tnulins-edge flaps, and Afferential leading-edge flaps. Becaure of Ae *"Jf 
damnine ofAe micraft at low altitude high speed flight coi^itions, the , 

eain^s Sliced (gain scheduled) in Ais area. A limited authority roll rate feedb^k is used 
Ki^e^ stebility in the area where the CAS roll rate feedback gam is reduced. 

The control stick position sensor signal, after it is shaped by a parabolic gradient, provides 
the desired roll response of the augmented aircraft. 

The roll CAS normal mode command path is comprised of an air data schooled g^ 
(Functions 4, 7, and 13). Function 4 is also the roll rate fc^back gain ^hedule. ^e gai 
varies wiA dynamic and static pressure to provide acceptable loop stability and roll 
rcsponre characteristics throughout the envelope. Notch filtera ^at^ ‘"f 
feedback path attenuate motion sensor inputs due to aeroclastic benAng of the aircraft 
smiturc.^A notch filter is located in Ae stick position signal to prevent struc^ mode 

couplmg through Ae stick position sensor. The frequencies ^ filter is 

determined during Ground Vibration Tesung. An adAdonal fou^ order roll-off filter s 
located in Ae feedback paA to provide adAtional attenuation of Ac structuiA mode inputs 
^ high frequency noise, the presence of which was detected during Iron Bird tesung. 

WiA wine stores the roll rate feedback gain is increased for roll rates in excess of 1 10 
r^KSind tmd the stick command gain (Function 13) is moAfied m order to limit 
dStairiS^uS toappmximatcly 1 50 degrees per «coni Wi^t wmg slews, the 
roll rale feedback gain is increased above 190 degrees per second to aid in limiting 
m»tauS^l mte*^ approximately 220 degrees per second. flVing stores ore no, an 

option in the flSbas simulation). 

The maximum roll rate capability with or without stores is limit^ to 15() de^s per 

wrmal load factor to reduce the vettical tail loads and reduce merual 

coupling. 
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The aileron command gain is scheduled with angle of attack and air data (Functions 35 and 
36). The angle-of-attack gain schedule eliminates excessive adverse sideslip during rolling 
maneuvers for high angles of attack at low dynamic pressure flight conditions. The air data 
gain schedule reduces the aileron command to zero at high dynamic pressure flight 
conditions where wing flexibility reverses the effective aileron aerodynamic moments. The 
air data schedule also avoids hinge moment limits. The aileron command is also limited by 
an angle of attack and air data schedule (Function 41 ) to optimize roll coordination. 

The roll CAS signal which drives the differential stabilator is multiplied by a gain and then 
is limited by a gain schedule (Functions 6 and 101) and rate limited at 35 degrees per 
second to reduce the coupling between the roll and pitch axis at high dynamic pressure 
where the stabilator approaches hinge moment and bending moment limits. Stabilator 
bending rnwrient alleviation is provided via a differential stabilator .surface command gain 
scheduled with load factor (Function 101). The differential stabilator command is also 
limited by the stabilator pitch command (pitch priority) and by an angle-of-attack schedule 
(Function 41). The ailertm and differential stabilator command signals also command 
rudder (for roll c(X)rdination) through the rolling surface-to-rudder interconnect (RSRI). 
The RSRI is gain .scheduled with air data and angle of attack (Functions 30 and 38). 

The differential trailing-edge flap command gain is scheduled with air data and angle of 
attack (Functions 31 and 34). The air data .schedule permits the maximum differential 
deflection available without encountering hinge moment limiting or excessive vertical tail 
loads. The angle-of-attack schedule reduces the available deflection above 5 degrees AOA 
and eliminates it above 10 degrees AOA. Differential flaps are not necessary to meet roll 
performance requirements above the.se angles of attack. 

The differential leading-edge flaps are gain .scheduled with air data and load factor 
(Function 93) to provide improved roll rate in the low to mid altitude transonic speed range. 
The gain is reduced for increasing load factor, increasing altitude, and decreasing Mach in 
the tran.sonic speed range. 

Differential ailerons are used for active suppression of a wing stores aeroelastic oscillation. 
Tlte forward .sensor assembly lateral accelerometer signal is filtered to extract the oscillation 
component, which then drives the ailerons. The active suppression system is switched on 
by an air data schedule if the heavy wing stores di.screte is set by the Mission Computer. 
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DEMODULATOR FILTER 



Figure 9.3. Lateral CAS - (a) Main feedback paths 
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Figure 9.3. Coniinued - (b) Aileron and differential stabilator commands 
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Figure 9.3. Continued - (c) Differential flaps 
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Figure 9.3. Concluded - (d) Servo command synthesis 
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Table 9.3. Lateral 

(roll) gains in 8.3.3 inner loop CAS 

Function 

AFU Mode 

Spin Mode 

Tin — 

TTTt 

0.0 

RK2 

3.22 

6.732 

RK5 

F7*(F4*F13) 

2.5 

RK6 

F4 

0.0 

RKIO 

1.33 

0.0 

RKll 

1.0 

300.0 

RK16 

F36 

F36 

RK17 

F6 

F6 

RK18 

F31 

0.0 

RK22 

0.0 

0.0 

YK18 

See directional axis, Table 9.5 


Table 9.4. Lateral (roll) filters in 8.3.3 inner loop CAS 


Filter 

Transfer Function 

Where Used 

R1 

1 

c = 0.80 

0) = 90.0 (r/s) 

TM order rolloff 
for roll rate 
sensor 

(s/ 0))2 + 2C(s/cd) + 1 

R2 

(s/0))2 + 2Cn(s/0)) + 1 

Cn = 0.05 
Cd = 0.50 
ci) = 31.0 (r/s) 

5 Hz notch for 
stick lateral 
position sensor 

(s/co)2 + 2Cd(s/w) + 1 

R3 

(s/ 0))2 + 2Cn(s/co) + 1 

Cn = 0.07 
Cd = 0.70 
0 ) = 44.0 (r/s) 

7 Hz notch for 
roll rate sensor 

(s/co)2 + 2Cd(s/0)) + 1 

R4 

(s/t 0)2 + 2^n(s/“) + 1 

Cn = 0.03 
Cd=0.70 
0) = 107 .0 (r/s) 

17 Hz notch for 
roll rate sensor 

(s/( 0)2 + 2Cd(s/«) + i 
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(c) Roll function 6 - Differential stabilator gain schedule, AFU mode, (Rj, P^) 


(125, 10) (325, 10) 


Ps>628 PSF 


P«<242 PSF 
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(d) Roil function 7 - Lateral command schedule gain, AFU mode, (Q^t Ps) 

Figure 9.4 Continued. 
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(e) Roll function 13 - Lateral command schedule gain, AFU mode, (Rj, Pj) 



(f) Roll function 31 - Differential trailing-edge flap gain schedule, AFU mode, (Rj, Pg) 

Figure 9.4 Continued. 
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AOA (deg) 


(g) Roll function 34 - Differential trailing-edge flap gain schedule, AFU mode, (a) 



AOA (deg) 


(h) Roll function 35 - Lateral forward l(X)p gain schedule, AFU mode, (a) 

Figure 9.4 Continued. 





■Gain 


(i) Roll function 36 - Aileron gain schedule, AFU mode, (Qc, Ps, Ri) 



(j) Roll function 39 - Rudder pedal to roll CAS interconnect gain schedule, AFU mode, (a) 

Figure 9.4 Continued. 
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(22.0, 0.25) Rj<0.27 


(16.7, 0.1)VV(19.0. 0.1) R>0.39 


AOA (deg) 


(k) Roll function 41 - Rolling surface limit schedule, AFU mode, (a, Rj) 


(0.434, 0.06) (0.56,0.06) (0.65.0.06) (0.74,0.06) 


P_-2200PSF / p„=1500PSF. 


p -1000PSF 
S: / 


Pg-2200PSF 


(0.294,0.0) (0.42,0.0) (0.51,0.0) (0.6, 0.0) 


(1) Roll function 93A - Differential leading-edge flap gain schedule component, 

AFU mode, (Rj, Pjjj)) 

Figure 9.4 Continued. 








(m) Roll function 93B - Differential leading-edge Bap gain schedule component, AFU mode, (Ps) 



(n) Roll function 93C - Differential leading-edge flap gain schedule component, 

AFU mode, (Nz^p) 

Figure 9.4 Continued. 
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(o) Roll function 101 - Differential stabilator load alleviation schedule, 
AFU mode, (Rj, Pj., 

Figure 9.4 Concluded. 
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9.1.3 Directional Auto Flap Up CAS 

Figure 9.5 (a-b) is a simplified version of the diagram "Figure 3 / OFF VlO.l" on page 1-99 of 
Nff)C report A4107, Vol I, REV. J. Paths having to do with outer loops and the direct electric link 
(DEL) reversion mode have been deleted. The definition of the directional gains (YKxx) in terms of 
yaw functions (FUNxx) are given in Table 9.5. Table 9.6 defines the directional filters used in 
figure 9.5. The yaw functions required to define quantities in both tables 9.1 and 9.2 are defined 
in figure 9.6 (a-1). The following quote is from pages 1-93, 1-94 of MDC report A4107, Vol 1, 
REV. J. ^mments inserted by the authors are in italics. Grammatical hyphens have been added. 

DiiecticHial control is accomplished by a control-by-wire (C!BW) system which commands 
the twin rudders. The yaw CAS is the primary mode of operation with a DEL (Direct 
Electric Link) mode from the rudder pedal force transducers to the rudder actuators 
available for CAS off operation. A rolling surface-to-rudder interconnect (RSRI) is 
incorporated to minimize the sideslip accompanying roll maneuvers. The block diagram is 
shown in Figure 3 (Figure 95 of this document). Directional functions and logic diagrams 
are given in Figures 3. 1 through 3.27 (Figure 9.6 of this document). Directional digital 
filters are summarized in Table 1-5 (Table 9.6 of this document). 

The normal mode of control is closed loop employing stability axis yaw rate {r cos(a) - p 
sin(a)) and lateral acceleration feedback control. 

The rudder pedal force transducer signal, after it is shaped by a dead band followed by a 
parabolic gradient (Function 14), is air data scheduled (Function 10) to prevent a command 
which would exceed the vertical tail load limits. Air data and angle-of-attack gain schedules 
(Function 17 and 1 14) eliminate aircraft departures for full pedal inputs. Full surface 
authority is available at high angles of attack and low Mach numbers through the rudder 
pedal-to-rolling surface and rolling surface-to-rudder interconnects described below. 

The CAS feedback loop contains a third order roll-off filter to attenuate high frequency 
noise and structural feedback modes. The computed stability axis yaw rate feedback signal 
is multiplied by a gain (Function 96) and then shaped by a washout network. The yaw rate 
(multiplied by cosine of angle of attack) feedback component is used to augment the Dutch 
roll damping. The remaining signal component, proportional to roll rate multiplied by the 
sine of angle of attack, is used to help provide sideslip reduction during high angle of attack 
maneuvering flight. Lateral acceleration feedback aids in reducing sideslip and provides 
turn coordination. The lateral acceleration and yaw rate feedback gains are air data 
scheduled (Functions 90 and 96) to increase directional stability at high altitude supersonic 
speeds and to optimize Dutch roll characteristics at low to mid altitude and speeds, 
respectively. Additional lateral acceleration air data and angle of attack gain schedules 
(Functions 1 12 and 1 13) improve directional stability above Mach 0.7 and 15 degrees angle 
of attack. Roll rate multipli^ by pitch rate is fed to the directional control system to reduce 
the effects of inertial couping (Function 108) 

The forward loop consists of a gain scheduled with air data (Function 45). The gain 
reduces rudder pedal sensitivity in the high dynamic pressure region and optimizes the 
closed loop directional damping the roll ctxrrdination characteristics. 

The rolling surface-to-rudder interconnect (RSRI) is used to control rudder deflections to 
minimize the sideslip accompanying roll maneuvers. The RSRI command is the total 
commanded differential aileron plus differential stabilator. The RSRI signal is passed 
through a lead-lag network and multiplied by variable gains, which are functions of angle 
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of attack and air data (Functions 30 and 38). The RSRl allows doordinated rolls to be 
flown with the feet on the floor. 

A rudder pedal-to-rolling surface interconnect is provided to reduce sideslip and angle of 
attack excursions due to rudder deflection at high angles of attack. The interconnect gain is 
scheduled with angles of attack and is schedule to zero at low angles of attack 
(Function 39). 
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Figure 9.5. Concluded - (b) Rudder servo command synthesis 


Table 9.5. Directional (yaw) gains in 8.3.3 inner loop CAS 


Function 


AFU Mode 


0.00072*F10 

0.234*F10 


Spin Mode 


0.00144 
0.466 


Table 9.6. Directional (yaw) filters in 8.3.3 inner loop CAS 


Filter 


IS + 1 


(s/ 0))2 + 2^(s/co) + 1 


s 


xs + 1 


Xf^ s + 1 


XqS + 1 


xs + 1 


X = 1/54.6 

Oil order rolloff in 
lateral accel. 
feedback path 

C = 0.60 
0) = 37.7 (r/s) 

2nd order rolloff in 
lateral accel. 
feedback path 

X= 1.0 

washout used in 
yaw rate gyro path 

Xn = 0.75 
Xd = 0.50 

RSRl lead-lag 

X = 1/20.45 

ISI order rolloff in 
lateral accel. 
feedback path 





























(c) Yaw function 30 - RSRI gain schedule, AFU mode, (Q^, Ps) 


Gain 

(VoltsA/olt) 



(d) Yaw function 38 - RSRI gain schedule, AFU mode, (a, Rj, Ps) 
Figure 9.6 Continued. 
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(e) Yaw function 42 - RSRI nonlinear gradient, AFU mode 
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(f) Yaw function 45 - Directional forward loop gain schedule, AFU mode, (Qq, Pg) 

Figure 9.6 Continued. 
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Gain 
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(g) Yaw function 90 - Lateral acceleration feedback gain schedule, AFU mode, (Rj, Pj) 


Gain 

(volts/volt) 



(h) Yaw function 96 - Yaw rate gain schedule, AFU mode, (Q^, Ps) 

Figure 9.6 Continued. 
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(i) Yaw function 108 - Directional inertial gain schedule, AFU mode, (Qc) 


Gain 

(volts/volt) 


(0.28,4.0) (0.76,4.0) 



(j) Yaw function 1 12 - Lateral acceleration feedback gain, AFU mode, (Rj) 

Figure 9.6 Continued. 








9.2. Flight Control System Implementation 

The flSbas simulation flight control system (FCS) subroutines were created by starting with 
existing code which modeled flight controls in the Differential Maneuvering Simulator real-time 
F18 simulation (dmsflS). The code structure and variable nomenclature used in the dmsflS 
simulation were left largely unchanged. The interface portions of the control subroutines were 
modified to be compatible with the Advanced Continuous Simulation Language (ACSL) structure. 
When calling FORTRAN subroutines from an ACSL simulation, arguments may be passed 
through a parameter list or through the global ACSL common block, ZZCOM. Using the global 
common block is risky because tdl the ACSL defined variables arc in ZZCOM and duplicate names 
can easily exist in the called subroutine. In the interest of execution speed, real-time subroutines 
use a multiple common block structure to pass arguments instead of parameter lists. The existing 
real-time ctanmon blocks were replaced with parameter lists containing the feedback and driver 
values required by the FCS subroutines. Sensor and servoactuator models, which appeared in the 
real-time FCS subroutines, were removed and coded as continuous models in the ACSL main 
program. Minor modifications were required for the CDC FORTRAN V to VAX FORTRAN 
conversion. For example, multiple assignment statements can occur on a single line of code in 
CDC FORTRAN V (A=B=C), but not in VAX FORTRAN. Logical variables pertaining to real- 
time systems (i.e. the real-time monitor, cockpit lights, strip chart recorders, etc.) were removed. 
Library functions and subroutines were modified as required for compatibility. Extraneous code 
which was present for an alternate, expanded, flight control system was removed. A general 
cleanup of unused variables and dead end code was performed. 

Five subroutines comprise the flight control loop. An executive routine, FI 8FCS, is called from 
the discrete block, FCS, in the main ACSL routine. The call is scheduled at a regular time interval 
specified by the ACSL variable, TSFCS. The default interval value is 1/40 (.025) second which 
represents a compromise between the various sampling rates of 80, 40, 20 and 10 samples per 
second used by the actual aircraft flight control computer. Sensor signals and pilot control 
positions are passed to the control routines through the subroutine call parameter list. Surface 
position commands are then computed and returned to the ACSL main routine through this 
parameter list. Figure 10.1 presents the overall flow of the executive routine and shows only the 
main driver and f^back viues and the resulting surface deflections. Specific options and special 
conditions were not shown in this diagram to reduce the clutter. A complete list of passed 
parameters for subroutine F18FCS is given in tables 9.7 and 9.8. 

The executive routine, F18FCS, computes logical variables indicating mode and filter coefficients. 
The mode logical variables are: spin mode, stores, the supermaneuver switch, automatic 
speedbrake retraction, and limiter override commands. The filter coefficients are for the dynamic 
pressure digital filters and the angle of attack signal digital filters. The remaining four control 
routines are called from FI 8FCS in the following order: 

(1) the longitudinal control model, FI 8LONG, 

(2) the lateral control model, FI SLAT, 

(3) the directional control model, FI 8DIR, 

(4) and an auxiliary routine, F18AUX. 

The basic inputs for the longitudinal axis control routine, F18LONG, are: 

(1) the pilot longitudinal stick command, 

(2) normal accelerometer feedback, 

(3) pitch rate gyro feedback, 

(4) angle-of-attack sensor feedback, 

(5) dynamic pressure sensor scheduling input, 

(6) and an aileron to stabilator interconnect command from the lateral control routine. 
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LONGrrUDINAL 



Figure 9.7. Implementation of r\\cf1Hhas flight control laws 
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The basic outputs for the longitudinal axis control routine are: 

( 1 ) stabilator sur^e commands, 

(2) the leading and trailing-edge flap commands, 

(3) and the symmetric thrust vectoring command. 

The basic inputs for the lateral axis control routine, F18LAT, are: 

( 1 ) the pilot lateral sdck command, 

(2) roll rate sensor feedback, 

(3) normal accelerometers feedback, 

(4) lateral accelerometer feedback, 

(5) angle of attack sensor feedback, 

(6) dynamic pressure sensor scheduling input, 

(7) and a rudder to rolling surface crossfeed command from the directional routine. 

The basic ou^uts for the lateral axis control routine are : 

(1) aileron surface commands, 

(2) the interconnect command to the stabilator and flaps, 

(3) and a rolling surface to rudder interconnect command. 

The basic inputs for the directional axis control routine, F18DIR, are: 

(1) the pilot pedal command, 

(2) laie^ accelerometer feedback, 

(3) roll rate gyro feedback, 

(4) yaw rate gyro feedback, 

(5) angle of attack sensor feedback, 

(6) and a rolling surface to rudder interconnect command from the lateral routine. 

The basic outputs for the directional axis control routine are : 

( 1 ) differential rudder command, 

(2) symmetric rudder command, 

(3) and a differential thrust vectoring command. 

The auxiliary routine blends and schedules the thrust vectoring commands and transforms 
symmetric and differential surface commands into left and right specific surface commands. The 
auxiliary routine originally contained servoactuator and hinge moment models. Hinge moments 
were not being computed at the time of this documentation and the servoactuatws existed in the 
continuous portion of the simulation. 


Table 9.7. Input Parameters for subroutine F18FCS 


flSbas 

Local 

Definition 

ALFD 

ALPDOT 

d/dt of angle trf attack, rad/sec 

ALPS 

ALPS 

angle-of-attack sensor output, degrees 

AYS 

NYSENS 

lateral accelerometer output, + right, g’s 

AZS 

NZSENS 

vertical accelerometer output, + up, g's 

A81 

A81 

left engine outlet area, in2 

A82 

A82 

right engine outlet area, in^ 

BETD 

BETADOT 

d/dt sideslip angle, rad/sec 

BEnXj 

BETADEG 

sideslip angle, degrees 

BLCG 

YBASIC 

buttock line coordinates of the c.g., inches 

CPL 

CPL 

left throttle command, degrees 

CPR 

CPR 

right throttle command, ctegrees 
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Tabic 9.7. Input Parameters for subroutine F18FCS (continued) 

* I k JU ' ^. 1 ' — — lA— — ^ 


hSbas 

Local 

bel‘inition 

CSALF 

COSALP 

cosine of angle of attack 

DAASTR 

DELATM 

calculated antisymmetric aileton trim deflection, degrees 

DATMBT 

DATMBUT 

aileron trim button setting, degrees 

DLGPCT 

DLGPCT 

landing gear state, |0,1 1, 0=gear up, l=gear down 

DRASTR 

DELRTM 

calculated antisymmetric rudder trim position, degrees 

DSSY 

STBPP 

symmetric stabilator deflecdon, degrees 

DSSYTR 

DELSTM 

calculated symmetric stabilator trim position, degrees 

DSTMBT 

DSTMBUT 

stabilator trim button setting, degrees 

FCRSET 

FCRESET 

lateral trim reset logical for spin, T/F 

FCSV83 

FCSV833 

set pitch forward loop gain logical, T/F 

FSCG 

XBASIC 

fuselage station coordinates of the c.g., inches 

F7BYPS 

F7BYPA 

schedule lateral gain logical, T/F 

G 

G 

acceleration due to gravity, ft/sec^ 

GLIMTR 

GLtMlTR 

g limit on longitudinal stick logical, T/F 

GLORID 

GLORIDE 

g limit override logical, T/F 

HSTORE 

HSTORE 

heavy stores logical, T/F 

ISENS 

IDEALS 

integer, (0) = ideal sensors, (1) = MDC sensor model 

KFS 

KFS 

stick spring force constant, n.d. 

LTHVEC 

TREVERS 

thrust vectoring engaged, T/F 

LTRFLG 

TRMMING 

logical, true if simulation in trim loop, T/F 

MACH 

MACH 

Mach number 

MNSPIN 

MANSPIN 

logical, manual spin mode, T/F 

NTVOPT 

NTVOPT 

integer, thrust vectoring pitch control blend 

ORIDE 

ORIDE 

logical, override angle-of-attack limiter, T/F 

PS 

PSENS 

output of roll rate sensor, Vsec 

PCA 

DELTAP 

lateral pilot stick input, inches 

PCR 

DELTRP 

rudder pedal force input, lbs 

PCS 

DELTSP 

longitudinal pilot stick input, inches 

PCSTR 

DLSPTM 

calculated trim position for longitudinal stick, in 

PD 

PEXDT 

roll acceleration, rad/sec^ 

PSAS 

PSAS 

logical, pitch stability augmentation system engaged, T/F 

PSTS 

PSTATIC 

Sensed static pressure, lbs 

QS 

QSENS 

output of pitch rate sen.sor, Vsec 

QBAR 

QBAR 

dynamic pressure, Ibs/in^ 

gets 

QCIS 

sensed compressible impact pressure, lbs 

QD 

QDOT 

pitch acceleration, rad/sec^ 

RS 

RSENS 

output of yaw rate sensor, Vsec 

RA2DG 

DEGRAD 

(180/71) 

RD 

RDOT 

yaw acceleration, rad/sec^ 

RSAS 

RSAS 

logical, roll stability augmentation system engaged, T/F 

SETVDC 

SETTVDC 

differential thrust vectoring set, degrees 

SETVSC 

SETTVSC 

symmetric thrust vectoring set, degrees 

SNALF 

SINALP 

sine of angle of attack 
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Tabic 9.7. Input Parameters for subroutine F18FCS (concluded) 

* ' ^ 


riSbas 

Locai 

PniaitiTn"' 

STORES 

STORES 

logical, stores switch engaged, T/F 

T 

T 

simulation time, sec 

TSFCS 

H 

update frame time for flight control subroutines, sec 

TT 

TT 

total thrust, lbs 

TVDBOT 

TVDBOUT 

logical, thrust vectoring deadband removed, T/F 

VT 

VTOTAL 

total airspeed, fi/sec 

WLCG 

ZBASIC 

waterline coordinate of c.g., inches 

WT 

WEIGHT 

weight of airplane, lbs 

YSAS 

YSAS 

logical, yaw stability augmentation system engaged 


Table 9.8. Output Parameters for subroutine F18FCS 

'■ B5BBBaaas=B 


i^iSbas 

Local 

Definition 

CAL 

DELTACL 

left aileron command, degrees 

CAR 

DELTACR 

right aileron command, degrees 

CFL 

DELTFCL 

left trailing-edge flap command, degrees 

CFR 

DELTFCR 

right trailing-edge flap command, degrees 

CNL 

DELTNCL 

left leading-edge flap command, degrees 

CNR 

DELTNCR 

right leading-edge flap command, degrees 

CRL 

DELTRCL 

left rudder command, degrees 

CRR 

DELTRCR 

right rudder command, degrees 

CSB 

DELTSBC 

speedbrake command, degrees 

CSL 

DELTSCL 

left stabilator command, degrees 

CSR 

DELTSCR 

right stabilator command, degrees 

CTVLL 

CTVLL 

left lower vane (thrust vectoring) command, degrees 

CTVLU 

CTVLL 

left upper vane (thrust vectoring) command, degrees 

CTVRL 

CTVLL 

right lower vane (thrust vectoring) command, degrees 

CTVLU 

CTVLL 

right upper vane (thrust vectoring) command, degrees 

DLTVCL 

DLTVCL 

left overall vane (thrust vectoring) command, degrees 

DLTVCR 

DLTVCR 

right overall vane (thrust vectoring) command, degrees 


The call to the FI 8 flight control system subroutine (F18FCS) appears twice in the main ACSL 
program (fISbas). A flowchart of the main program which illustrates the trimming and main 
program loops is shown in figures 2.2 through 2.4 in section 2 of this repiort. 

The call to subroutine F18FCS in the initial. ..end code section of the simulation (figure 2.3) is 
required by the trimming procedure. F18FCS is only called if the ”FCS engaged” logical (LFCS) 
is true. The call to F18FCS in initial. ..end ctxle block is required because the ACSL translator 
sorts the FCS discrete block and places it after the derivative section where it is can then be 
checked for an event crossing, (T=TFCS in this ca.se). Actuator models represent physical 
systems (the servoactuator motors that drive the aircraft surfaces) and appear in the continuous 
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(derivative) portion of the program code. If the call to F18FCS were not present in the initial 
section of the ACSL program, the following sequence would occur: 

(a) The trim roudne would drive a pilot input. 

(b) Actuator responses would be computed before the discrete flight control 
system call. 

(c) The surface control command from the previous loop would drive the 
actuator. 

(d) No surface deflection response to the trial pilot input would be computed. 

(e) A singularity error would be noted in the trim routine. 

(0 ‘ESTIMATED [A] MATRIX SINGULAR' would be returned at ACSL 
prompt. 

(g) The simulation would remain untrimmed. 

As a result of entering the ACSL runtime command 'START, the simulation process begins 
execution at the top of the initial section. The top portion of the initial section code is executed 
once at the beginning of each run (i.e., each time START is entered at the ACSL prompt) and is 
composed primarily of type declarations (LCXllCAL, INTEGER) and constant statements. 
Constant statements arc implemented by the ACSL translator as Fortran DATA statements and the 
valws arc entered once into the computer memory when the program is linked. A few direct 
assignment statements, notably 1TFLG=1, also occur in top portion of the initial section code and 
these values arc reset with each run. The continuation statement, 1100..CONT1NUE, marks the 
t(^ of the trim loop. Code appearing after this continuation statement is executed with each loop of 
the trim sequence until the pass on which the variable ITFLG > 2 is true. Successful trim is 
indicted by IRFLG=3. The subroutine INITD is called following the 1100 statement. This call 
was included to clear the event list and prevent problems with the crossing value setting of the 
afterburner switch discrete . Initial values for sensor data are computed prior to the initial flight 
cmitrol system call. The first call to the flight control system computes surface deflection 
commands based on the initial trim values of pilot control ccMnmands. These initial trim pilot 
cmtrol values are read on the first pass through the initial section by the subroutine LOAD from the 
trim data file 'F18TR'. The trial pilot control values are used to compute trim are: 

(a) longitudinal stick position, PCSTR (£ilot Control Stabilator TRim>. 

(b) lateral stick position, PCATR, (Eilot Control Aileron TRim). 

(c) and pedal force, PCRTR (£ilot Control Rudder TRim^ 

The logical variable LFCS must be set true for the flight control system to be called. A false value 
of LFCS will cause the flight control system to be bypassed. The LFCS false condition is present 
to enable the simulation to be run from direct surface deflection commands. Trim is achieved when 
the first NXTR elements of the XTRIM array (see table 2.1 in section 2.) are driven to values that 
cause first NXTR elements of the YTRIM array (table 2.2) to be zero to within some tolerance. 

The specific combination of XTRIM values used to achieve trim is dependent on which type of 
dim condition is required (a straight and level trim should not require pedal and lateral stick 
adjustment, unless some asymmetry was introduced to the aerodynamics). This combination is set 
up in the trim condition file ('F18TR'). The prt)cedurc for setting up a trim case is more fully 
described in the section 2.2 of this report. 
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10. CONCLUSION 

A mathematical model and associated computer program (fISbas) to simulate a twin-tailed high 
performance fighter airplane (McDonnell Douglas F/A-18) were described. The simulation was 
intended to support advanced control law research for application to high performance aircraft. 

The modeled flight envelope is extensive, allowing investigations in the high-angle-of-attack 
portion of the flight regime. While no comparisons with F/A-18 flight data are offered, the 
simulation is sufficiently detailed to support proof-of-concept studies. 

The simulation program was written in the Advanced Continuous Simulation Language (ACSL). 
Wherever possible, continuous system dynamics were modeled in the DERIVATIVE portion of the 
ACSL program structure. This policy allows important features of ACSL, such as the automatic 
generation of linear models and error checking on the numerical integration process, to be utilized. 
Furthermore, ACSL, by implementing continuous and discrete event dynamics via the DERIVATIVE 
and DISCRETE block conventions, allows a coherent implementation of multirate digital control 
laws that should not be lightly discarded. 

The simulaticHi math model includes the nonlinear, six degree-of-freedom rigid-body equations, an 
engine model, sensors, and first-order actuators with rate and position limiting. A simplified form 
of the F/A-18 digital control laws (version 8.3.3) is implemented, including only the Auto Flap Up 
(AFU) flight mode used for air combat maneuvering. Aerodynamic forces juid moments are 
calculated from a wind-tunnel-derived database using table lodc-ups with linear interpolation. The 
aerodynamic database has an angle-of-attack range of -10® to +90® and a sideslip range of -20° to 
+20®. The effects of elastic deformation were incorporated in a quasi-static-elastic manner. Elastic 
degrees of freedom are not actively simulated. In the engine model, the throttle-commanded 
steady-state thrust level and the dynamic respon.se characteristics of the engine are based on airflow 
rate as determined from a table look-up. Afterburner dynamics are switch^ in at a threshold based 
on the engine airflow and throttle position. 

The source code can be obtained from the authors. 


129 



REFERENCES 


Advanced Continuous Simulation Language (ACSL) Reference Manual, Edition 4.2, Mitchell and 
Gauthier Associates, Concord, MA, 1987. 

Arbuckle, P.D., Buttrill, C.S., and Zeiler, T.A.; "A New Simulation Model Building Process for 
Use in Dynamics Integration Research"; AIAA Paper No. 87-2498-CP; Proceedings of 
the 1987 Flight Simulation Technologies Conference-, Monterey, CA, August, 1987. 

Buttrill, C.S., Zeiler, T.A., and Arbuckle, P.D.; "Nonlinear Simulation of a Flexible Aircraft in 
h^euvering Flight": AIAA Paper No. 87-2501 -CP; Proceedings of the 1987 Flight 
Simulation Technologies Conference-, Monterey, CA, August, 1987. 

Etkin, Dynamics of Atmospheric Flight-, John Wiley & Sons, Inc.; New York, 1972 

Format Subcommittee; "The SCI Continuous System Simulation Language (CSSL)", SCI 

Simulation Software Committee, Simulation Councils, INC, Simulation, Vol. 9, No. 6 
December, 1967. 

Gainer, T.G. and Hoffman, S.; "Summary of Transformation Equations and Equations of Motion 
Used in Free-Flight and Wind-Tunnel Data Reduction and Analysis", NASA SP-3070 


MDC A4107, "F/A-18 Flight Control Electronic Set Control Laws", Volume I, Issue date: 19 
March 76, Revision Date: 1 July 88, Revision letter: J, McDonnell Aircraft Company 

MDC A7247, "F/A-18 Stability and Control Data Report, Vol I: Low Angle of Attack", Issue date 
31 August 1981, Revision date 15 November 1982, Revision Letter B, McDonnell Aircraft 
Company 


MDC A7247, "F/A-18 Stability and Control Data Report, Vol II: High Angle of Attack", Issue date 
31 August 1981, McDonnell Aircraft Company 

MDC A8575, "F/A-18 Basic Aerodynamic Data", Issue date 31 March 1984, McDonnell Aircraft 
Company 


MDC A7813 - Volume I, "F/A-18A Flight Control System Design Report, Volume 1, System 

Description and Theory of Operation", Revisit>n date: 28 September 1984, Revision letter 
A, McDonnell Aircraft Company 

MDC A7813 - Volume II, "F/A-18A Flight Control System Design Report, Volume II, Flight 
Control System Analysis - Inner Loops", Issue date: 15 June 1984, McDonnell Aircraft 
Company 

MDC A7813 - Volume 111, "F/A-18A Flight Control System Design Report, Volume III, Flight 
Control System Analysis - Automatic Flight Modes", McDonnell Aircraft Company 

MDC A8575, "F/A-18 Basic Aerodynamic Data", Issue date 31 March 1984, McDonnell Aircraft 
Company 

Roskam, J,; Airplane Flight Dynamics and Automatic Flight Control-, Roskam Aviation and 
Engineering Company; Route 4, Box 274, Ottawa, Kansas; 1979 


130 



Zeiler, T.A., and Buttrill, C.S,; "Dynamic Analysis of an Unrestrained, Rotating Structure 
Through Nonlinear Simulation"; AIAA Paper 88-2232-CP; Proceedings of the 
AIAAIASMEIASCEIAHS 29di Structures, Structural Dynamics and Materials Conference, 
Williamsburg, VA, 1988. 


131 



SYMBOLS 


Symbol 

Definition 

Variable Program 
Name Unit 

a.r.c. 

aerodynamic reference center, in 



b 

wing span reference length used for force and moment 
coefficients, 37.42 ft 

BWRF 

CSL 

Bl^ 

buttock line coordinate of accelerometers, + right, in. 
(default = 10.5) 

NYL(X:(2) 

CSL 

BLcg 

buttock line coordinate of c.g., + right, in. (default = 0.0) 

FSCG 

CSL 

BLengL 

buttock line coordinate of engine/airframe interface point for left 
engine, + right, - 18.9, in. 

BLENGL 

CSL 

BLengj^ 

buttock line cootdinate of engine/airframe interface point for 
right engine, + right, 1 8.9, in. 

BLENGR 

CSL 

BLrf 

buttock line of aerodynamic reference center, + right, 0.0 in. 

BLRF 

CSL 

c 

wing chord reference length used for force and moment 
coefficients, 1 1 .52 ft 

CWRF 

CSL 


basic drag coefficient, controls at zero, rigid airplane, f(a,M) 

CDO 

SFAERRF 


total drag coefficient, steady flow, along -x axis in the stability 
frame, ^ 

CDREFO 

SFAERRF 

^D5f 

drag coefficient per unit symmetric trailing-edge flap, deg"^ 




f(a,M) 

CDDF 

SFAERRF 

^D8n 

drag coefficient per unit symmetric leading-edge flap, deg'^ 




f(a,M) 

CDDN 

SFAERRF 

eg. 

center-of-gravity, in 



^lo 

basic roll coefficient, controls at zero, rigid airplane, f(P,a,M) 

CIO 

SFAERRF 


roll coefficient per non-dimensional roll rate, 




f(a,M) 

CIP 

SFAERRF 
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C|a8h 

ClAS„ 


L& 





roll coefficient per non-dimensional yaw rate 



f(a.M) 

total roll coefficient, steady flow, about x axis in the body 
frame, Ci|^^ ^ 


ClR 


CIREFO 


change in roll coefficient per degree of differential trailing-edge 
flaps 

change in roll coefficient per degree differential stabilator 


CIFDDF 

DCIDDH 


change in roll coefficient per degree of differential leading-edge 
flaps 

total lift coefficient, along -z axis in the stability frame 


CIFDDN 


basic lift coefficient, controls at zero, rigid airplane, f(a,M) 


CLO 


lift coefficient per non-dimensional pitch rate, 9 Cl / 



f(a,M) 


total lift coefficient, steady flow, = 


CLQ 

CLREFO 


total lift coefficient per non-dimensional rate-of-change in a. 


ac ' 


dCt/a Jv 


, f(a,M) 


CLADRF 


rigid lift coefficient per non-dimensional rate-of-change in a. 


dCL/a 


/• _ 

' ac 

2V 


'rigid airplane), f(a,M) 


lift coefficient per unit symmetric railing-edge flap, rigid 
airplane, deg‘^ f(a,M) 

lift coefficient per unit symmetric leading-edge flap, rigid 
airplane, deg’^, f(a,M) 


CLAD 


CLDF 


CLDN 


SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

USAERRF 

USAERRF 

SFAERRF 

SFAERRF 
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Cm 

total pitch coefficient, about y axis in the body or stability frame 

Cmo 

basic pitch coefficient, controls at zero, rigid airplane, f(a,M) 

Cmq 

pitch coefficient per non-dimensional pitch rate, rigid airplane, 
f(a,M) 

Cmsf 

total pitch coefficient, steady flow, = 

Cm& 

total pitch coefficient per non-dimensional rate-of-change in a. 




rigid pitch coefficient per non-dimensional rate-of-change in a. 


^Cjn / 3^^j(rigid airplane), f(a,M) 

CmSf 

pitch coefficient per unit symmetric trailing-edge flap, rigid 
airplane, deg'*, f(a,M) 


pitch coefficient per unit symmetric leading-edge flap, rigid 
airplane, deg'*, f(a,M) 


basic yaw coefficient, controls zero, rigid airplane, f(P,a,M) 

s 

yaw coefficient per non-dimensional yaw rate. 
f(a,M) 


yaw coefficient per non-dimensional yaw rate, 
f(a,M) 

^nsf 

total yaw coefficient, steady flow, about z axis in the body 
frame, 

^«ASf 

change in yaw coefficient per degree differential trailing-edge 
flaps, (= 0.0) 
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CMO SFAERRF 

CMQ SFAERRF 
CMREFO SFAERRF 

CMADRF USAERRF 

CMAD USAERRF 

CMDF SFAERRF 

CMDN SFAERRF 
CNO SFAERRF 

CNP SFAERRF 

CNR SFAERRF 
CNREFO SFAERRF 

CNFDDF SFAERRF 



^"ASh 

change in yaw coefficient per degree differential stabilator. 



f(«» 5 h» 

DCNDDH 

SFAERRF 

‘^»AS„ 

change in yaw coefficient per degree differential leading-edge 
flaps, (= 0.0) 

CNFDDN 

SFAERRF 

CYq 

basic side force coefficient, controls at zero, rigid airplane. 




f(P,a,M) 

CYO 

SFAERRF 

% 

side force coefficient per non-dimensional roll rate. 




dCy/af^l^) f(a,M) 

CYP 

SFAERRF 

Cvr 

side force coefficient per non-dimensional yaw rate. 





CYR 

SFAERRF 

^Ysf 

total side force coefficient, steady flow, along y axis in the 
body or stability frame, q 

CYREFO 

SFAERRF 

^YA5f 

change in side force coefficient per degree differential trailing- 
edge flaps, (=0.0) 

CYFDDF 

SFAERRF 

CyaSh 

change in side force coefficient per degree differential stabilator 

DCYDDH 

SFAERRF 

^Ya5„ 

change in side force coefficient per degree differential leading- 



LlUfl 

edge flaps, (= 0.0) 

CYFDDN 

SFAERRF 

F 

force, lbs, eqn. (4.1) 



FSa 

fuselage station coordinate of the accelerometers, + aft, in. 
(default = 306.75) 

NYLOC(l) 

CSL 

FScg 

fuselage station coordinate of c.g., + aft, in. (default = 455.0) 

FSCG 

CSL 

FSeng 

fuselage station coordinate of engine/airframe interface point for 
right or left engine, + aft, 687.5, in. 

FSENGR 

CSL 

FSrf 

fuselage station of aerodynamic reference center, + aft, 458.56 
in. 

FSRF 

CSL 
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total x-body component of thrust-induced forces, lbs 

FVRE(l) 

CSL 

’’ye 

total y-body component of thrust-induced forces, lbs 

FVRE(2) 

CSL 


total z-body component of thrust-induced forces, lbs 

FVRE(3) 

CSL 

g 

acceleration due to gravity, 32.174 ft/sec^ 

G 

CSL 

Hp 

altitude appropriate for current density, p, in the 1962 Standard 
Atmosphere tables, feet 

HRF 

CSL 

^XX 

roll moment of inertia, slug-ft2 

DCX 

CSL 

Ixz 

roll-yaw cross coupling moment of inertia, slug-ft2 

IXZ 

CSL 

lyy 

pitch moment of inertia, slug-ft2 

lYY 

CSL 

^zz 

yaw moment of inertia, slug-ft2 

IZZ 

CSL 

Ixz 

(1,3) entry of the 3x3 direction cosine matrix [L| that 
transforms a vector in the local (Earth) frame to the body frame. 



^yz 

(2,3) entry of the 3x3 direction cosine matrix [L] that 
transforms a vector in the local (Earth) frame to the body frame. 



^zz 

(3,3) entry of the 3x3 direction cosine matrix |LJ that 
transforms a vector in the local (Earth) frame to the body frame. 



m 

mass, slugs, eqn. (4-1) 

MASS 

CSL 

M 

Mach number 



P 

+x component of rotational vekx’ity in the btxiy frame, rad/sec 

P 

CSL 

Ps 

static pressure (Ibs/ft^) 



q 

+y component of rotational velocity in the body frame, rad/sec 

Q 

CSL 

Q= 

compressible pressure 1 .4 *Ps*M 2/2.0 (lb.s/ft2) 



r 

+z component of rotational velocity in the body frame, rad/sec 

R 

CSL 

Ri 

compressible/static pressure QJP^ (lbs/ft2) 




elasiic/rigidmuliiplieronC|p,(=I.O) FRCIP sfaerrf 
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FRCIDA SFAERRF 




elastic/rigid factor applied to aileron roll increment, = R 

clastic/rigid effectiveness factor applied to roll due to rudder 
increment, f(M,Hp) FRClDR 

elastic/rigid factOT applied to stabilator lift coefficient increment, 
f(MJip) FRCLDH 


FLT 






flight-derived effectiveness factor applied to aileron lift 


coefficient increment, f(M,Hp) FFCLDA 

elastic/rigid effectiveness factor applied to aileron pitch 

coefficient increment, f(M,Hp) FFCMDA 

elastic/rigid factor applied to stabilator pitch coefficient 

increment, f(M,Hp) FRCMDH 

elastic/rigid multiplier on Cp^, f(M,a) FRCNP 

elastic/rigid factor applied to aileron yaw coefficient increment, 
f(M,Hp) FRCNPA 



elasticMgid effectiveness factor applied to yaw due to rudder 
increment, f(M,Hp) FRCNDR 

elastic/rigid multiplier on Cy , f(M,H,.) FRCYP 

Ip K 



elastic/rigid factor applied to aileron side force coefficient 
increment, f(M,Hp) 

elastic/rigid effectiveness factor applied to side force due to 
rudder increment, f(M,Hp) 


racYDA 


FRCYDR 


SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 
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effectiveness factor applied to roll due to rudder accounting for 
stabilator interference, f(a,8H) 

RCIDRH 

SFAERRF 


effectiveness factor applied to yaw due to rudder accounting for 
stabilator interference, f(a,5n) 

RCNDRH 

SFAERRF 


effectiveness factor applied to side force due to rudder 
accounting for stabilator interference, f(a,5j|) 

RCYDRH 

SFAERRF 

s 

wing reference area used for force and moment coefficients, 
400 

SWRF 

CSL 

t 

time, seconds 

T 

CSL 

THl 

total thrust produced by left engine, lbs 

THL 

CSL 

THr 

total thrust produced by right engine, lbs 

THR 

CSL 

TH*;L 

x-body component of force produced by left engine at the 
engine/airframe interface point (687.5,18.9,1(K).0), lbs 

THXL 

CSL 

TH,;R 

x-body component of force produced by right engine at the 
engine/airframe interface point (687.5,18.9,100.0), lbs 

THXR 

CSL 

THy;L 

y-body component of force prtxluced by left engine at the 
engine/airframe interface point (687.5,18.9,100.0), lbs 

THYL 

CSL 

'THy:R 

y-body component of force produced by right engine at the 
engine/airframe interface point (687.5,1 8.9, 1(X).0), lbs 

THYR 

CSL 

TH,;L 

z-body component of force produced by left engine at the 
engine/airframe interface point (687.5,18.9,100.0), lbs 

THZL 

CSL 

TH,;R 

z-body component of force produced by right engine at the 
engine/airframe interface point (687.5,18.9,100.0), lbs 

THZR 

CSL 

u 

+x component of inertial translational velocity in the body 
frame, ft/sec 

U 

CSL 

V 

+y component of inertial translational vekx:ity in the body 
frame, ft/sec 

V 

CSL 

V 

total airspeed, ft/sec 

VT 

CSL 

w 

+z component of inertial translational vekKity in the body 
frame, ft/sec 

W 

CSL 
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WLa 

waterline coordinate of accelerometers, + up, in. 
(default = 91.3) 

NYLOC(3) 

CSL 

WLcg 

waterline coordinate of c.g., + up, in. (default = 102.8) 

WLCG 

CSL 

WLeng 

waterline coordinate of engine/airframe interface point for right 
or left engine, + up, 100.0, in. 

WLENG 

CSL 

WLrf 

waterline coordinate of aerodynamic reference center, + up, 
100.0 in. 

WLRF 

CSL 


distance from c.g. to the a.r.c. along X-body axis, + forward, 
= - (FSrf- FScg)/12. ft 

XARF 

CSL 

Xcif 

distance from c.g. to the engine/airframe interface point along 
x-body axis, + forward, = - (FSeng - FScg)/12, ft 

XERF 

CSL 

yarf 

distance from c.g. to the a.r.c. along Y-btxly axis, + right, = 
(BLrf- BLcg)/12, ft 

YARF 

CSL 

y^L 

distance from c.g. to the left engine/airframe interface point 
along y-body axis, + right, = (BLcngL “ BLcg)/12, ft 

YERFL 

CSL 

ycrfR 

distance from c.g. to the right engine/airframe interface point 
along y-body axis, + right, = (BLcng^ - BLcg)/12, ft 

YERFR 

CSL 

Zaif 

distance from c.g. to the a.r.c. along Z-body axis, + down, = 
-(WLrf-WLcg)/12,ft 

ZARF 

CSL 

Zcrf 

distance from c.g. to the engine/airframe interface point along z- 
body axis, + down, = -(WLcng - WL(.g)/12, ft 

YERF 

CSL 

a 

angle of attack, degrees 

ALPDEG 

ALDGRF 

SFAERRF 

CSL 

P 

sideslip angle, degrees 

BETDEG 

BEDGRF 

SFAERRF 

CSL 


increment in drag coefficient due to lift, f(CL^j,M) 

DCDCL 

SFAERRF 


increment in drag coefficient due to left stabilator position, rigid 
airplane, f(a,5H^,M) 

DCDDHL 

SFAERRF 
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AC|)g increment in drag coefficient due to right stabiiator position, 
rigid airplane, f(a,5Hj^,M) 

ACd^lg increment to drag coefficient for fully deployed landing gear, 
f(a) 


AC, 




AC 


i>5rt 


increment in drag coefficient due to antisymmetric rudder, 

(= 0 . 0 ) 

increment in drag coefficient due to symmetric (toe-in) rudder, 

f(a, 6f/y^) 


ACpgs^ increment to drag coefficient for full 60° deflection in speed 
brakes, f(a,M) 


ACj 


asm 


ACi 


AC, 


roll increment due to airplane asymmetries, f(a) 

increment to C|^ due to elastic deformation, f(M,Hp) 

F 

p increment to C)p due to elastic deformation, f(M,Hp) 


ACjg 


Al 


AC] 


*6 


Ar 


increment in roll coefficient due to left aileron position, rigid 
airplane, f(a,5^j^,M) 

increment in roll coefficient due to right aileron position, rigid 
airplane, f(a,5^j^,M) 


AC|g^ increment in roll coefficient for full 20 degree deployment of 

symmetric trailing-edge flaps, rigid airplane, used with sign of 

P. f(IPI,a,M) 

ACjg increment in roll coefficient for full 25 degree deployment of 

symmetric leading-edge flaps, rigid airplane, used with sign of 

P, f(|p|,a,M) 


DCDDHR 

DCDDLG 

DCDDR 

DCDDRT 

DCDSB 

DCIASM 

DCIRF 

DCIBF 


DC 1 DAL 


DCIDAR 


DCIDF 


DCIDN 


SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 
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increment in roll coefficient due to antisymmetric rudder, rigid 
airplane, used with sign of 5 r^^, f(a,|5R^^|,M) 

DCIDR 

SFAERRF 


increment to roll coefficient for full 60° deflection in speed 
brakes per degree of sideslip, deg'*, f(a,|P|,M) 

DCISB 

SFAERRF 


increment to basic lift coefficient due to elastic deformation, 
f(a,M,Hp) 

DCLOF 

SFAERRF 


increment to due to elastic deformation 

DCLQF 

SFAERRF 


elastic increment to Cl^ , f(a,H) 

DCLADF 

USAERRF 


increment to Cl5j: elastic deformation, deg'* 

,f(a,M,Hp) 

DCLDFF 

SFAERRF 

*‘^ 5 „ 

increment to Ci ~ due to elastic deformation, deg'* 
,f(a,M,Hp) 

DCLDNF 

SFAERRF 

W 

increment in lift coefficient due to left aileron, rigid airplane. 

DCLDAL 

SFAERRF 

^Ar 

increment in lift coefficient due to right aileron, rigid airplane, 
f(a,5AL.M) 

DCLDAR 

SFAERRF 

^^LX 

®Hl 

increment in lift coefficient due to left stabilator position, rigid 
airplane, f(a,5^l^,M) 

DCLDHL 

SFAERRF 

^^LX 

increment in lift coefficient due to right stabilator position, rigid 
airplane, f(a,5Hj^,M) 

DCLDHR 

SFAERRF 

'^‘^Slg 

increment to lift coefficient for fully deployed landing gear, 
f(a) 

DCLDLG 

SFAERRF 
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DCLDR SFAERRF 


AC, 


^Sr 


AC 


lSj^t 


AC, 




AC 


Ac: 


m© 

F 

mq 

F 


increment in lift coefficient due to antisymmetric rudder, (=0.0) 

increment in lift coefficient due to symmetric rudder (rudder 

toe-in), f(a,5Rj,5f/A) dcldrt 

increment to lift coefficient for full 60 degree deflection in speed 
brakes, f(a,M) dclsb 

increment to basic pitch coefficient due to elastic deformation, 
f(a,M,Hp) DCMOF 

DCMQF 


increment to C^^ due to elastic deformation 


elastic increment to Cj„^, f(a,H) 


DCMADF 




Sf 


AC 


m 


8n 




LEX 


AC, 


ms 


Al 


AC, 


ms 


AC, 


m5, 


Hl 


DCMDDF 


Ar 


increment to Cj^^g^ due to elastic deformation, deg'^ 

,f(a,M,Hp) 

increment to C^g due to elastic deformation, deg'^ 

,f(a,M,Hp) 

increment in pitch due to the leading-edge wing extensions, 
f(a,M) 

increment in pitch coefficient due to left aileron, rigid airplane, 

deg'^ f(a,8^^,M) dcmdal 

increment in pitch coefficient due to right aileron, rigid airplane. 


DCMDNF 


DCMLEX 


deg'*, f(a,5^|^,M) 


DCMDAR 


increment in pitch coefficient due to left stabilator position, rigid 
airplane, f(a, ,M) DCMDHL 


SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

USAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 
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increment in pitch coefficient due to right stabilator position, 
rigid airplane, f(a,5^ ,M) 

increment to pitch coefficient for fully deployed landing gear, 
f(a) 

increment in pitch coefficient due to antisymmetric rudder, 
f(a,l5Rl,M) 


increment in pitch coefficient due to symmetric rudder (rudder 
toe-in), f(a,5RT,5f/A) 


DCMDHR 


DCMDLG 


DCMDR 


DCMDRT 


AC, 


m5< 


SB 


AC 


AC 


nr 

F 

np 


AC 


"8 


Al 


increment to pitch coefficient for full 60 degree deflection in 
speed brakes, f(a,M) 

increment to due to elastic deformation, f(M,Hp) 
increment to C^p due to elastic deformation, f(M,Hp) 

increment in yaw coefficient due to left aileron position, rigid 
airplane, f(a,8Aj^,M) 


DCMSB 

DCNRF 

DCNBF 


DCNDAL 


AC 


"8 


Ar 


increment in yaw coefficient due to right aileron position, rigid 
airplane, f(a,8Aj^,M) 


DCNDAR 



increment in yaw coefficient for full 20 degree deployment of 

symmetric trailing-edge flaps, rigid airplane, used with sign of 
p. f(|p|,a,M) 


DCNDF 



increment in yaw coefficient for full 25 degree deployment of 

symmetric leading-edge flaps, rigid airplane, used with sign of 
p, f(|pl,a,M) DCNDN 


SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 

SFAERRF 
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increment in yaw coefficient due to antisymmetric rudder, rigid 
airplane, used with sign of f(«.|SRj^j,I.M) 

DCNDR : 

SFAERRF 

^^"8sb 

increment to yaw coefficient for full 60 degree deflection in 
speed brakes per degree of sideslip, deg'*, f(a,||i|,M) 

DCNSB 

SFAERRF 


increment to Cy^. due to elastic deformation, f(M,Hp) 

DCYRF 

SFAERRF 

ACyp 

increment to Cyp elastic deformation, f(M,Hp) 

DCYBF 

SFAERRF 

ACy5^^ 

increment in side force coefficient due to left aileron position. 



rigid airplane, f(a,5^^,M) 

DCYDAL 

SFAERRF 


increment in side force coefficient due to right aileron position. 



rigid airplane, f(a,8^j^,M) 

DCYDAR 

SFAERRF 

aCysj 

increment in side force coefficient for full 20 degree deployment 



of symmetric trailing-edge flaps, rigid airplane, used with sign 
of p, f(lPl,a,M) 

DCYDF 

SFAERRF 

ACyx 

increment in side force coefficient for full 25 degree deployment 




of symmetric leading-edge flaps, rigid airplane, used with sign 
of p, f(|p|,a,M) 

DCYDN 

SFAERRF 

-^Cys^ 

increment in side force coefficient due to antisymmetric rudder, 
rigid airplane, used with sign of f(tt,l5j^^^|,M) 

DCYDR 

SFAERRF 


increment to side force coefficient for full 60 degree deflection 
in speed brakes per degree of sideslip, deg'*, f(a,lPl,M) 

DCYSB 

SFAERRF 

e 

pitch angle, radians 

THE 

CSL 

p 

density of the air, slugs/ft^ 

RHO 

CSL 

<I> 

roll angle, radians 

PHI 

CSL 
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'F 

yaw angle, radians 

PSI 

CSL 

^cngL 

left engine cant angle in x-y plane, pos nose right, (-1.98*) 

AZEL 

CSL 

^cngR 

right engine cant angle in x-y plane, pos nose right, (1.98*) 

AZER 

CSL 


Subscripts 

aerodynamic, with respect to the air mass 
due to air loads 
due to asymmetries 
engine induced 
due to gravity 

quantity computed at the aerodynamic reference center (a.r.c.) 
component due to winds, turbulence, or gust 
component along the X3 (body frame) axis 
component along the ¥3 (body frame) axis 
component along the Z3 (body frame) axis 


a 

A 

asm 

E,cng 

G 

rf 

wg 

X 

y 

z 


Notation 

{ ) column vector 

[] matrix or row vector 

f]' transpose of [ 1 

xF? the variable x restricted to the range, b < x < a, typically used to describe a function 

argument that is limited. 

Acronyms 

ACSL Advanced Continuous Simulation Language 

AFU Auto-Flap Up mode of the F/A- 1 8 Flight Control System (8.3.3 PROM) 

AGCB Aircraft Guidance and Controls Branch at NASA LaRC 
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BL 

CAS 

CDC 

DMS 

dmsflS 

FIT 

FS 

fISbas 

flSharv 

HARV 

LaRC 

LEX 

MDC 

mdcfIS 

PROM 

PSF 

RSRI 

SPIN 

STORES 

WL 


Buttock line coordinate, + right 
Control Augmentation System 
Control Data Corporation 

DifTerential Maneuvering Simulator twin-dome facility at NASA LaRC 

Real-time simulation program of an F/A-18 developed at LaRC for use in the DMS 
facility. 

Functional Integration Technology 
Fuselage station coordinate, + aft 

Simulation model and computer program described in this d(x:ument 

Simulation model and computer program resulting from extensive modifications to 
the fIShas program and in support of the HARV research program 

High-Angle-of-Attack Research Vehicle 

Langley Research Center 

Leading Edge Extension 

McDonnel Douglas Corporation 

Batch simulation program of an F/A- 1 8 developed by MDC and obtained for use at 
LaRC around 1986. 

Programmable Read Only Memory 

pounds per square foot 

Roll surface to rudder interconnect 

Spin mode of AFU in F-18 8.3.3 PROM set 

Logical for wing stores on or off (Always off) 

Waterline coordinate, -i- up 
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APPENDIX A - SIMULATION VARIABLES 
Simulation Variable Naming Plan 

E)^elq)ment of the simulation model involved the efforts of several engineers and programmers 
arra has included receipt of data, documentation, and several Fortran subroutines from MDC and 
other rources. This effort has given the authors a keen sense of the desirability of a "standard" 
variable nomenclature for coding simulation models. Without a "standard" variable nomenclature, 
r^as generally impossible to follow the coded logic without having the documentation in hand. 
Tills section dest^bes our efforts at maintaining a consistent variable naming strategy. The 
creation of meaningful variable names was made challenging by the restriction in ACSL version 9 
to SIX character variable names. ACSL version 10 supports longer variable names. 

The authors developed a system for naming variables. This should not be confused with issuing a 
variable dictiOTaiy for project personnel to use. The philosophy behind the variable nomenclature 
system was that any given variable name could be constructed from a set of predefined mnemonics 
(base names; ^fix and suffix modifiers), given an appropriate set of construction rules. This 
philosophy IS based upon ideas enunciated by Mitchell and documented in Arbuckle (1987) A set 
of mnemonics developed by the authors is listed in tables A.l through A.3. The construction rules 
for creating a variable name are as follows: 

1) In keeping with the ANSI Fortran 77 standard and the requirements of ACSL, all 
variable names are six characters or less. 

2) Prefix and suffix modifiers may only be used in combination with a base name. 

3) Multiple prefixes may only be used if their nesting levels are different. The highest 
nesting level prefix will be leftmost, followed by the base name. 

4) Multiple suffixes may only be used if their nesting levels are different. The highest 
nesting level suffix will be rightmost. 

5) Base name truncation may be used to meet the requirement that all names be six 
characters or less. If truncation is required, remove the right character of the base 
name. 


As an exainple of using the nomenclature system and construction rules, consider a variable 
defining P in units of degrees at the aerodynamic reference point. The base name for B is BET, the 
suffixes are DG and RF, respectively. The suffix DG has a nesting level of two, and the suffix RF 
has a nesting level of three (see table A. 3). Constructing the variable name using rule number five 
results in BETDGRF, which is a seven-character name and violates rule number one. Using rule 
number six, remove ihe right character of BET to give us BEDGRF as the variable name which 
satisfies the construction rules. 

When simulation model components were derived from preexisting code, the authors were 
reluctant to totally discard the variable name structure already in place. Not all of the variable 
names in the simulation were constructed according to the rules of this section. 



Table A.l Base Variable Naming Plan 


Base Name 

Description 

ALF 

angle of attack, a, radians 

BET 

sideslip angle, P, radians 

C 

control surface command, degrees 

Cxx 

direction cosine, xx = {XX, XY, XZ, YX, ... ,ZY, ZZ} 

D 

control surface deflection, degrees 

DG2RD 

degrees to radians conversion factor (n/180), rad/* 

En 

quaternion, n = 0,1, 2, 3 

GAM 

longitudinal flight path angle, radians 

H 

-z component of position in Eanh frame (altitude above 
mean sea level), ft 

L 

logical or integer switch 

MACH 

Mach number 

P 

X body component of rotational velocity, rad/sec 

PHI 

Euler bank angle, \|/, radians 

PSI 

Euler yaw angle, <{>, radians 

Q 

y body component of rotational velocity, rad/sec 

QBAR 

dynamic pressure, q, lb/ft2 

R 

z body component of rotational velocity, rad/sec 

RD2DG 

radians to degrees (180/n:), */rad 

RHO 

air density, slugs/ft^ 

TH 

total thrust, lb 

THE 

Euler pitch angle, 0, radians 

U 

X body frame component of translational velocity, ft/sec 

V 

y body frame component of translational velocity, ft./sec 

VT 

total velocity magnitude (inertial), ft/sec 

W 

z txxly frame component of translational velocity, ft/sec 

X 

X component of position in earth frame, ft 

Y 

y component of position in earth frame, ft 


Base Name E)escriprion 

CS trigonometric cosine of x 

R ratio 

SN OTgonometric sine of x 



?u??H^ev3'TesCTTpaon^ 

A A '■ * ij 1 J 


UWM 

rai 

intermediate value of a state 

2 

2 

quantity has been squared 

A 

1 

applies to aileron 

A 

2 

aerodynamic 

B 

2 

bias 

CG 

2 

airplane center-of-gravity 

D 

1 

d/dt, time rate of change 

E 

2 

engine, thrust related 

F 

1 

applies to trailing-edge flaps 

1 

2 

inertial 

IC 

5 

initial condition, variable is in the state variable i.c. array 

L 

4 

lift 

LL 

2 

lower limit 

LM 

2 

limit 

LN 

5 

dummy variable used only to act as a possible input for 



ACSL ANALYZ functions 

MN 

2 

minimum 

MX 

2 

maximum 

P 

2 

probe 

PLA 

1 

power lever angle 

R 

2 

rudder 

R 

4 

right 

RA 

2 

variable has units of radians or rad/sec 

RF 

3 

reference quantity 

S 

1 

stabilator 

S 

2 

sensor 

SB 

1 

speed brake 

SG 

4 

signal 

SL 

2 

sea level 

SY 

3 

symmetric 

TR 

5 

value calculated for trim, may not be a state variable 

UL 

2 

upper limit 

WG 

2 

aerodynamic disturbance, wind or gust 

X 

2 

X body frame component 

Y 

2 

y body frame component 

Z 

2 

z body frame component 

ZR 

2 

"zero", reference to an i.c. or trim state 


Simulation Dictionary Option 

ACSL provides for maintaining a simulation dictionary. If the logical LVDEF in the simulation is 
set true at run time, a call to LIS i'L)(29) is made in the terminal block of the translator source code. 
Calling LISTT)(29) causes a list of variable names and definitions to be read from unit 29 and 
written to the "PRN" output unit along with the values of all the variables in the simulation [ACSL 
(1987)]. All variables and their values are listed on the "PRN" output unit. The variables not 
found on the dictionary (unit 29) are simply listed without definitions. The ACSL executive 
program is not aware of any variables internal to called Fwtran subroutines. 
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List of Simulation Variables 


A81 

A82 

ABL 

ABLLL 

ABLLR 

ABR 

ABSWL 

ABSWR 

ABULL 

ABULR 

ALDGRF 

ALF 

ALFB 

ALFD 

ALFDG 

ALFDP 

ALFLM 

ALFP 

ALFPIC 

ALFS 

ALFSO 

ALFSl 

ALFTR 

AX 

AY 

AYS 

AYSO 

AYSl 

AYW 

AZ 

AZEL 

AZER 

AZS 

AZSO 

AZSl 

AZW 

BEDGRF 

BET 

BETD 

BETDG 

BETTR 

BLCG 

BLENGL 

BLENGR 

BLRF 

BWRF 

CIRFSF 

CAASLN 

CAASP 

CAL 

CAR 


left engine nozzle exit area, in^ 

right engine nozzle exit area, in^ 

marker for left afterburner "state update" discrete block 

left engine afterburner thrust lower rate limit, Ibs/sec 

right engine afterburner thrust lower rate limit, Ibs/sec 

marker for right afterburner "state-update" discrete block 

.true, if left engine afterburner is activated 

.true, if right engine afterburner is acdvated 

left engine afterburner thrust upper rate limit, Ibs/sec 

right engine afterburner thrust upper rate limit, Ibs/sec 

angle of attack at the aerodynamic reference center (a.r.c.), degrees 

angle of attack at the c.g., i^ians 

air data probe local angle-of-attack bias 

d/dt of ALF, rad/sec 

angle of attack at the c.g., degrees 

d/dt of ALFP, degrce/sec 

limited angle of attack, used in Pitot tube sensor model, degrees 

air data probe angle of attack, degrees 

air data probe angle-of-attack initial condition, degrees 

sensed angle of attack, equal to ALFSO or ALFSl based ISENS = 0 or 1 

actual or true angle of attack, equal to ALFDG, degrees 

output of air data models, degrees 

angle of attack required for trim, radians 

acceleration in x-body direction, g's 

acceleration in y-body direction, g's 

sensed lateral acceleration, AY, equal to AYSO or AYSl based on ISENS = 0 or 1, 
g's 

true lateral acceleration, = AY, g's 

output of lateral accelerometer sensor mrxlel, g's 

wind axis lateral acceleration, ft/sec^ 

acceleration in z-body direction, g's 

azimuth angle of left engine as installed, 1.98 degrees 

azimuth angle of right engine as installed, 1.98 degrees 

sensed normal acceleration, equal to AZSO or AZSl based on ISENS, g's 

true normal acceleration, equal to AZ, g's 

output of normal accelerometer sen.sor model, g's 

wind axis normal acceleration, ft/sec^ 

sideslip angle at the aerodynamic reference center (a.r.c.), degrees 

sideslip angle at the c.g., radians 

d/dt of BET, rad/sec 

sideslip angle at the c.g., degrees 

value of BET required for trim, calculated in trim search, radians 

buttock line coordinate of the center of gravity, inches 

buttock line coordinate of the engine/airframe interface, left engine, -18.9 in. 

buttock line coordinate of the engine/airframe interface, right engine, 18.9 in. 

buttock line coordinate of the aerodynamic reference center, 0.00 inches 

reference wing span used for nondimensional aero coefficients, 37.42 ft 

total steady- flow roll coefficient about x-body at a.r.c., n.d. 

command^ antisymmetric aileron, linearization input only, degrees 

dynamic check perturbation signal to antisymmetric aileron command, degrees 

commanded left aileron, degrees 

commanded right aileron, degrees 
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CASYLN 

CASYP 

CATL 

CATR 

CDRFSF 

CFASLN 

CFASP 

CFL 

CFR 

CFSYLN 

CFSYP 

QNT 


CLRFSF 

CLTR 

CMRFSF 


CNASLN 

CNASP 

CNL 

CNR 

CNRFSF 


CNSYLN 

CNSYP 


C02VT 

CPASLN 

CPASP 

CPL 

CPLTR 

CPR 

CPRTR 

CPSYLN 

CPSYP 

CPSYSG 

CRASLN 

GRASP 

CRL 

CRR 

CRSYLN 

CRSYP 

CSALF 


commanded antisymmetric aileron, linearization input only, degrees 
dynamic check perturbation signal to symmetric aileron command, degrees 
left power lever angle command from auto throttle, degrees 
right power lever angle command from auto throttle, degrees 

• 

drag coefficient at aerodynamic reference center, steady flow ( a = 0 ), n.d. 
cmnmanded antisymmetric trailing-edge flap, linearization input only, degrees 
dynamic check perturbation signal to antisymmetric trailing-edge flap command, 
degrees 

command to left trailing-edge flap actuator, degrees 

command to right trailing-^ge flap actuator, degrees 

commanded symmetric trailing-edge flap, linearization input only, degrees 

dynamic check perturbation signal to symmetric trailing-edge flap command, 

degrees 

simulation communication interval, time interval at which data is written to the 
master time history file or available at the screen, seconds 

lift coefficient at aerodynamic reference center, steady flow ( a = 0 ), n.d. 

Target trim value for CLRFSF if "12" is in active part of lYSEL array. 

pitch (about y-body axis) coefficient at aerodynamic reference point, steady flow 

( a = 0 ), n.d. 

commanded antisymmetric leading-edge flap, linearization input only, degrees 
dynamic check perturbation signal to antisymmetric leading-edge flap command, 
degrees 

command to left leading-edge flap actuator, degrees 
command to right leading-edge flap actuator, degrees 

yaw (about z-bcxiy axis) coefficient at aerodynamic reference point, steady flow 
( a = 0 ), n.d. 

commanded symmetric leading-edge flap, linearization input only, degrees 
dynamic check perturbation signal to symmetric leading-edge flap command, 
degrees 

intermediate quantity, ^ , sec 

commanded antisymmetric power lever angle, linearization input only, degrees 
dynamic check perturbation signal to antisymmetric power lever angle command, 
degrees 

left power lever angle command from the pilot, degrees 

left power lever angle command from the pilot required to for trim, degrees 

right power lever angle command from the pilot, degrees 

right power lever angle command from the pilot required to for trim, degrees 

commanded symmetric power lever angle, linearization input only, degrees 

dynamic check perturbation signal to symmetric power lever angle command, 

degrees 

commanded symmetric power lever angle, signal array used in IDYNCK=3 option 
wherein an arbitrary input stream can be constructed 

commanded antisymmetric rudder (both t.e.left), linearization input only, degrees 

dynamic check perturbation signal to antisymmetric rudder, degrees 

command to left rudder actuator, degrees 

command to right ruder actuator, degrees 

commanded symmetric rudder, linearization input only, degrees 

dynamic check perturbation signal to symmetric rudder command, degrees 

cosine of angle of attack at the c.g., cos(a) 
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CSALRF 

CSASLN 

CSASP 

CSAZEL 

CSAZER 

CSB 

CSBET 

CSL 

CSPHI 

CSPH02 

CSPS02 

CSR 

CSSYLN 

CSSYP 

CSTHE 

CSTHIV 

CSTHLM 


CSTH02 

CSTVOA 

CTVLL 

CTVLU 

crvRL 

CTVRU 

CWRF 

CXRFSF 

CXX 

CXY 

CXZ 

CYRFSF 

CYX 

CYY 

CYZ 

CZRFSF 

CZX 

CZY 

CZZ 

DAOL 

DAOR 

DAAS 

DAASLN 

DAASTR 

DADL 

DADR 

DAL 

DALIC 

DALLN 

DALLNl 


cosine of angle of attack at the a.r.c., cosCa^f) 

commanded antisymmetric stabilator, linearization input only, degrees 

dynamic check perturbation signal to antisymmetric stabilator command, degrees 

cosine of azimuth angle of left engine as installed, cos(AZEL) 

cosine of azimuth angle of right engine as installed, cos(A2^R) 

command to speed brake actuator, degrees 

cosine of sideslip angle at the c.g. 

command to left stabilator actuator, degrees 

cosine of the roll angle PHI 

cos( PHITR/2) 

cos(PSITR/2) 

command to right stabilator actuator, degrees 

commanded symmetric stabilator, linearization input only, degrees 

dynamic check perturbation signal to symmetric stabilator command, degrees 

cosine of the pitch angle, cos(0) 

= 1./ CSTHLM 

= cos(0) if lcos(0)l is larger than EPS, otherwise equal to EPS with the sign of 

cos(0). EPS = 10-*0 
cos(THETR/2) 

cosine of the thrust vectoring vane (mentation angle 
lower left thrust vector vane deflection command, degrees 
upper left thrust vector vane deflection command, degrees 
lower right thrust vector vane deflection command, degrees 
upper right thrust vector vane deflection command, degrees 

aerodynamic wing reference chord, c, 1 1 .52 ft 

aerodynamic axial force (x-body) coefficient at a.r.p., steady flow ( a = 0 ), n.d. 
direction cosine, x Ixxly component of unit x Earth frame vector 
direction cosine, x b<xiy component of unit y Earth frame vector 
direction cosine, x txxly component of unit z Earth frame vector 

aercxiynamic side force (y-lxxly) coefficient at a.r.p., steady flow ( (x = 0 ), n.d. 
direction cosine, y b(xiy component of unit x Earth frame vector 
direction cosine, y body component of unit y Earth frame vector 
direction cosine, y bexly component of unit z Earth frame vector 

• 

aercxiynamic normal force (z-bexiy) coefficient at a.r.p., steady flow ( a = 0 ), n.d. 

direction cosine, z b(xiy component of unit x Earth frame vector 

direction cosine, z body component of unit y Eanh frame vector 

direction cosine, z btxly component of unit z Earth frame vector 

position output of left aileron actuator, degrees 

position output of right aileron actuator, degrees 

antisymmetric aileron deflection, degrees 

component of total antisymmetric aileron deflection, linearization input only, 
degrees 

antisymmetric aileron deflection resulting from trim calculations, degrees 
d/dt of left aileron, degrees/sec 
d/dt of right aileron, degrees/sec 

total left aileron position, includes linearization input, degrees 
left aileron deflection initial condition, degrees 

component of total left aileron deflection, linearization input only , degrees 
= DALLN (if .not.LLNCSY) or (DASYLN - DAASLN) (if LLNCS Y). 
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DAR 

DARIC 

DARLN 

DARLNl 

DASY 

DASYLN 

DASYTR 

DATMBT 

DFOL 

DFOR 

DFAS 

DFASLN 

DFASTR 

DFDL 

DFDR 

DFL 

DFLIC 

DFLLN 

DFLLNl 

DFR 

DFRIC 

DFRLN 

DFRLNl 

DFSY 

DFSYLN 

DFSYTR 

DG2RA 

DLG 

DLGPCT 

DLTVCL 

DLTVCR 

DNOL 

DNOR 

DNAS 

DNASLN 

DNASTR 

DNDL 

DNDR 

DNL 

DNLIC 

DNLLN 

DNLLNl 

DNR 

DNRIC 

DNRLN 

DNRLNl 


total right aileron position, includes linearization input, degrees 
right aileron deflection initial condition, degrees 

component of total right aileron deflection, linearization input only, degrees 
= DARLN (if .not.LLNCSY) or (DASYLN + DAASLN) (if LLNCSY). 
symmetric aileron deflection, degrees 

component of total symmetric aileron deflection, linearization input only, degrees 

symmetric aileron deflection resulting from trim calculations, degrees 

aileron trim button, pilot input, (l.,0.,-l) 

position output of left trailing-edge flap actuator, degrees 

position output of right trailing-edge flap actuator, degrees 

antisymmetric trailing-edge flap deflation, degrees 

component of total antisymmetric trailing-edge flap deflection, linearization input 
only, degrees 

antisymmetric trailing-edge flap deflection resulting from trim calculations, degrees 
d/dt of left trailing-edge flap, degrees/sec 
d/dt of right trailing-^ge flap, degrees/sec 

total left trailing-edge flap position, includes linearization input, degrees 
left trailing-edge flap deflection initial condition, degrees 

component of total left trailing-edge flap deflection, linearization input only, degrees 

= DFLLN (if .not.LLNCSY) or (DFSYLN - DFASLN) (if LLNCSY). 

total right trailing-edge flap position, includes linearization input, degrees 

right trailing-edge flap deflection initial condition, deg^s 

component of total right trailing-edge flap deflection, linearization input only. 


degrees 

= DFRLN (if .not.LLNCSY) or (DFSYLN + DFASLN) (if LLNCSY). 
symmetric trailing-edge flap deflection, degrees 

component of totd symmetric trailing-edge flap deflection, linearization input only, 
degrees 

symmetric trailing-edge flap deflection resulting fran trim calculations, degrees 
degrees to radians conversion, (7t/180) 

landing gear deployment; assumes values in the interval [0.0, 1.0] where 0.0 
represents full up, and 1.0 represents full down. Default is full up. 

% DLG, (AOA probe) 

left ovei^l thrust vectoring vane deflection command, degrees 
right overall thrust vectoring vane deflection command, degrees 
position output of left leading-edge flap actuator, degrees 
position output of right leading-edge flap actuator, degrees 
antisymmetric leading-edge flap deflection, degrees 

component of total antisymmetric leading-edge flap deflection, linearization input 
only, degrees 

antisymmetric leading-edge flap deflection resulting from trim calculations, degrees 
d/dt of left leading-edge flap, degrees/sec 
d/dt of right leading-edge flap, degrees/sec 

total left leading-edge flap position, includes linearization input, degrees 

left leading-edge flap deflection initial condition, degrees 

component of total left leading-edge flap deflection, linearization input only, 

dcsrccs 

= DNLLN (if .not.LLNCS Y) or (DNSYLN - DNASLN) (if LLNCSY). 
total right leading-edge flap position, includes linearization input, degrees 
right leading-edge flap deflection initial condition, degrees 
component of total right leading-edge flap deflection, linearization input only. 


degrees 

= DNRLN (if .not.LLNCSY) or (DNSYLN + DNASLN) (if LLNCSY). 
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DNSY 

DNSYLN 

DNSYTR 

DPASLN 

DPASTR 

DPLDL 

DPLDR 

DPLL 

DPLLIC 

DPLLN 

DPLLNl 

DPLR 

DPLRIC 

DPRLN 

DPRLNl 

DPSYLN 

DPSYTR 

DRQL 

DROR 

DRAS 

DRASLN 

DRASTR 

DRDL 

DRDR 

DRL 

DRLIC 

DRLLN 

DRLLNl 

DRR 

DRRIC 

DRRLN 

DRRLNl 

DRSY 

DRSYLN 

DRSYTR 

DSOL 

DSOR 

DSAS 

DSASLN 

DSASTR 

DSB 

DSBD 

DSBIC 

DSBLN 

DSDL 

DSDR 

DSL 

DSUC 

DSLLN 


symmetric leading-edge flap deflection, degrees 

component of total symmetric leading-edge flap deflection, linearization input only, 
degiMs 

symmetric leading-edge flap deflection resulting frwn trim calculations, degrees 
component of tot^l antisymmetric power lever angle input, linearization only, 
degrees 

antisymmetric power lever angle deflection resulting from trim calculations, degrees 

d/dt of left throttle boost output signal, degrees/sec 

d/dt of right throttle boost output signal, degrees/sec 

left throttie boost output signal, degrees 

left throttle boost output signal initial condition, degrees 

added to DPPL, linearization input only, degrees 

= DPLLN (if .not.LLNCS Y) or (DPSYLN - DPASLN) (if LLNCS Y). 

right throttle boost output signal, degrees 

right throttle boost output signal initial condition, degrees 

added to DPPR, linearization input only, degrees 

= DPRLN (if ,not.LLNCSY) or (DPSYLN + DPASLN) (if LLNCS Y). 

component of total symmetric throttle boost output signal, linearization input only, 

degrees 

symmetric power lever angle deflection resulting from trim calculations, degrees 

position output of left rudder actuator, degrees 

position output of right rudder actuator, degrees 

antisymmetric rudder deflection, both t.e. left, degrees 

component of total antisymmetric rudder deflection, linearization input only, 

degrees 

antisymmetric rudder deflection resulting from trim calculations, degrees 
d/dt of left rudder, degrees/sec 
d/dt of right rudder, degrees/sec 

total left rudder position (pos t.e.left), includes linearization input, degrees 
left rudder deflection initial condition, degrees 

component of total left rudder deflection, linearization input only, degrees 
= DRLLN (if .not.LLNCS Y) or (DRSYLN - DRASLN) (if LLNCS Y). 
total right rudder position (pos t.e.left), includes linearization input, degrees 
right rudder deflection initial condition, degrees 

component of total right rudder deflection, linearization input only, degrees 
= DRRLN (if .not.LLNCSY) or (DRSYLN + DRASLN) (if LLNCSY). 
symmetric rudder deflection, both t.e. in, degrees 

component of touil symmetric rudder deflectitm, linearization input only, degrees 
synrimetric rudder deflection resulting from trim calculations, degrees 
position output of left stabilator actuator, degrees 
position output of right stabilator actuator, degrees 
antisymmetric stabilator deflection, degrees 

component of total antisymmetric stabilator deflection, linearization input only, 
degrees 

antisymmetric stabilator deflection resulting from trim calculations, degrees 

speed brake deflection, degrees 

(Vdt DSB, degrees/sec 

initial condition for DSB, degrees 

component of total speed brake deflection, linearization input only, degrees 
d/dt of left stabilator, degrees/sec 
d/dt of right stabilator, degrees/sec 

total left stabilator position, includes linearization input, degrees 
left stabilator deflection initial condition, degrees 

component of total left stabilator deflection, linearization input only, degrees 
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DSLLNl = DSLLN (if .not.LLNCSY) or (DSSYLN - DSASLN) (if LLNCSY). 

DSR total right stabilator position, includes linearization input, degrees 

DSRIC right stabilator deflection initial condition, deg^s 

DSRLN component of total right stabilator deflection, linearization input only, degrees 
DSRLNl = DSRLN (if .not.LLNCSY) or (DSSYLN + DSASLN) (if LLNCSY). 

DSSY symmetric stabilator deflection, degrees 

DSSYLN component of total symmetric stabilator deflection, linearization input only, degrees 
DSSYTR symmetric stabilator deflection resulting from trim calculations, degrees 
DSTMBT pilot stick pitch trim button 

DTVfflX lower left thrust vectoring vane deflection, actuator output, degrees 

DTVW^U upper left thrust vectoring vane deflection, actuator output, degrees 

DTVffllL lower right thrust vectoring vane deflection, actuator output, degrees 

DTVORU upper right thrust vectoring vane deflection, actuator output, degrees 

DTVDLL d/dt of lower left thrust vectoring vane deflection, degrees/sec 

DTVDLU d/dt of upp>er left thrust vectoring vane deflection, degrees/sec 

DTVDRL d/dt of lower right thrust vectoring vane deflection, degrees/sec 

DTVDRU d/dt of upper right thrust vectoring vane deflection, degrees/sec 

DTVL total left engine thrust vectoring vane deflection, degrees 

DTVLL left lower thrust vectoring vane deflection, degrees 

DTVLLN added to total left engine thrust vectoring vane deflection, linearization input only, 
degrees 

DTVLU left upper thrust vectoring vane deflection, degrees 
DTVR total left engine thrust vectoring vane deflection, degrees 

DTVRL right lower thrust vectoring vane deflection, degrees 

DTVRLN a^ed to total right engine thrust vectoring vane deflection, linearization input only, 
degrees 

DTVRU right upper thrust vectoring vane deflection, degrees 

DX 1 interme^ate variable in pressure sensor function calculations 

DX2 intermediate variable in pressure sensor function calculations 

EO quaternion, n.d. 

EOO EO prior to normalization 

EOOO Direct integral of EOD = d/dt EO 

E02 E0**2, used in direction cosine calculations 

EOD d/dt of EO, = ,5*(-El*P-E2*Q-E3*R), rad/sec 

EOIC Initial condition for EO, derived from trim Euler angles 

El quaternion, n.d. 

ElO El prior to normalization 

ElOO Direct integral of EID = d/dt El 

E12 El ** 2 , used in direction cosine calculations 

EID d/dt of El , = .5*( E0*P-E3*Q+E2*R), rad/sec 

ElIC Initial condition for El, derived from trim Euler angles 

E2 quaternion, n.d. 

E20 E2 prior to normalization 

E200 Direct integral of E2D = d/dt E2 

E22 E2**2, used in direction cosine calculations 

E2D d/dt of E2, = .5*( E3*P+E0*Q-E1 *R), rad/sec 

E2IC Initial condition for E2, derived from trim Euler angles 

E3 quaternion, n.d. 

E30 E3 prior to normalization 

E3(X) Direct integrd of E3D = d/dt E3 

E32 E3**2, us^ in direction cosine calculations 

E3D ^dt of E3, = .5*(-E2*P+E 1 *Q+EO*R), rad/sec 

E3IC initial conation for E3, derived from trim Euler angles 
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EPS 

EVENT 

F7BYPS 

FCRSET 

PCS 

FCSV83 

FGKAOA 

FGRESL 

FGRESR 

FHIL 

FHIR 

FHMl 

FRML 

FRMR 

FSCG 

FSENG 

FSRF 

FTD 

FTDAT 

FVR 

FVRA 

FVRE 

FVRGR 

FVRI 

FVRWG 

G 

GAM 

GAMDG 

GAMTR 

GAMZR 

GDTHR 


GINV 

GKIL 

GKIR 

GK3L 

GK3R 

GKAOA 

GLIMTR 

GLOAD 

GLORID 

H 

HD 

HGC 


used to prevent divide by zero, lO"*® 

time at which dynamic check perturbation is applied, sec 

FCS logical for lateral gain schedule bypass, default = .F„ (DO NOT 

(31ANGE) 

FCS logical to reset lateral trim integrator, default = .F., (DO NOT CHANGE) 
marker for flight control system discrete block 

FCS logical for pitch forward loop gain (8.3.3 PROM set), default = .T., (DO 
NOT CHANGE) 

= Mach number, but limited so that: 0.9 < FGKAOA < 1.4. Used as independent 
variable for sensor angle-of-attack gain value look-up. 

gross force produced by left engine prior to inclusion of thrust vectoring effects, 
lbs. 

gross force produced by right engine prior to inclusion of thrust vectoring effects, 
lbs. 

left engine thrust vectoring vane deflection angle, degrees 

right engine thrust vectoring vane deflection angle, degrees 

= Mach number, but limited so that: 0.0 < FHMl < 1.65. Used as independent 

variable in static pressure sensor I(K)k-up. 

left engine ram drag, lbs 

right engine ram drag, lbs 

fuselage station coordinate of the center of gravity, inches 

fuselage station coordinate of the engine/aiiframe interface, either engine, 687.5 in. 

fuselage station coordinate of the aerodynamic reference center, inches 

event marker for flight test data control input discrete block 

array used to hold 1 record of flight data read off a file, 50 elements. 

total generalized force vector for rigid-body dynamics 

aerodynamic component of FVR 

thrust induced component of FVR 

gravity induced component of FVR 

inerti^ component of FVR 

aerodynamic disturbance (wind and gust) compxxient of FVR 
acceleration due to gravity, 31.174 ft/sec^ 
longitudinal flight path angle, radians 
longitudinal flight path angle, degrees 

target longitudinal flight path angle used in trim search, radians 
GAM - GAMTR, driven to zero by trim search 

scaling parameter placed on a random number applied to next trim search solution 
candidate to prevent one component of the solution vector from reaching a solution 
earlier and remaining fixed. An unchanging element of solution vector will lead to a 
singularity in the sensitivity matrix estimation. 

1/g, used to scale various variables 

gain on left engine thrust error, VAC. 

gain on right engine thrust error, VAC. 

gain on left engine core dynamics. Second order filter. 

gain on right engine core dynamics. Second order filter. 

angle-of-attack sensor gain value. 

FCS logical used in g-limit on stick input, default=.T. (DO NOT CHANGE) 

desired load factor used in trim search when TCASE=3, g’s 

FCS logical used for g load override (TK5, TK8), default = .F. (DO NOT 

CHANGE) 

vertical distance from sea level to aircraft c.g., ft 
d/dt of H, ft/sec 

vertical distance from ground to the a.r.c., ft. 
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HIC 

HRF 


HSTORE 

I 

lALG 

( 1 ) 

( 2 ) 

(3) 

( 4 ) 

(5) 

( 6 ) 

(7) 

(9) 

IDYNCK 

( 0 ) 

( 1 ) 

( 2 ) 

(3) 

II 

IN2FT 

ISBNS 

( 0 ) 

( 1 ) 

ITFLG 

ITR 


ITRMX 

IXSEL 


IXX 

IXZ 

lYSEL 


lYY 

IZZ 

J 

KFS 

LAMDA 

LCAASP 


LCASYP 


initial condition for H, ft 

1962 standard atmosphere altitude in feet corresponding to a given an air density, 
p, in slugs/ft3 

FCS logical for heavy stores, default = .F. (DO NOT CHANGE) 
summing index used throughout simulation 
integer switch to select ACSL integration algorithm: 

Adam s-Moul ton, variable step, variable order 
Gear's stiff; variable step, variable order 
Runge-Kutta isi order or Euler 

Runge-Kutta 2iul order (default for this simulation) 

Runge-Kutta 4di order (default for ACSL) 

User supplied subroutine (INTEG) 

Runge-Kutta-Fehlberg 2Dd order 

Runge-Kutta-Fehlberg, 51h order 

Integer switch used to invoke dynamic check capability 

No disturbance input 

One sided pulse 

Doublet 

User defined signal, read from array SIG 
index used for sensor data tables 
conversion factor, 1/12 (ft/in) 

FCS flag for sensor model 

sensor effects off, ideal measurements 

MDC sensor model (default) 

integer switch used to control program flow when trimming 

loop count used in inner loop of trim search. When a new candidate solution vector 
is calculated, the resulting derivatives are then calculated ITRMX times. This 
allows implicit loops in the trim loop to settle out. Implicit loops ^ caused by the 
fact that accelerations are inputs the the control system that determine control 
surface position. An alternate way to break the loop would be to include actuator 
states in the trim vector. The trim search works best when the trim vector is kept 
small. 

number of passes made through equations of motion in inner trim loop. See ITR. 

Default = 6. 

integer pointer array of size NXTRMX used to select elements of trim solution 
vector, XTRIM. The first NXTR elements of IXSEL select from 29 possible 
independent variables. XTRIM is formed in subroutine XFORM. 
aircraft moment of inertia about x-body axis, slug-ft^ 
aircraft xz cross product of inertia, slug-ft^ 

integer pointer array of size NXTRMX used to select elements to drive to zero in 

the trim output vector, YTRIM. The first NYTR elements of lYSEL select from 15 

possible output variables. YTRIM is formed in subroutine YFORM. 

aircraft moment of inertia about y-body axis, slug-ft2 

aircraft moment of inertia about z-body axis, slug-ft2 

simulation summing index 

longitudinal stick force spring, = 7.0 Ibs/in, input to F18FCS 
difference between ground track angle and heading angle, radians 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric aileron 
command. 

(.F.) no action (default) 
logical switch 



LCFASP 

LCFSYP 

LCNASP 

LCNSYP 

LCPASP 

LCPSYP 

LCRASP 

LCRSYP 

LCSASP 

LCSBP 

LCSSYP 

LDBA 

LDBMI 

LDBOUT 


(.T.) A dynamic check perturbation is added to the symmetric aileron command. 
(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric trailing-edge 
flap command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the symmetric trailing-edge flap 
command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric leading-edge 
flap command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the symmetric leading-edge flap 
command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric p>ower lever 
angle command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the symmetric power lever angle 
command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric rudder 
command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the symmetric rudder command. 
(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the antisymmetric stabilator 
command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the sjjeed brake command. 

(.F.) no action (default) 
logical switch 

(.T.) A dynamic check perturbation is added to the symmetric stabilator 
command. 

(.F.) no action (default) 

logical switch - used for debugging aerodynamic calculations. 

(.T.) Extended debugging output is written to unit PRN for each call to the 
aerodynamic coefficient calculations. 

(.F.) no action (default) 

logical switch - used for inertial effects calculation debugging. 

(.T.) Extended debugging output is written to unit PRN for each calculation of 
the inertial components of the generalized force vector and mass matrix. 
(.F.) no action (default) 
logical switch - primary simulation debugging logical 
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LDBSS 

LDBTR 

LEULER 

LFCS 

LFTD 

LLDTR 

LUMDF 

LLNCSY 


LPCAP 

LPCRP 

LPCSP 

LQSE 


(.T.) Extended ACSL debugging output is written to unit PRN after the trim 
search or when time reaches TSTP. 

(.F.) no action (default) 

logical switch • debugging logical for the subroutine that calculates the accelerations 
by solving the vector equation [M](a} = (f) 

(.T.) Debugging output is written to unit PRN from subroutine XRBVDC. 

(.F.) no action (default) 

logical switch • used for debugging failed trim attempts. 

(.T.) Each new solution trial, XTRIM, and the resulting YTRIM vector are 
written to unit=6 (typically the screen) and also unit=PRN 
(.F.) no action (default) 
logic^ switch; 

(.T.) Euler angles rates are integrated directly. 

(.F.) quaternion derivatives are integrated, (default) 
logical switch - enables flight control system - read in from file "fl8tr.dat" 

(,T.) FCS on (default) 

(.F.) FCS off 

logical switch - enables use of flight data to drive control surface deflections 
(.T.) Control deflections set to data in FI DAT. Actuators bypassed. 

(.F.) Normal simulation operation, (default) 
logical switch 

(.T.) When START command is executed, initial conditions for all possible 

XTRIM variables are read in from file "fl8tr.dat". The default on the 
first execution of START is .T. Reset to .F. after "fl8tr.dat" 
is read in. 

(.F.) no action 
logical switch 

(.T.) Assign hardwired values to array of upper and lower bounds on elements of 
XTRIM variables. Trim algorithm confines search to this space. The 
hardwired values are defin^ in array DEFLIM in subroutine TRIM. 

(.F.) Use limit values read from "fl8tr.dat" file, (default) 
logical switch, used to control the 'LN' exogenous inputs that can be created to 
support the 'ANALYZ' command In ACSL. Dictates whether the exogenous inputs 
available to the ANALYZ option, which are used to calculate _LN1 inputs, are 
symmetric/antisymmetric quantities or right/left quantities. 

(.T.) _LN 1 = _S YLN + _ASLN (default) 

(.F.) _LN1=_LN 

where _LN 1 is summed with actual control surface position in degrees. Examples 
of _LN1 variables and DSRLNl (left stabilator), DSLLNl (right stabilator), and 
DRLLNl (left rudder), 
logical switch 

(.T.) Put dynamic check perturbation on lateral stick input. 

(.F.) no action (default) 
logical switch 

(.T.) Put dynamic check perturbation on rudder pedal input. 

(.F.) no action (default) 
logical switch 

( T.) A dynamic check perturbation is added to the longitudinal stick. 

(.F.) no action (default) 

logical switch - passed to subroutine SFAERRF and USAERRF - used to 
enable/disable the quasi-static -elastic (QSE) modeling option in the aerodynamic 
model 
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LRDFTI 


LRDWR 


LRTE 


LTHCGL 

LTHCGR 

LTHDMS 


LTHVEC 


LTR 


LTRFLG 


LVDEF 


LWRFTO 


LWRTR 


MACH 

MACHRF 

MACHTR 

MARR 

MASS 

MINT 

MIRR 

MMRR 

MNSPIN 


( T ) normal QSE aerodynamics. All flex/rigid ratios and increments apply, 
(default) 

(.F.) All flexMgid ratios and increments are set to limiting value as altitude 
becomes large or density goes to zero. 

logical switch - used to drive control surface positions according to a data file - as 
would occur with flight test data 

(,T ) Subroutine FTDRDF is called. Data is read from file="mhftdi.dat". 

( F ) no action (default) 

internal integer switch - intermediate flow control switch set in the simulation and 
used to control action in call to subroutine RDWRTR - can not be set by runtime 
commands 

(1) causes fl8tr.dat trim file to be read in (default) 

(2) causes f 1 8tr.dat trim file to be written on 

logical switch - passed to subroutine SFAERRF - used when conducting 
comparisons to recover the form of the aerodynamic model found in the original 
DMS simulation 

(.T.) activates form of subroutine SFAERRF that agrees with 
dmsfIS math model (default) 

(.F.) activates changes discussed in section 5.1 - aerodynamic model 

left engine thrust moment about body frame x-axis at c.g., ft-lbs 

right engine thrust moment about body frame x-axis at c.g., ft-lbs 

logical switch - used to duplicate the DMS throttle position 

(.T.) use DMS scaled throttle position 

(.F.) normal MDC FI 8 throttle position (default) 

logical switch - used to invoke preliminary thrust-vectored engine model 

(.T.) preliminary thrust-vectored engine mixlel (default) 

(.F.) non-thrust-vectored engine model 
logical switch 

(.T.) Trim search is initiated after runtime START command. If trim fails, TSTP 
set = 0. (default) 

(.F.) No trim attempted. Integration proceeds with current i.c.’s until T=TSTP. 
internal logical switch - used to indicate current TRIM status 
(.T.) simulation in trim search loop 
CF.) not in trim search loop 

logical flag - used to obtain extended output with variable definitions form 
dictionary 

(.T.) All simulation variable values printed out onto unit=PRN with dictionary 
definitions taken from unit=29. 

(.F.) no action (default) 

logical switch - used for FLT data comparison 

(.T.) subroutine FTDWRF called - data written to file="ftifo.dat" 

(.F.) no action (default) 

logical switch - used to write new file "fl8tr.dat" 

(.T.) Write the next calculated trim values to file "fl8tr.dat". Set LWRTR to .F. 
(.F.) no action (default) 

Mach number, n.d. 

Free-stream Mach number at a.r.c., n.d. 

Mach number to be achieved in trim search if "1" is in the active IXSEL list, n.d. 
aerodynamic component of rigid-rigid quadrant of mass matrix 
aircraft mass, slugs 

minimum integration step size (default value supplied by ACSL) 
inertial component of rigid-rigid quadrant of a mass matrix 
total rigid-rigid quadrant of a mass matrix 
FCS logical spin switch. Default = .F. (DO NOT CHANGE). 
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MTHCGL 

MTHCGR 

MUDGTR 

MUTRMX 


MXITTR 

MXSTP 

NE 

NSTP 


NTHCXjL 

NTHCGR 

NTVOPT 

NXTR 

NXTRMX 

NYLOC 

NYSSG 

NYTR 

NZSSG 

NZZR 

ORIDE 

P 

PIL 

PIR 

P2L 

P2R 

PAGSIZ 

PCA 

PCAFTD 

PCASG 

PCATR 

PCR 

PCRFTD 

PCRSG 

PCRTR 

prs 

PCSFTD 

PCSSG 

PCSTR 

PD 

PDG 

PGAIN 

PHI 

PHIO 

PHID 

PHIDG 


left engine thrust moment about body frame y-axis at c.g., 

right engine thrust moment about body franie y-axis at c.g., ft-lbs. 

bank angle in degrees sought in trim search if "29 included m active DCSEL list. 

Will produce coordinated turn if TCASE = 2. 

In Hm& flSbas trim strategy (section 2.2), M^TR^ is equivalent to v^bk 
"MUMAX" in ACSL TRIM algorithm (see section 5.2.1, ACSL REPbKbrMCH 

MANUAL REV. 4.1) - default = 10 ^ , u «n 

maximum number of iterations for simulation trimming scheme - default - Jto 

maximum integration step size - default = 0.00625 seconds 

square root of sum of squares of quaternion intermediate states, used to normalize 

quaternions at each time step to prevent numerical problems - n.d 

minimum number of integration steps in a communication interval. Detault - . 

This allows integration step size to be controlled by MXSTP as per 

recommendation in ACSL manual. f lu 

left engine thrust yaw moment about body z-axis at c.g. (TV), 

right engine thrust yaw moment about body z-axis at c.g. (YV), fLlbs 

integer switch - PCS flag for thrust vector pitch control blend (l=ON). Default - 


number of independent variables for simulation trim search (< NXTRM^ 
maximum number of independent variables for simulation trim search (10) 

3 element array describing location of ny accelerometer in FS/BL/WL coordinates, 


inches 

lateral acceleration sensor signal, g's ^ ktv'td\ 

number of dependent variables for simulation mm algorithm (-NA 1 Kj 
normal acceleration sensor signal, g's * i • 

normal acceleration 1-g steady state removed. NZZR - - 

PCS logical to override angle-of-attack limiter. Default = .F. (DONOl 

CHANGE). 


body axis roll rate, rad/sec 

gain in left core engine dynamics. order filter. 

gain in right core engine dynamics. 20^1 order filter 

gain in left core engine dynamics. 2^ order filter 

gain in right core engine dynamics. order filter 

Integer that defines page size as number of lines. (45 for NOS, 55 for 

pilot lateral stick position, inches 

pilot lateral stick position from flight test data, inches 

pilot lateral stick position signal array for IDYNCK = 3 

pilot lateral stick position required for trim, inches. 

pilot rudder pedal force, lbs 

pilot rudder pedal force from flight test data, lbs 

pilot rudder pedal force signal array for IDYNCK=3 

pilot rudder pedal force required for trim, lbs 

pilot longitudinal stick position, inches 

pilot longitudinal stick position flight test data, inches 

pilot longitudinal stick position signal array for IDYNCK=3 

pilot longitudinal stick po.sition required for trim, inches 

body axis roll acceleration, rad/sec^ 

body axis roll rate , degrees/sectmd 

gain on dynamic check perturbation command 

Euler roll angle, radians 

Euler roll angle intermediate state, integral of PHID 
Euler roll angle rate of change, rad/sec 
Euler roll angle, degrees 


VAX) 
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PHILN 

PHITR 

PHIZR 

PIC 

PLAL 

PLAR 

PLDL 

PLDR 

PRN 

PS 

PSO 

PSl 

PSAS 

PSI 

PSIO 

PSID 

PSIDG 

PSILN 

PSITR 

PSTATC 

PSTS 

PSTSO 

PSTSl 

PT 

PTl 

PTURB 

PW 

PWDG 

Q 

QBAR 

QBARS 

QBARSB 

QBARSC 

QC 

QCDA 

QCDB 

QCDC 

QCDD 

QCDE 

QCDF 

QCDG 

QCDH 

QCIS 

QCISO 

QCISl 

QCX 

(¥> 

QDG 

QIC 

QS 

QSO 


Euler roll angle linearization dummy variable, radians 

Euler roll angle required for trim, radians 

Euler roll angle referenced to "PHITR", radians 

body axis roll initial condition, rad/sec 

left engiiK total throttle position (power level angle), degrees 

right engine total throttle position (power level angle), degrees 

left throttle boost output error term. 

right throttle boost output error term 

high volume output goes to unit=PRN. Passed to Fortran subroutines that write to 
an output file. Default = 9. 

sensed roll rate , Vsec, = PSO (if ISENS=0) or = PSl (if ISENS=1) 

actual roll rate, Vsec 

output of rate sensor model, Vsec 

FCS logical for pitch SAS on. Default = .T. (DO NOT CHANGE). 

Euler yaw angle, radians 

Euler yaw angle intermediate stale, integral of PSID 
Euler yaw angle rate of change, ra^sec 
Euler yaw angle, degrees 

Euler yaw angle rate linearization dummy variable, radians. [Constant] 

Euler yaw angle required for trim, radians 

static atmospheric pressure, = f(H), PSF 

output of static atmospheric pressure sensor, PSF 

value assigned to PSTS if ISENS=() (Ideal sensors) 

value assigned to PSTS if ISENS=1 (Sensor model #1) 

intermediate variable in compressible pressure computation 

intermediate variable in compressible pressure computation 

perturbation value used for pulse or double signal if IDYNCK=1 OR 2 Default=0 

wind axis roll rate, rad/sec 

wind axis roll rate, Vsec 

body axis pitch rate, rad/sec 

aircraft dynamic pressure, q, Ibs/ft^ 
intermediate variable, =qS, lbs 
intermediate variable, = qSb, lbs 

intermediate variable, = qSc, lbs 

compressible impact pressure, PSF 

function used in compressible impact pressure computation 

function used in compressible impact pressure computation 

function used in compressible impact pressure computation 

function used in compressible impact pressure computation 

function used in compressible impact pressure computation 

function used in compressible impact pressure computation 

function used in compressible impact pre.ssure compu tat ion 

function used in compressible impact pressure computation 

output of compressible pressure sensor model, PSF 

compressible pressure-ideal sen.sor, PSF 

compressible pressure-sensor model #1, PSF 

intermediate value in sensed compressible pressure computation 

body axis pitch acceleration, rad/sec^ 

body axis pitch rate, '/sec 

body axis pitch rate initial condition, rad/sec 

sensed pitch rate, Vsec, = QSO or =QS 1 based on "ISENS" 

actual pitch rate, Vsec 
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QSl 

QW 

QWDG 

R 

RA2DG 

RD 

RDG 

RHO 

RHOSL 

RIC 

RRHO 

RS 

RSO 

RSI 

RSAS 

RTABL 

RTABR 

RTELLL 

RTELLR 

RTEULL 

RTEULR 

RW 

RWDG 

SBERR 

SBNLRL 

SETVDC 

SETVSC 

SIG 


SIGCPS 

SIGNAL 

SIGPCA 

SIGPCR 

SIGPCS 

SNALF 

SNALRF 

SNAZEL 

SNAZER 

SNBET 

SNPHI 

SNPH02 

SNPS02 

SNTHE 

SNTH02 

SNTVOA 

SPSFXD 

STORES 


output of rate gyro model, Vsec 

wind axis pitch rate, rad/sec 

wind axis pitch rate, ‘/see 

body axis yaw rate, ‘/see 

radians to degrees conversion factor, Vrad 

body axis yaw acceleration, rad/sec^ 

body axis yaw rate, Vsec 

air ^nsity, slugs/ft^ 

air density at sea level, 0.00237688 slugs/ft^ 

body axis yaw rate initial condition, rad/sec 

ratio of actual air density to RHOSL, n.d. 

sensed yaw rate,.‘/sec, = RSO OR =RSl based on "ISENS 

actual yaw rate, Vsec 

output of rate gyro model, Vsec 

FCS logical for SAS on. Default = .T. (DO NOT CHANGE), 
limited left afterburner thrust rate of change, Ibs/sec 
limited right afterburner thrust rate of change, Ibs/sec 
left engine core spool-up/-down rate lower limit, = - 5500 Ibs/sec 
right engine core spool-up/-down rate lower limit, = - 5500 Ibs/sec 
left engine core spool-up/-down rate upper limit, = 5500 Ibs/sec 
right engine core spool-up/-down rate upper limit, = 5500 Ibs/sec 
wind axis yaw rate, rad/sec 
wind axis yaw rate, Vsec 

error term between speed brake position and command, degrees 

speed brake no load rate limit, = 60/2.5, degrees/second 

set thrust vectoring differential command, input to subroutine F18FCS, 

default = -101. [Constant] moc/^o 

set thrust vectoring symmetric command, input to subroutine F18FCS, 

default = -101. [Constant] 

array of input values, in order of increasing time. Use SIG for open loop control 
surface deflections independent of of simulated pilot input 
SIG(1 ->100) = Signal amplitude values 

SIG(101->200) = Time breakpoint of corresponding amplitude values 

value of throttle signal at current time, = CSPYSG(T) 

value of generic signal at current time, = SIG(T ) 

value of lateral stick signal at current time, = PCASG (T) 

value of rudder pedal signal at current time, = PCRSG (T) 

value of longitudinal stick signal at current time, = PCSSG (T) 

sin( a = angle of attack at c.g.) 

sin( Oref = angle of attack at a.r.c.) 

sine of azimuth angle of installed left engine n.d. 

sine of azimuth angle of installed right engine n.d. 

sin( P = sideslip angle at c.g.) 


sin((J>) 

sin( PHITR / 2 ), n.d. 

sin( PSITR / 2 ), n.d. 

sin( THETA), n.d. 

sin( THETR / 2 ), n.d. 

sine of thrust-vectoring vane orientation 

sensed static pressure gain value 

FCS logical for stores logic (default = F. 


angle, n.d. 

) (DO NOT CHANGE). 
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SWRF 

T 

T1 

T2 

T3 

TABL 

TABUC 

TABR 

TABRIC 

TANTHE 

TAUGMN 

TAUGMX 

TAUL 

TAULS 

TAUR 

TAURS 

TCASE 


TEMPQ 

TESTVL 

TESTVR 

TH 

THE 

THEO 

THED 

THEDG 

THELN 

THETR 

THL 

THLO 


THLLN 

THR 

THRO 


THRCL 

THRCR 

THRFL 

THRFR 

THRLN 

THXL 

THXR 

THYL 

THYR 

THZL 

THZR 

TIDLE 


aerodynamic (wing) reference area, 400 
time, sec 

temporary variable used in rotational acceleration equations = IXX*P - IX7*R 
temporary variable used in rotational acceleration equations, = IYY*0 

temporal variable used in rotational acceleration equations, = - IXZ*P + I77*R 
left afterburner thrust, lbs "i » ^ r -r k 

left afterburner thrust initial condition, lbs 
right afterburner thrust, lbs 
right afterburner thrust initial condition, lbs 
tan(6) 

minimum engine thrust with afterburner engaged, lbs 
maximum engine thrust with afterburner engaged lbs 
left engine time constant for VAH first order lag dynamics 
sa>^ value of TAUL from last time in core engine dynamics 
right engine time constant for VAH first order lag dynamics 
saved value of TAUR from last time in core engine dynamic 
integer switch to determine trim solution sought 
(1) no action - used for normal wings-level trim 

digTOs', angle in 

nequired for symmetric puli-up at load factor in g’s selected by 

coefficient = QBAR / 100 

test value for left engine core dynamic nonlinear switch logic 
test value for right engine core dynamic nonlinear switch logic 

total engine thrust for thrust vectored engine lbs 

Euler pitch angle, radians 

Euler pitch angle intermediate state, integral of THED 
Euler pitch angle rate of change, rad/sec 
Euler pitch angle, degrees 

Euler pitch angle linearization dummy variable, radians 

Euler pitch angle rate required for trim, radians 

left engine net thrust, lbs 

bft engine net thrust intermediate state, lbs 

= TTL if in military power range 

= TABL if in afterburner 

left engine thrust linearization dummy variable lbs 

right engine net thrust, lbs 

right engine net thrust intermediate state, lbs 

= TTR if in military power range 

= TABR if in afterburner 

thrust command sent to left engine, lbs 

thrust command sent to right engine, lbs 

left engine reference thrust. < 0 if in afterburner dynamics lbs 

right engine reference thrust. < 0 if in afterburner dynamics lbs 

nght engine thrust linearization dummy variable, lbs 

left engine net thrust component in x-direction (NTV), lbs 

nght engine net thrust component in x-direction (NTV), lbs 

left engine net thrust component in y-direction (NTV), lbs 

nght engine net thrust component in y-direction (NTV), lbs 

left engine net thrust component in z-direction (NTV), lbs 

nght engine net thrust component in z-direction (NTV) lbs 

engine thrust at idle throttle setting, lbs 
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TMIL 

TNETL 

TNETR 

TRTOL 

TSFCS 

TSFTD 

TSTP 

TT 

TTL 

TTR 

TVDBOT 

TVLUC 

TVLUIC 

TVOA 

TVRUC 

TVRUIC 

U 

UA 

UA2 

UARF 

UD 

UDA 

UDWG 

UIC 

UWG 

V 

VA 

VA2 

VACL 

VACR 

VAEOL 

VAEOR 

VAEDL 

VAEDR 

VAEL 

VAELIC 

VAER 

VAERIC 

VAFL 

VAFLIC 

VAFR 

VAFRIC 

VAHL 

VAHLS 

VAHR 

VAHRS 

VAKIDL 

VAKIDR 

VAKDDL 


maximuin military thrust, lbs 

miscellaneous output - left engine net thrust (TV), lbs 

miscellaneous output - right engine net thmst (TV), lbs 

convergence requirement for successful trim - default = 0.00005 

time interval at which flight control system discrete is invoked, seconds 

time interval between flight test data points, sec 

simulation stop time - default = 0.0 seconds 

total thrust from left and right engine thrust, lbs 

sum of left engine thrust terms VAKL & VAFL (NTV), lbs 

sum of right engine thrust terms VAKL & VAFL (NTV), lbs 

FCS logical . Used to removed thrust vectoring deadband. Makes trimming 

possible when thrust vectoring is selected. Default = .T. 

CTVLL initial condition 

CTVLU initial condition 

thrust -vectoring vane orientation angle, degrees 

CTVRL initial condition 

CTVRU initial condition 

body axis inertial velocity in x-btxly direction, ft/sec 
body axis atmospheric velocity in x-body direction, ft/sec 

UA squared, ft^/sec^ 

UA at a.r.c., ft/sec 

body axis inertial acceleration in x-direction, ft/sec^ 
d/dtofUA,ft/sec2 

d/dt of aerodynamic disturbance component in x-body axis, ft/sec^ 

body axis inertial velocity in x-dircction i.c., ft/sec 

x-body component of aenxlynamic disturbance velocity, ft/sec 

body axis inertial velocity in y-body direction, ft/sec 

body axis atmospheric velocity in y-body direction, ft/sec 

VA squared, ft^^ec^ 

left engine error term. = THRCL - VAFL 

right engine error term. = THRCL - VAFL 

rate-limited left engine spool- up-/-down intermediate state, Ibs/sec 

rate-limited right engine spool-up-/-down intermediate state, Ibs/sec 

artificial fast derivative inserted to achieve explicit formulation, the integral of which 

isVAEOL 

artificial fast derivative inserted to achieve explicit formulation, the integral of which 
is VAEOR 

rate-limited left engine spool-up-/-down , Ibs/sec 

rate-limited left engine spool-up-/-down initial condition, Ibs/sec 

rate-limited right engine spool-up-/-down , Ibs/sec 

rate-limited right engine spool-up-/-down initial condition, Ibs/sec 

main component of left engine core dynamics, lbs 

main component of left engine core dynamics initial condition, lbs 

main component of right engine core dynamics, lbs 

main component of right engine core dynamics initial condition, lbs 

output of left engine core dynamics nonlinear switch, lbs 

saved value of VAHL from last time in core engine dynamic 

output of right engine core dynamics nonlinear switch, lbs 

saved value of VAHR ftom last time in core engine dynamic 

intermediate left core engine dynamic, saved past value, Ibs/sec 

intermediate right core engine dynamic, saved past value, Ibs/sec 

left core engine spool up input, Ibs/sec^ 
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VAKDDR 

VAKDL 

VAKDR 

VAKL 

VAKUC 

VAKR 

VAKRIC 

VARF 

VD 

VDA 

VDBA 

VDWG 

VGS 

VHIDL 

VHIDR 

VHIL 

VHILIC 

VHILS 

VHIR 

VHIRIC 

VHIRS 

VIC 

VKDLIC 

VKDRIC 

VS 

VT 

VTD 

vnNv 

VTRF 

VTRINV 

VTTR 

VWG 

W 

WA 

WA2 

WARF 

WD 

WDA 

WDWG 

WFL 

WFR 

WIC 

WIDTH 

WLCG 

WLENG 

WLRF 

WT 

WWG 

XALF 

XARF 


right core engine spool up input, Ibs/sec^ 
intermediate left core engine dynamic, Ibs/sec 
intermediate right core engine dynamic, Ibs/sec 
output of left engine core 20d order dynamics, lbs 
initial condition of integrator for VAI^ 
output of right engine core 2iKl order dynamics, lbs 
initial condition of integrator for VAKR 
VA at a.r.c., ft/sec 

body axis inertial acceleration in y-direction , ft/sec^ 
d/dt of VA, ft/scc^ 

deadband applied to thrust vectoring vane deflection, radians 

d/dt of aerodynamic disturbance component in y-body axis, ft/sec^ 
aircraft ground speed, ft/sec 
VHIL rate of change, Ibs/sec 
VHIR rate of change, Ibs/sec 

left engine state created by ]S1 order lag mechanization, lbs 

left engine core thrust at initial condition, lbs 

saved value of VHIL' from last time in core engine dynamic, lbs 

right engine state created by 1^1 order lag mechanization, lbs 

right engine core thrust at initial condition, lbs 

saved value of 'VHIR' from last time in core engine dynamic., lbs 

body axis inertial velocity in y-direction i.c., ft/sec 

initial condition on first integral in left core engine dynamic, Ibs/sec^ 

initial condition on first integral in right core engine dynamic, Ibs/sec^ 

speed of sound, fi/sec 

aircraft total atmosphere velocity, ft/sec 

aircraft total atmosphere acceleration, ft/sec^ 

coefficient to normalize by total velocity, sec/ft 

total airspeed at a.r.c., ft/sec 

coefficient to normalize by total velocity at a.r.c., sec/ft 
aircraft total atmosphere velocity after successful trim, ft/sec 
y-body component of aerodynamic disturbance velocity, ft/sec 
body axis inertial velocity in z-body direction, ft/sec 
body axis atmo^heric velocity in z-body direction, ft/sec 
WA squared, ft^/sec^ 

WA at a.r.c., ft/sec 

body axis inertial acceleration in z-direction, ft/sec2 
(d/dt)ofWA, ft/sec2 

d/dt of aerodynamic disturbance component in z-body axis, ft/sec^ 
left engine fuel flow 
right engine fuel flow 

l^y axis inertial velocity in z-direction i.c., ft/sec 

time span to hold step input dynamic check perturbation command,.default = 1 0 
seconds 

water line coordinate of the center of gravity, inches 

water line coordinate of the engine/airframe interface, right/left engine, 100.0 in. 
water line coordinate of the aerodynamic reference center, = 100.0 inches 
aircraft weight, lbs 

z-body component of aerodynamic disturbance velocity, ft/sec 
air data probe alpha with source noise added, degrees 
x-direction distance from c.g. to a.r.c., ft 
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XBAR noimalized x-location of accelerometer from c.g., n.d. 

XD aircraft inertial ground track velocity in track direction , ft/sec 

XDDWG d/dt of XDWG, ft/scc^ j- 

XDWG aerodynamic disturbance velocity component in x-earth frame direction, fVsK 

XDWGIC aerodynamic disturbance velocity component in x-earth frame direction - imtial 

condition, ft/sec 

XHMl sensed static pressure breakpoint vector . 

XRFCG ’x’ component of the location of right (left) engine/airfirame interface point in c.g.- 
centered^body frame, = -(FSENG - FSCG)/12„ ft 
XTHCG total thrust force at c.g. in body x-direction (TV), Iw 

XTHCXjL left engine thrust force at c.g. in body x-direction Ids 

XTHCGR right engine thrust force at c.g. in body i 

XTRLM 2 by 30 array with upper and lower limits on XTRIM independent variables 

XTVL left engine thrust force in body x-direction (TV), lbs 

XTVR right engine thrust force in body x-direction (TV), lbs 

YARF y-dircetion distance from c.g. to a.r.c., ft 

YBAR normalized y-location of accelerometer from c.g., n.d. 

YD aircraft inertial ground track velocity in crosstrack direction, it/sec 

YDDWG d/dt of YDWG, ft/sec2 . 

YDWG aerodynamic disturbance velocity component along earth frame y-axis, tys^ 

YDWGIC aerodynamic disturbance vekKity component along earth frame y-axis - initiaJ 

condition, ft/sec 

YGKA02 sensed angle-of-attack gain table 

YGKAOA sensed angle-of-attack gain table . , . , . . 

YRFCX5L V component of the location of left engme/airframe interface point in c.g.- 
centered body frame, = (BLENGL - BLCG)/12., ft. 

YRFCGR V component of the location of right engine/airframe interface point in c.g.- 
centered body frame, = (BLENGR - BLCG)/1^, ft. 

YSAS FCS logical for yaw SAS on. Default = .T. (E^ NOT CHANGE) 

YTHCXjL left engine thrust force at c.g. in body y-di^tion (TV), lbs 

YTHCXJR right engine thrust force at c.g. in btxly y-direction (TV), lbs 
YTVL left engine thrust force at nozzle in y-body direction CTV). 'Ds 

YTVR right engine thrust force at nozzle in y-body direction (TV), lbs 

ZARF z-direction distance from c.g. to a.r.c., ft 

ZBAR normalized z-l(K:ation of accelerometer from c.g., n.d. 

ZDDWG (d/dt)ofZDWG,ft/sec2 , . r • 

ZDWG aerodynamic disturbance velocity comjxinent along earth frame z-axis, tys^ 

ZDWGIC aerodynamic disturbance velocity component along earth frame z-axis - iniual 

condition, ft/sec . , . _r • . „ „ 

ZRFCG ’z’ component of the location of right (left) engine/airframe interface point in c.g. 

centered body frame, = -(WLENG - WLCG)/12., ft. 

ZTHCX3L left engine thrust force at c.g. in body z-direction IDs 
ZTHCGR right engine thrust force at c.g. in l^y z-direction (TV), lbs 
ZTVL left engine thrust force in body z-direction (TV), lbs 

27TVR right engine thrust force in body z-direction (TV), lbs 

ZZABL d/dt of left engine thrust in AB before rate limiting, Ibs/sec 
77.A RR d/dt of right engine thrust in AB before rate limiting, Ibs/sec 
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APPENDIX B - EXAMPLE CASES 
Trim to Level Flight (TCASE=1) 

The following wings level (TCASE=1) trim cases are demonstrated by an example terminal 
session: 

(1) M=.6, H=10,000 (ft),flight path angle (y) = 0. 

(2) M=.8, H=10,000 (ft),flight path angle (y) = 0. 

(3) a = 5*. H=10,000 (ft), flight path angle (y) = 0. 

(4) a = 10*. H=10,000 (ft),flight path angle (y) = 0. 

Commands entered at the keyboard are in bold courier. The program and VMS responses are 
in plain Courier. Comments inserted after the session to explain actions are in italics. 

AGCB9$ 8aetup_bas A.ssign data file.t to names read by program. 

$ assign [CSB.F18BAS .DATA_FILES]T5AER0.F18BAS tSaero 

$ assign [CSB.F18BAS .DATA_FILES]T5ENG.F18BAS tSeng 

$ assign [CSB .F18BAS .DATA_FILES] F18BAS .DCt for029 

$ set noverify 

AGCB9S type fl8tr.dat Show copy of f18tr.dat.dat file used to store 

trim information. 

F18 WITH 60% FUEL, TIP MISSLES, LEVEL TURN 


MASS, INERTIAS, AND CG LOCATION 


33310.0 

WT 

(LBS) 

WEIGHT 


23000.0 

I XX 

(SLUG*FT**2) 

INERTIA 

ABOUT X AXIS 

151293.0 

lYY 

(SLUG*FT**2) 

INERTIA . 

ABOUT Y AXIS 

169945.0 

IZZ 

(SLUG*FT**2) 

INERTIA . 

ABOUT Z AXIS 

“2971.0 

IXZ 

(SLUG*FT**2) 

XZ PLANE 

INERTIA PRODUCT 

455.0 

FSCG 

(IN) 

CG LOG - 

FUSELAGE STATION 

0.0 

BLCG 

(IN) 

CG LOC - 

BUTTOCK LINE 

102.8 

WLCG 

(IN) 

CG LOC - 

WATER LINE 

— 


FLIGHT CONDITIONS 

AND MODELING OPTIONS 

F 

LTHVEC 

1 THRUST VECTORING ON 


T 

LFCS 

1 FLT CONTROL 

SYTEM ON 


F 

LTHDMS 

1 USE DMS PLA 

SCHEDULE 


T 

LRTE 

1 USE R/T EQV 

AERO 


0.60000 

MACHTR 

(N.D.) 



10000.00 

HIC 

(FT) 



30.00 

MUDGTR 

(DEG) 



1.00000 

GLOAD 

(G) 






O’D TM 

SETUP IF LTR = .TRUE 





1 

TCASE 1 

1=LEVEL, 2=COORD TURN 

(MUDGTR), 3=PULL UP (GLOAD) 

4 

NUMBER 

OF DRIVER VARIABLES 


ELEMENT 

LIMITS 



VARIABLE NAME, UNITS 

14 

“2.500 

5.000 


PCSTR (INCHES ) 

16 

31.000 

127.000 


DPSYTR (% POWER) 

3 

-0.100 

1.500 


ALFTR (RADIANS) 

7 

-1,500 

1.500 


THETR (RADIANS) 


4 NUMBER OF DRIVEN VARIABLES 

ELEMENT VARIABLE NAME, UNITS 

1 UD (FT/SEC2) 
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3 

5 

7 


WD (FT/SEC2) 
QD (RAD/SEC2) 
GAMZR (RADIANS) 


INITIAL CONDITIONS 
0.60000002 
O.OOOOOOOOE+00 
0.43141015E-01 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
0.43145921E-01 
0 .OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
0.23838108E-01 
O.OOOOOOOOE+00 
60.480717 
O.OOOOOOOOE+00 
0.35326593E-01 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
O.OOOOOOOOE+00 
0 .OOOOOOOOE+OO 
3.2825463 
0 .OOOOOOOOE+OO 
3.4605167 
0. OOOOOOOOE+OO 
0 . OOOOOOOOE+OO 
30.000000 
0 . OOOOOOOOE+OO 


MACHTR 

(N.D.) 

BETTR 

(RADIANS) 

ALFTR 

(RADIANS) 

PIC 

(RAD/SEC) 

QIC 

(RAD/SEC) 

RIC 

(RAD/SEC) 

THETR 

(RADIANS) 

PHITR 

(RADIANS) 

PSITR 

(RADIANS) 

GAMTR 

(RADIANS) 

DTVL 

(DEGREES) 

DTVR 

(DEGREES) 

PCATR 

(INCHES ) 

PCSTR 

(INCHES ) 

PCRTR 

(LBS ) 

DPSYTR 

(% POWER) 

DPASTR 

(% POWER) 

DSSYTR 

(DEGREES) 

DSASTR 

(DEGREES) 

DASYTR 

(DEGREES) 

DAASTR 

(DEGREES) 

DRSYTR 

(DEGREES) 

DRASTR 

(DEGREES) 

DNSYTR 

(DEGREES) 

DNASTR 

(DEGREES) 

DFSYTR 

(DEGREES) 

DFASTR 

(DEGREES) 

CSB 

(DEGREES) 

MUDGTR 

(DEGREES) 

future 

expansion 


SEND OF DATA READING SECTION 
CASE SELECTIONS : 

TCASE 1 STRAIGHT & LEVEL STEADY STATE 

2 COORDINATED TURN STEADY STATE 

3 PULL-UP STEADY STATE 

TRIM VALUE SELECTIONS: 


TRIM 

DRIVER 

ARRAY 

TRIM OUTPUT ARRAY 

1 

MACHTR 

(N.D.) 

1 

UD 

(FT/SEC2) 

2 

BETTR 

(RADIANS) 

2 

VD 

(FT/SEC2) 

3 

ALFTR 

(RADIANS) 

3 

WD 

(FT/SEC2) 

4 

PIC 

(RAD/SEC) 

4 

PD 

(RAD/SEC2) 

5 

QIC 

(RAD /SEC) 

5 

QD 

(RAD/SEC2) 

6 

RIC 

(RAD/SEC) 

6 

RD 

(RAD/SEC2) 

7 

THETR 

(RADIANS) 

7 

GAMZR 

(RADIANS) 

8 

PHITR 

(RADIANS) 

8 

PHIZR 

(RADIANS) 

9 

PSITR 

(FIADIANS) 

9 

THE 

(RADIANS) 

10 

GAMTR 

(RADIANS) 

10 

LAMDA 

(RADIANS) 

11 

DTVL 

(DEGREES) 

11 

FYTOT 

(G S) 

12 

DTVR 

(DEGREES) 

12 

CL-CLTR 

(N.D. ) 
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13 

PCATR 

(INCHES ) 

14 

PCSTR 

(INCHES ) 

15 

PCRTR 

(LBS ) 

16 

DPSYTR 

(% POWER) 

17 

DPASTR 

(% POWER) 

18 

DSSYTR 

(DEGREES) 

19 

DSASTR 

(DEGREES) 

20 

DASYTR 

(DEGREES) 

21 

DAASTR 

(DEGREES) 

22 

DRSYTR 

(DEGREES) 

23 

DRASTR 

(DEGREES) 

24 

DNSYTR 

(DEGREES) 

25 

DNASTR 

(DEGREES) 

26 

DFSYTR 

(DEGREES) 

27 

DFASTR 

(DEGREES) 

28 

CSB 

(DEGREES) 

29 

MUDGTR 

(DEGREES) 

30 

future 

expansion 


13 CM 


(N.D.) 


AGCB9$ ac8l/run flBbas 


Execute file 'fl8has.exe'. 


The ACSL Licence will terminate in 3 months 

Ask your systems manager to contact Mitchell and Gauthier for new key 
Telephone (USA) 508/369-5115 

Switching CMD unit to 4 to read F18BAS.CMD 
ACSL>st art Stan case (I) TSTP =0.0 .w program slops 

immediately (rfter trim attempt. LI R true 
causes trim attempt to be made. 

»> READING VALUES TO INIT TRIM SEARCH (F18TR FILE) 

»> READING AERO DATA 

TIME ELAPSED = 7.238281 SECONDS 

»> READING ENGINE DATA 

TIME ELAPSED = 1.339844 SECONDS 

NORM OF Y ON FIRST PASS = 3 . 0825590E-03 

COLUMN 2 NOT REPEATABLE. ERROR = 6 . 519258022308350E-00 9 

TRIM ACHIEVED. ERROR = 4 . 4 3152 34E-0 6 

ITERATIONS REQUIRED = 27 

ACS^^iT t c^^rn=8 0 Set display column width=80 (reqiured only 

for this demo). 

ACSL>d machtr , hie, vt , ud, vd, wd, pd, qd, rd, thedg, alf dg, . . . 
ACSL>pcstr,dpsytr,phi 

HIC 10000.0000 VT 646.442000 

VD 0. WD 1.3487E-04 

QD-1.3995E-05 RD 0. 

ALB'DG 2.47177000 PCSTR-0 . 0687 6620 

PHI 0. 


MACHTR 0 . 60000000 
UD-2 . 3931E-04 
PD 0. 

THEDG 2.47138000 
DPSYTR 60.4787000 
ACSL>d ixsel 

IXSEL 14 

7 
0 
0 

ACSL>s niachtr=.8 
ACSL>8tart 

NORM OF Y ON FIRST PASS = 
TRIM ACHIEVED. ERROR = 
ITERATIONS REQUIRED = 
t exceeds tstp 


Display .selected XTRIM variables. 

16 3 

0 0 

0 0 

Change desired Mach for new trim. 
Start case (2) 


0.2624109 
1 . 0789349E-05 
23 



ACSL>d machtr, hie, vt , ud, vd, wd, pd, qd, rd, thedg, alfdg, . . . 
ACSL>pcstz, dpsytz 


MACHTR 0.80000000 
UD-2 .7995E-04 
PD 0. 

THEDG 1.32924000 
DPSYTR 80.5041000 
ACSL>8 ixsel=14, 16, 1, 7 


ACSL>a alftr=. 087266462 


HIC 10000.0000 
VD 0. 

QD 3.8318E-05 
ALFDG 1.32953000 


VT 861.923000 
WD 5.5112E-04 
RD 0. 

PCSTR-1 .00055000 


Change XTRIM select list by replacing 
ALFFR with MACHTR. This causes 
MACHTR to be varied while ALFIR is left 
fixed. Useful for trimming to a desired angle 
of attack. 

Set ALFTR in radians to produce 5 degrees 
a. 


ACSL>8t:art Start case (3) 

NORM OF Y ON FIRST PASS = 1.022277 

COLUMN 2 NOT REPEATABLE. ERROR = 4 . 951 1334054 84 1 43E-00 7 

COLUMN 4 NOT REPEATABLE. ERROR = 2 . 402740051365981E-007 

TRIM ACHIEVED. ERROR = 2 . 0998859E-06 

ITERATIONS REQUIRED = 16 

t exceeds tstp 

ACSL>d machtr, hie, vt , ud, vd, wd, pd, qd, rd, thedg, alfdg, . . . 
ACSL>pcstr , dpsytr, phi 


MACHTR 0.42064100 
UD-3. 4658E-05 
PD 0 . 

THEDG 5.00025000 
DPSYTR 55.0637000 
ACSL>8 alftr=. 1745329 


HIC 10000.0000 
VD 0. 

QD 7.1250E-06 
ALFDG 5.00000000 
PHI 0. 


VT 453.200000 
WD 7.4657E-06 
RD 0 . 

PCSTR 0.21334900 


Set ALFTR in radians to produce 10 degrees 


a. 


ACSL>start Start case (4) 

NORM OF Y ON FIRST PASS = 0.2984953 

TRIM ACHIEVED. ERROR = 2 . 27 4 92 62E-05 

ITERATIONS REQUIRED = 11 


t exceeds tstp 


ACSL>d machtr, hie, vt, ud, vd, wd, pd, qd, rd, thedg, alfdg, . . . 
ACSL>pc8tr, dpsytr, phi 


MACHTR 0.29906200 
UD-5.2243E-04 
PD 0. 

THEDG 10.0004000 
DPSYTR 60.9197000 


HIC 10000.0000 
VD 0 . 

QD 8.4760E-05 
ALFDG 10.0000000 
PHI 0. 


VT 322.210000 
WD 8.9943E-04 
RD 0 . 

PCSTR 0.50991900 


ACSL>Stop 


End ACSL run 
AGCB9$ 


End session 



Trim to Bank Angle (TCASE=2) 

The following trim to bank angle (TCASE=1) trim cases are demonstrated by ari example terminal 
session. The session is begun by reading in the TCASE=1 trim data file fl8tr.dat, then using run 
time commands to set up the bank angle trim. 

(1) M=.6, H= 10,000 (ft), wings level, TCASE=1 

(2) M=.6, H=10,000 (ft), TCASE=2, Bank angle (n) = 30* 

(3) M=.6, H= 10,000 (ft), TCASE=2, Bank angle (n) = 60* 


AGCB9$ acsl/run flSbas 


The ACSL Licence will terminate in 3 months 

Ask your systems manager to contact Mitchell and Gauthier for new key 
Telephone (USA) 508/369-5115 

Switching CMD unit to 4 to read F18BAS.CMD 
ACSL>start Start case (1) 

»> READING VALUES TO INIT TRIM SEARCH (F18TR FILE) 

»> READING AERO DATA 


TIME ELAPSED = 

»> READING ENGINE DATA 
TIME ELAPSED = 
NORM OF Y ON FIRST PASS = 
COLUMN 

TRIM ACHIEVED 


7.230469 


SECONDS 


1.339844 SECONDS 
3.0825590E-03 

2 NOT REPEATABLE. ERROR = 6 . 519258 022308 350E-0 0 9 

ERROR = 4 . 4315234E-06 

27 


ITERATIONS REQUIRED = 
t exceeds tstp 
ACSL>s tcwprn=80 

ACSL>d machtr, hie, tease, alfdg, thedg, pcstr, dpsytr 


MACHTR 0.60000000 
ALFDG 2.47177000 
DPSYTR 60.4787000 
ACSL>d ud, vd, wd, pd, qd, rd 
UD-2 .3931E-04 
PD 0 . 

ACSL>8 tcase=2 , llimd£= . t . , nxtr=8 


HIC 10000.0000 
THEDG 2.47138000 


VD 0 . 

QD-1 .3995E-05 


TCASE 1 

PCSTR-0 .06876620 


WD 1.3487E-04 
RD 0 . 


Convert from wings level to a trim at bank 
angle using the TCASE=2 option. Expand 
trim from 4x4 to 8x8 . Nlh'R will be 
automatically set to NXTR. 


Choose desired bank angle. 
Start case (2) 


ACSL>8 ixsel=15, 13, 2, 8, 7, 3, 14, 16 
ACSL>8 iysel=6, 4, 2, 11, 1, 3, 5, 7 

ACSL>8 nmdgtr=30 
ACSL>8tart 

NORM OF Y ON FIRST PASS = 7 . 2 1 31 425E-02 

MUMAX REACHED. CALCULATE NEW SENSITIVITY MATRIX. 

NORM OF Y ON FIRST PASS = 3 . 7534259E-02 

TRIM ACHIEVED. ERROR = 2 . 8 61 61 43E-05 

ITERATIONS REQUIRED = 46 

t tstp 

ACSL>d machtr, hie, phidg, mudgtr, az , alfdg, thedg, phidg, dpsytr 

HIC 10000.0000 
AZ-1 .15330000 
PHIDG 30.0299000 


MACHTR 0.60000000 
MUDGTR 30.0000000 
THEDG 2.42915000 
ACSL>d ud, vd, wd, pd, qd, rd 

UD 2.3536E-06 
PD-2 .2889E-04 


PHIDG 30.0299000 
ALFDG 2.81582000 
DPSYTR 61.2118000 


VD 6.9801E-05 
QD-1.9547E-06 


WD-1.3113E-05 
RD 1.6206E-06 
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Slari case (3) 


ACSL>* mudgtrs60 
ACSL>atart 

NORM OF Y ON FIRST PASS = 0.1437438 

COLUMN 8 NOT REPEATABLE. ERROR = 1 . 303851604461670E-008 

TRIM ACHIEVED. ERROR - 1 . 1091912E-05 

ITERATIONS REQUIRED -= 26 


t exceeds tstp 

ACSL>d machtr, hie, phidg, mudgtr, az, alfdg, thedg, phidg, dpsytr 

MACHTR 0.60000000 HIC 10000.0000 PHIDG 60.0811000 

MUDGTR 60.0000000 AZ-1 . 99337000 ALFDG 4.67035000 

THEDG 2.30884000 PHIDG 60.0811000 DPSYTR 76.6398000 

2 gsof normal load (stability frame) are produced by 60* of bank angle, AZ is in the 
body frame and is slightly different from NZ in the stability frame. 


ACSL>s lwrtr«.t 
ACSL>start 

NORM OF y ON FIRST PASS 
ALREADY TRIMMED. 


1.1091912E-05 


Set flag so next succesfful trim results are 
saved to flStrdat. 

Repeat case (3) to cause trim data to be 
written to file fI8tr,dat\ 


t exceeds tstp 
ACSL>atop 
End ACSL run 


End session 



Steady State Pull-Up (TCASE=3) 

The following steady-state pull-up (TCASE=3) trim cases arc demonstrated by an example terminal 
session. The session is begun by reading in the TCASE=1 trim data file fl 8 tr.dat, then using run 
time commands to set up the pull-up trim. 

(1) M=. 6 , H=10,000 (ft), wings level, TCASE=1, flight path angle ( 7 ) = 0. 

(2) M=. 6 , H= 10,000 (ft), TCASE=3, GLOAD (n^) = 2 g’s. 

(3) M=. 6 , H=10,000 (ft), TCASE=3, GLOAD (n^) = 4 g’s. 

AGCB9$ copy fl8tr_tcl.dat fl8tr.dat Copy fl8tr.dat file used for TCASE=1 trim 

cases into generic file read by flSbas 
simulation. 

AGCB9$ acsl/run flSbas 


The ACSL Licence will terminate in 3 months 

Ask your systems manager to contact Mitchell and Gauthier for new key 
Telephone (USA) 508/369-5115 


Switching CMD unit to 4 to read F18BAS.CMD 
ACSL>atart Start case (1) 

»> READING VALUES TO INIT TRIM SEARCH (F18TR FILE) 


»> READING AERO DATA 

TIME ELAPSED = 7.351563 SECONDS 

»> READING ENGINE DATA 

TIME ELAPSED = 1.492188 SECONDS 

NORM OF Y ON FIRST PASS = 3 . 0825590E-03 

COLUMN 2 NOT REPEATABLE. ERROR 6 . 51 9258022308350E-00 9 

TRIM ACHIEVED. ERROR = 4 . 4315234E-06 

ITERATIONS REQUIRED = 27 


t exceeds tstp 


ACSL>8 tcwprn=80 

ACSL>d tease, alfdg, thedg, machtr , hie 

TCASE 1 ALFDG 2.47177000 THEDG 2.47138000 

MACHTR 0.60000000 HIC 10000.0000 

AC S L >d ixsel Selected XTRIM variables 


IXSEL 14 

7 
0 
0 


16 

0 

0 


3 

0 

0 


ACSL>d iysel 
lYSEL 


1 

7 

0 

0 


Selected YTRIM variables 
3 
0 
0 


5 

0 

0 


ACSL>d az 

AZ-0. 99906600 

ACSL>s tease=3, gload=2 TCASE=^3 sets pitch rate automatically for a 

pull-up with load (g\s) = GLOAD. 

ACSL>start Start case (2) 

NORM OF Y ON FIRST PASS = 0.2461748 

COLUMN 2 NOT REPEATABLE. ERROR = 1 . 220281 404 8 1234 OE-0 07 

TRIM ACHIEVED. ERROR = 7 . 6116430E-06 

ITERATIONS REQUIRED = 10 


t exceeds tstp 

ACSL>d tease, al£dg,thedg,machtr, hie, az 

TCASE 3 ALFDG 4.67142000 

MACHTR 0.60000000 HIC 10000.0000 


THEDG 4.67179000 
AZ-1. 99336000 
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WD-2 . 6588E-04 


QD-2 .2835E-05 


ACSL>d ud, wd,qd 

UD 5.5424E-04 
ACSL>« gload»4 
ACSL>st«rt Start case (3) 

NORM OF Y ON FIRST PASS = 0.4928297 

COLUMN 3 NOT REPEATABLE. ERROR = 3 . 632492280303268E-007 

TRIM ACHIEVED. ERROR - 5 . 9509330E-07 

ITERATIONS REQUIRED = 20 

t exceeds tstp 

ACSL>d dpaytr,q Show required throttle (%) and pilch rate 

(ris). 

DPSYTR 112.504000 Q 0.14931300 

ACSL>d tease, alfdg, thedg, machtr, hie, az 

TCASE 3 ALFDG 8.99947000 THEDG 8.99940000 

MACHTR 0.60000000 HIC 10000.0000 AZ-3 . 9507 6000 

ACSL>StOp 
End ACSL run 
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